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Abstract 

We study three dimensional conformal field theories described by U (N) Chern- 
Simons theory at level k coupled to massless fermions in the fundamental repre- 
sentation. By solving a Schwinger-Dyson equation in lightcone gauge, we com- 
pute the exact planar free energy of the theory at finite temperature on as a 
function of the 't Hooft coupling A = N/k. Employing a dimensional reduction 
regularization scheme, we find that the free energy vanishes at |A| = 1; the con- 
formal theory does not exist for |A| > 1. We analyze the operator spectrum via 
the anomalous conservation relation for higher spin currents, and in particular 
show that the higher spin currents do not develop anomalous dimensions at lead- 
ing order in 1/N. We present an integral equation whose solution in principle 
determines all correlators of these currents at leading order in and present 
explicit perturbative results for all three point functions up to two loops. We also 
discuss a lightcone Hamiltonian formulation of this theory where a Woo algebra 
arises. The maximally supersymmetric version of our theory is ABJ model with 
one gauge group taken to be ^7(1), demonstrating that a pure higher spin gauge 
theory arises as a limit of string theory. 
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1 Introduction 

Chern- Simons theories are fascinating from several points of view. The attractive 
features of these theories, unique to three-dimensions, had been recognized early on, 
notably in [l] and [2]. They arise, for example, in the study of knot invariants js], the 
classification of two-dimensional rational conformal field theories [i] , and in the study of 
the quantum Hall effect [5]. More recently, super conformal Chern-Simons theories [6}[7] 
have been shown to play a crucial role in the AdS/CFT correspondence [s]. 

In this paper we will be interested in non-supersymmetric, but conformally in- 
variant, U{N) Chern-Simons theories with matter. As we shall explain below, these 
theories have an important distinguishing feature; they effectively admit large N non- 
supersymmetric lines of fixed points of the renormalization group. 

Lines of fixed points parameterized by a coupling constant are especially interesting 
from the viewpoint of the AdS/CFT correspondence. Such lines of d dimensional CFTs 
have the potential of interpolating between a simple field theoretic description at weak 
coupling and a relatively simple bulk gravitational description at strong coupling, as 
demonstrated by the famous supersymmetric examples [9|[8]. 

Now when d > 4, lines of fixed points appear to be rather exotic. The examples 
we know of (like the large Banks-Zaks fixed line of QCD) involve theories with a 
parametrically large number of flavours of matter fields. It is interesting, on the other 
hand, that effective fixed lines of large A^ Chern-Simons theories coupled to matter 
fields are plentiful and very easily constructed in 2+1 dimensions even with very simple 
matter contentQ 

Consider a level k, U{N) Chern-Simons theory coupled to a single fermion in any 
representation of the gauge group. The only gauge invariant power counting relevant 
or marginal operators in such a theory are the fermion kinetic term and mass term. 
A continuum quantum theory built from such a Lagrangian depends on two discrete 
parameters, k and A^, and a single continuous parameter, the physical mass m of the 
fermionic field. At energies E ^ m, the dynamics of this theory is scale invariant 
as well as nontrivial. Nontriviality is ensured by the fact that the discrete Chern- 
Simons coupling, which induces interactions among the fermions, cannot run and so 
is nonvanishing even at arbitrarily high energy. The nontrivial CFT that controls the 

^We thank O. Aharony for discussions on this point. 
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high energy behaviour of this system is most directly constructed by choosing the bare 
mass to set the physical mass m of this system to zero. 

The parameters k and labeling the CFT are discrete. However it is well known 
that the loop counting parameter in a U{N) Chern-Simons theory is the 't Hooft cou- 
pling A = ^ whenever all matter fields transform in representations whose dimension 
does not grow faster with N than A^^. In the large (and simultaneously large k) 
limit A is effectively a continuous parameter (exactly as in ABJM theory [s]). For this 
reason the discretum of CFTs described by integer values of k and N coalesces into a 
fixed line in the large A^ limit. 

We emphasize that a variety of such non-supersymmetric fixed lines exist in three 
dimensions. Choosing the fermions that transform in, say, the adjoint representation 
of U {N) gives one such fixed line of theories. Choosing fermions that transform in the 
bifundamental of U{N) x U{M) yields another example - one that preserves parity 
when N = M and the two Chern-Simons levels are equal and opposite |^ 

The study of fixed lines of large A^ Chern-Simons theories with matter and their 
bulk duals appears to be an interesting programme. For examples with a large amount 
of supersymmetry this programme was spectacularly initiated by ABJM [s] and carried 



forward in several follow up papers 10 12 . In this paper we initiate a detailed study 



of perhaps the simplest of the non-supersymmetric fixed lines - the theory of a single 
fundamental fermion coupled to a U{N) level k Chern Simons theory turns out to be 
surprisingly tractable at all values of A, for reasons we now explain. 

The Chern-Simons coupled gauge field has A^^ components. Superficially, the large 
A^ limit of these theories is governed by the summation over a complicated web of planar 
graphs. The complexity is illusory, as pure Chern-Simons theory has no propagating 
degrees of freedom - the only propagating degrees of freedom in our system are the 
fundamental fermions. Consequently, the theory we investigate is a vector model with 
A^ degrees of freedom. Large A^ limits of vector theories are much simpler than their 
matrix counterparts, and sometimes prove to be exactly solvable. Indeed, both in 
terms of diagrammatics and canonical structure, the theory we study bares a close 
resemblance to 't Hooft's solution of two-dimensional QCD in the large- A^ limit using 
light-cone gauge |T3) 

The simplicity of the vector model in the large A^ limit, combined with the crucial 
choice of light-cone gauge, allows us to derive several exact results (valid to leading 



^The ease of construction of conformal field theories in d = 3 intriguingly suggests that it is 
particularly easy to construct non supersymmetric quantum theories of gravity in d = 4. 
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It has been conjectured recently 14 that the appropriate two dimensional analogs of large N 



vector models are VFjv minimal models. Its duality with higher spin gauge theories and in particular 



Vasiliev's system in AdS^ are further explored in 15 - 19 
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order in but to all orders in A.) One of our key results is the exact expression 
for the free energy of our theory as a function of A and temperature. In obtaining 
this result we find that our line of fixed points exists only if < A < 1. We also 
enumerate the spectrum of "single trace" primary operators and demonstrate that 
their scaling dimensions are not renormalized as a function of A. We present an integral 
equation whose solution in principle determines all correlation functions of all these non 
renormalized operators - though we leave a detailed study of this integral equation and 
its solutions for future work. Based on the analysis of the spectrum of primary operators 
and explicit computation of some three point functions at low orders in perturbation 
theory, we also make several statements about the as-yet-unidentified holographic dual 
to our theory, which, as we discuss, must be a higher spin gauge theory. In the rest of 
this introduction we will summarize the key results of our paper in more detail. 



1.1 Outline of the paper and summary of results 



A summary of our main results is as follows. 

In section [2| we first calculate the exact fermion propagator for the theory on in 
light-cone gauge. Our result, which is valid to leading order in the ^ expansion but is 



exact in A is listed in (2.23). We then turn to the finite temperature partition function 
of our theory, i.e., the partition function of the theory on x S"^. We demonstrate 
that this partition function is completely determined by the self energy in the exact 
fermion propagator on x S*^. This self energy obeys a nonlinear Schwinger-Dyson 



integral equation (2.66), which, by an apparent accident of numbers, turns out to admit 



a remarkably simple exact solution (2.80). Using this solution, we find that the free 



energy of our theory as a function of temperature and A, in a box of volume V2 (which 
is taken to be very large) is given by 



NV2T^ 
67r 



A 



A 



6 / dy y\n{l + e-y) 



where c is the unique real solution to the equatior]0 



c = 2A In ( 2 cosh 



:i.2) 



(1.2) has no real solutions for |A| > 1; indeed, our fixed line of theories exists only in 



the interval |A| G [0,1). 

The fact that the free energy is proportional to —NV2T^- is an immediate conse- 
quence of the extensivity, conformality and large counting (disc diagrams dominate 



Though it is not manifest the free energy in eq.(l.l) is an even function of A, as is required by 
the invariance of the free energy under parity. 
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the free energy at large N), and could have been asserted on general grounds. The 
nontrivial part of (1.1) is the function of A, 



■ .1 - A 



+ 6 / dy y\n [I + e-y) 



that multiplies the factor —V2T^N. About A = this function has an analytic expan- 
sion in even powers of A. As |A| increases in [0, 1), this function decreases monotonically 
from the free value J;:C(3) to zero (see Fig [l] in section [2|. h{X) is a measure of the 
number of available states in our system; the fact that it monotonically decreases im- 
plies a thinning of degrees of freedom at stronger coupling that is taken to an extreme 
as A — )■ 1, when h{\) vanishes. Recall that the ABJM theory exhibits a similar phe- 
nomenon as A is scaled to oo. 

In the strict large limit we expect all transport properties of the finite tem- 
perature phase to be governed by a collision-less Boltzmann transport equation (the 



analogue of equation 9.66 of 22 in the critical 0{N) sigma model in 3 dimensions, 
see Section [s]) In other words the finite temperature system appears to behave like 
a collection of free fermions with a A and temperature dependent dispersion relation. 
This conclusion is also suggested by the fact that the exact finite temperature fermion 



two point function (2.80) has an extremely simple analytic structure as a function 
the Lorentzian frequency u. The propagator has no cuts on the real u axis The only 
singularity in this propagator in the complex u plane, are the poles located at 



In section |3| we study the spectrum of those operators whose dimension stays 
fixed in the large A^ limit. All such operators are constructed as a product of "single 
trace" operators, where by "single trace" we mean operators such as ip^ipi, which are 
formed out of the contraction of a single fundamental fermionic index with a single 
antifundamental fermionic index. We demonstrate that the spectrum of single trace 
operators is not renormalized as a function of A at leading order in A^. The set of 
single trace primaries is given, as in the free theory 



28 



by a single "current" j/*-* 



of dimension s + 1 at every spin s = 1, 2, . . . , oo together with the single dimension two 
scalar ipip. 

In the strict large A^ limit, all the operators listed above except for the scalar are 
primaries of short representations of the conformal algebra as we now explain. Short 
representations always have null states: the null states for a dimension s + 1, spin s 
current are simply the states formed out of the divergence of the current. While the 

^We thank K. Damle and S. Dutta for discussions on this point. 
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(s) 

divergence of Jli-i...^^ does vanish in the free theory, at nonzero couphng the currents 
obey an equation of the schematic form 

d.J^^^ ^\jJ+yjJ (1.3) 



Although the RHS of (1.3) is nonvanishing, it is a multitrace contribution, and so 
contributes only at subleading order in when inserted into a two point function. In 
other words the currents J^*-* are effectively conserved, hence protected, within two 
point functions, at leading order in large A^. The operators J^^"* develop anomalous 



dimensions at first subleading order in (see subsection 3.4 for details). 

A complete solution of our theory requires an algorithm to compute all correlation 
functions of J^^^ at every value of A. In section |4] we demonstrate that the three point 
functions of the spin operators, with all free indices chosen to lie in the x~ directio"^ 



are determined once we solve the integral equation (4.4) that determines "quantum 



corrected" versions of the corresponding operators. In section |5] we also derive the 



integral equation (5.2) whose solution determines the full quantum effective action of 
the theory in terms of the fermionic fields. This effective action may be used to compute 
arbitrary correlators of the spin operators with all polarizations in the x~ direction. 
We leave a study of the relevant integral equations, and of the correlation functions 
of currents with arbitrary polarizations, to future work. In this paper (see section 
|4]) we have contented ourselves with a study of three point functions of the current 
operators at tree, one-loop and two-loop level obtained using perturbative techniques. 
In particular, we demonstrate that certain parity odd tensor structures of three point 



functions recently found in 35 36 arise at one-loop order. 

In section |6| motivated by the natural appearance of bilocal variables in the large- 
A^ limit (see, e.g., (46]), we briefly discuss a light-cone Hamiltonian formulation that 
enables one to express the theory entirely in terms of bilocal variables, M{x,y), that 
satisfy a Woo algebra. Using the method of coadjoint orbits, we then derive an action 
for these bilocal variables and study its large- A^ saddle point, which turns out to be 
particularly simple. 

What is the bulk dual description of the line of CFTs studied in this paper? In order 



to address this question, we first note that it has been conjectured in 27 , extending 25 



see also 26 for earlier related work), that our field theory in the free A — >■ limit 



admits a dual description as the parity-preserving "type-B" Vasiliev theory 24 . We 
will assume the correctness of this conjecture in what follows. Clearly, then, the bulk 
description of our finite A theory is some deformation of the type B Vasiliev theory. It is 
also clear that deformation must continue to preserve higher spin symmetry in the bulk 



refers to a null direction, and we are working in the light-cone gauge = 0. 
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at the classical level; any bulk deformation that explicitly breaks Vasiliev's higher spin 
symmetry classically would require the introduction of new bulk fields (representing 
the longitudinal polarizations). The fact that the spectrum of single trace operators in 
our theory is not renormalized as a function of A implies that the dual theory does not 
have new fields corresponding to longitudinal polarizations, and so presumably does 
not break the higher spin invariance classically. It follows that the dual description is 
some higher spin theory in AdS^ at every value of A, in which higher spin symmetry is 
broken at loop level in the bulk. 

Vasiliev has written down a family of higher spin theories in AdS4^, parameterized by 



a single function of one variable 24 . The cubic interactions in this family of theories, 
for instance, are governed by a phase (as well as the overall coupling constant). The 
parity preserving type B theory (dual to the parity preserving free fermion theory) is 
a special one in this family, where the phase is equal to it/ 2 in some conventions, see 
e.g. I27I. The only other parity preserving theory, known as type A model, is such 



that the phase vanishes, and is conjecturally dual to the free/critical bosonic vector 



model 25 . Theories with generic phase different from or 7r/2 do not preserve parity. 
It is natural to suggest that the dual description of the line of theories studied in this 
paper is given by a one parameter (the parameter is A) set of such generalized Vasiliev 
theories with a A dependent phase. However, the evidence for this suggestion appears 
mixed (see section [?]). A preliminary bulk study of 3 point functions of operators 
computed using the phase generalized Vasiliev theory, suggests that the result should 
be given by the schematic form 

sm^e{JJJ)F + cos^e{JJJ)B (1.4) 

where {J J J) p is the current-current-current 3 point function in the large N free Fermi 
theory and {JJJ)b is the corresponding 3 point function in the large N free boson 
theory. As we describe in section |4] below, the explicit perturbative 3 point functions 
described in the paragraph above appear to agree impressively with this prediction at 
two loops if we set 6 = |(1 — A) + C(A'^)|^ On the other hand this identification, together 



with (1.4), implies that all one loop contributions to current-current-current three point 



functions vanish. As far as we can tell this is not the case. Conformal invariance allows 



parity violating structures for the three point functions of currents 36 , and these 
structure appear to be generated, in field theory, with nonzero coefficient at order A 
(see section |4|). We are not sure what to make of this disagreement. If we accept it 



^One may speculate that this relation between 9 and A is in fact exact, and that our theory in 
the A — >■ 1 limit could be a theory of weakly coupled bosons! There is an analog in two dimensions: 



the Wn minimal model 14 has a 't Hooft coupling A, which also ranges from to 1. In the A — >■ 
limit it becomes a theory of free fermions, while in the A — > 1 limit it becomes a theory of free bosons 
(modulo subtleties with twisted sectors). 
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in the most straightforward way, it would suggest that our theory is not dual to the 
general phase Vasiliev model and that our line of Chern-Simons theory is dual to a 
yet-to-be-constructed higher spin theory that reduces, in the A ^ limit, to the parity 
preserving Vasiliev theory. It would be especially interesting to identify the bulk dual 
to our field theory in the strong coupling limit A — )• 1. However we leave the resolution 
of these issues (and several others) to future work. 

In subsection |7.3| we also note that the maximally supersymmetric extension of the 



theory studied in this paper is precisely the ABJ model 12 with one of the gauge 
groups taken to be U{1). It follows that the higher spin bulk dual to the maximally 
supersymmetric version of the model studied in this paper is a particular limit of string 
theory. 

In summary, in this paper we have initiated the task of determining an exact so- 
lution to a simple but nontrivial fixed line of quantum field theories, namely U{N) 
Chern-Simons theories coupled to fermionic fundamental matter. We have been able 
to exactly compute the finite temperature partition function. It is possible that corre- 
lation functions in this theory are also exactly computable both at zero and at finite 
temperature. Exact results for current-current and stress tensor-stress tensor two 
point functions, as functions of temperature, would be particularly fascinating to ob- 
tain. Our work also suggests several obvious generahzations, for instance to the theory 



of fundamental bosons 37 or supersymmetric Chern-Simons theory. It also throws up 
a sharp question for the dual theory: what precise deformation of the parity-preserving 
Vasiliev theory is our system dual to? We hope to return to all these fascinating ques- 
tions in the near future. 

Note added: O. Aharony, G. Gur-Ari and R. Yacoby have recently demonstrated 
that massless scalar fundamental matter coupled to Chern-Simons theory also leads to 



a large N fixed line of CFTs 37 . Using arguments very similar to those employed in 
part of section 3 of our paper, they have also argued that the higher spin currents in 
their theory are not renormalized as a function of the coupling constant. We thank 
O. Aharony for sharing his results with us prior to publication and also for several 
discussions on these issues over the last year. 



2 Free Energy on at Finite Temperature 

In this section we will evaluate the free energy of our theory (a single species of massless 
fundamental fermions coupled to a Chern-Simons gauge field) in the t' Hooft large N 
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limit. The Euclidean action for our theory is 

S='^Jtt (^AdA + '^-A'^ + J ^YD.^fj, (2.1) 

Our theory is taken to be at temperature T, and lives on a spatial whose regulated 
volume we denote as V2. 

In order to evaluate the free energy of our theory, we first fix a gauge. We work in 
the lightcone gauge A- = 0. This gauge is defined in terms of an analytic continuation 
from Lorentzian space. For this reason A_ = does not imply A+ = 0. Alternatively 
one can work in Lorentzian space and analytically continue to Euclidean space (after 



integrating out the gauge fields) in (2.24). The gauge boson self interaction term 



vanishes in this gauge, a feature that enormously simplifies analysis. Computations 
in pure Chern-Simons theory on in this Euclidean continuation of the light-cone 



gauge were done previously in 42 , 43 



In our analysis below we will encounter divergent integrals that need to be regulated. 
We choose to regulate all integrals in the scheme of dimensional reduction. More 
specifically we evaluate all integrals as follows. We evaluate 7 traces, e contraction 
etc in d = 3. This process leaves us with a set of scalar integrals. We then evaluate 
the resultant integrals by analytic continuation from d = 3 — e dimensions. This 
regularization scheme is widely employed in the previous literature on Chern-Simons 
matter theories (see e.g. [20||2l||6|[7]). It is manifestly Lorentz invariant, and also 
respects gauge invariance at least up to two- loops (see (21]). We will assume without 
proof in what follows that our regularization scheme is indeed gauge invariant. If this is 
indeed the case then the theory defined by this regularization scheme must be Lorentz 
invariant even though we work in a gauge that breaks Lorentz invariance. We will find 
some evidence for the Lorentz invariance of our final results giving some a posteriori 
evidence for our assumption that our regularization scheme respects gauge invariance. 

As we will explain below, the finite temperature free energy of our theory is com- 



pletely determined by the fermion self energy on x (see (2.52)). In order to 
evaluate the free energy we proceed as follows. As a preliminary step to our analysis, 
we first determine the exact fermion propagator of our theory on M^. We then deter- 
mine the exact fermion propagator on x S*^. Finally, we proceed to use this result 
to determine the free energy of our theory. 



^This feature is also true of the more straightforward "temporal" or axial gauge A3 — 0. Feynman 
diagrams in this gauge are, however, plagued by logarithmic divergences that we have found difficult 
to interpret and deal with. In contrast the divergences in the lightcone gauge employed in this 
subsection are relatively tame, and are easy to interpret and deal with. We thank S. Bhattacharyya 
and J. Bhattacharya for extensive discussions on perturbation theory and its divergences in Axial 
gauge. 
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2.1 Exact fermion propagator on 



In the planar limit, the exact fermion propagator in light-cone gauge receives contri- 
butions from a very simple class of diagrams. Indeed, the gluon vertex vanishes in our 
gauge, and moreover at large we can neglect fermion loops since the fermions are 
in the fundamental representation. Therefore, only rainbow diagrams contribute to 
the exact propagator. These can be resummed exactly by solving a Schwinger-Dyson 
equation which may be depicted diagrammatically as 




The planar fermion self energy. 



This picture may be transcribed into equations as follows. Let the exact fermion 
propagator be given by 



(V'(p)mV'(-?)") 



Then 



where 



S(p) 



N 
~2 



(27r)3 



5" 



X {2'Kf5{p-q). 



(2.2) 



(2.3) 



-q)) = {2^f5{p-q)G,M^'^' 

is the bare gluon propagatorj^ Here M^are is the mass term that appears in the bare 
Lagrangian. In what follows we will adjust M^are to ensure that the physical fermion 
mass vanishes (this choice corresponds to tuning the theory to the conformality) . Note 
that with our definition of S in (2.2), the sum of IPI diagrams is given by — S. 



The equation (2.3) applies in fact in any ghost free gauge in which all gluon inter- 
actions vanish. In this section we will solve (2.3) in the lightcone gauge. 



2.1.1 Exact solution of the gap equation 



In the this subsection we will solve the gap equation (2.3) in light-cone gauge (see 
Appendix |A] for a detailed listing of our conventions, propagators etc in lightcone 



^The factor of ^ in (2.3) has its origin in the i on the RHS of (A.3). 
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gauge): 



(2^ 



r7'" 



(2.4) 



Let us first better understand the matrix structure of tlie self energy S. Let A 
represent an arbitrary 2x2 matrix 

A = All + + A_7~ + ^37^ 

For use below and in later sections we define the matrix valued functions of A, Hj^{A) 
and H^{A) as 

H+{A) = 73^7+ - 7+^7^ = 2 {An+ - v4_ /) (2.5) 

H^{A) = 7^A7- - 7-^7^ = 2 (-A/7~ + A+I) (2.6) 



The gap equation (2.4) may be rewritten as 



S(p) = -z27rA 



f27r)3 



i-fi'qf, + Mbare + S(g) 



(2.7) 



In the discussion which follows we will sometime abbreviate the notation S(p) to S for 
brevity. Now let 

S = lE^Y + S// - Mbarel (2.8) 



Using 



1 



-i7'^(p^ + E^) + E, 



z7M(g^ + S^) + S, (p + E)2 + E2 



together with (2.5), we may rewrite (2.7) as 

d^q 



E(p) = -iAnX 



7+E, + zJ(g + E(g))_ 



1 



(27r)3 (g^ + E^(g))(gM + SM(g)) + E,(g)2 (p - g)^ 



(2.9) 



Plugging (|2.8|) into the LHS of (|2.9|) and equating the coefficients of linearly indepen- 

(2.10) 



dent matrices, it follows immediately that 

E_ = Eg = 0. 
And that 

d^q E/ 



E+(p) = -AttX J 

Ulip) - Mbare = 4vrA J 



1 



(27r)3((g + E(g))2 + E,(g)2) {p - q)+ 
d^q g_ 1 



(2.11) 



(27r)3((g + E(g))2 + E,(g)2) (p - qY 
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What can we say about the dependence of and on pi First note that 

the that the RHS of the two equations in (2.11 ) is independent of . It follows that E is 



a function only of the in plane momenta and p^ but is independent of p^. Rotational 
invariance in the 12 plane and requirement of conformality (i.e. the requirement that no 
mass scale enter the physical propagator) then together completely fix the momentum 
dependence of E_|_ and E^: 



S/(p) = foPs 



P 9o 



where 



Ps 



pI+pI 



V2\p-'\ = V2\ 



p 



(2.12) 



(2.13) 



and /o and 5^0 are dimensionless numbers (that are functions of the coupling constant 
A). 



Plugging ( |2.13D and ( |2.10D into (|2.11|) we find 

4n\ 



9o 



fo\p\- Mbare 



P 

AnX 



(2^ 



Qsfo 



ql + qK^ + 9o + f^) 



{p - qY 



(2.14) 



(27r)3 ql + q2(l + g^ + f^) - q) + 
We will now proceed to determine the numbers Qq and /q as functions of A. 



We first note that the integrals on the RHS of (2.14) are (power counting) linearly 
divergent. In order to proceed we need to regulate these divergences. We adopt a 
regulator that is manifestly Lorentz invariant as well as plausibly gauge invariant. Our 
regularization procedure is simply the following: we analytically continue all diagrams 
to 3 — e dimensions. Two of these dimensions span the 1-2 plane. We integrate over 
the remaining 1 — e dimensions using the formula 



X 



AH 



(2.15) 



Here A is a number of order unity which is easily computed. However, as we will see 
below, none of the integrals we compute in this paper have a ^ type divergence (this 
corresponds to an absence of logarithmic divergences, were we to use a momentum 
cutoff). The only effect of the dimensional regularization cut off procedure, employed 
in our paper, is to discard linear (and below also cubic) divergences in a gauge and 
Lorentz invariant manner. For that reason, when we take e — t- at the end of the 
calculation, we effectively set A to unity. Hence we immediately set A to unity instead 
of carrying it around in our computation. The regularization procedure we employ here 
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is essentially the dimensional reduction scheme used in 21 , adapted to our light-cone 
gauge. 



Using ( |2J[5| ) in ( pl4| ) yields 

90 = - 



27tX 



fo 



V^ + 9o + fSJ (27r)2g|(p-g) 



foPs - Mbare = SttA- 



(2.16) 



In order to do the integrals we move to polar coordinates in the 12 plane. The integrals 
we need to evaluate are 



d\ 



q- 



d'q ql- 2 



qrdqTde P^'^'^-^^ 







g2 + p2 _ 2p,q, cos 9 
q^dql'^e- P-^-^-^'^'^ 







Is +P'i ~ "^Psqs cos 6 



{2nyp- qlip - q)+ {27r)^PsJo 
Using contour techniques it is not difficult to verify that for q > p 

de 2tt 



g2 _|_ p2 _ 2pq cos 9 q^ ~ p"^ 
2^ rf^cos^ p 27r 







It follows that 



27r 



dd 



g2 _|_ p2 _ 2pq cos 9 qq^ — p"^ 
Psqs cos 9 - q"^ 



2n 



d9 



2n 



d9 



2tt 



d9 



Qs + pI 


— 2psqs cos 9 


PsQs 


cos 9 — ql 


Qs + pI 


— 2psqs cos 9 


Psqs 


— ql cos 9 


Is + 


— 2paqa cos 9 


Psqs 


— ql cos 9 


Qs + pi 


— 2psqs cos 9 



{qs < Ps) 
-2tx {qs > Ps) 
271 (qs < Ps) 
{qs > Ps). 



It follows that (2.16) reduces to 

^0 = - 

foPs - Mbare = ' 



A/o 



'1 + 90 + 



A 



1 + ^0 + fo J p. 



X 



'1 + 90 + 



Ps 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



It follows from (2.20) that Mhare = (this is a consequence of our use of dimensional 
regularization; Mbare is linearly divergent in a cut off regulator, as is clear from the 
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second of (2.20)). The remaining equations reduce to 

A/o 



^0 

/o 



'1 + 90 + 
A 



(2.21) 



V^ + 9o + fi 



The solution to (2.21) is remarkably simple 



fo 
90 

90 + fo 



A 

-A^ 




(2.22) 



In other words, the self energy receives contributions from one and two loop graphs 
but not at any higher order in perturbation theory! This completes our solution of the 
gap equation. 



We emphasize that our final result (2.22) depends in a crucial way on our choice 



of regularization scheme. Were we, for instance to regulate the integrals (2.14) by 





_ (p-i) 

modifying the gauge boson propagator with a Gaussian damping factor e ^ then 
we would have found go = ^ (from the integral over the angle in the vector p — q). It 
is of course quite clear that a crude cut off on the gauge boson momentum does not 
preserve either gauge or Lorentz invariance. We hope on the other hand that the more 
sophisticated dimensional regularization scheme preserves both these symmetries. This 
assumption, which is as yet unproved, is the main weakness in the analysis presented 
in this section. We will return to this point at the end of the section. 

In summary, the exact fermion propagator, at leading order in large A^, is given by 



1 



ipsl^ + ip-l~ + ^(1 - A2)p+7+ + Xp^ 



x{27i)^{p-q). (2.23) 



(Here m, n are colour indices and we have suppressed the spinor indices). 



2.1.2 Rederivation of the gap equation as a Schwinger-Dyson equation 



The starting point of our analysis in this subsection is the path integral representation 
of the partition function 

where the action 5* is the Euclidean space lightcone gauge action for our field theory, 
listed explicitly in (A. 6). The gauge fields appear quadratically in in (A. 6) and may 
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be integrated out. Integrating out the gauge field from (A. 6) yields the path Integra' 



10 



where the action S is now given by 

S = i 



+ 



d?p — 

2-ni r d^p cP'r d?q 1 
~k J (27r)3 (27r)3 {2Tifq^ 



(2.24) 



Let us pause to note that, starting from the bifocal action (2.24), one can conve- 
niently derive the Schwinger-Dyson equation for the fermion self-energy (as in [46^ ) 
via 

5 







D^D^(6l6'{p'+p) 



-rip') e 



(2.25) 
(2.26) 



which gives the following relation involving the exact fermion propagator and four-point 
functions, 

{27Tf6'{p'+p) =tp,Y {ripWip')) 



+ 



2m 
~Y 
2m 



d?r d^q 1 
(27r)3 (27r)3 ^ 



7"^ (V^a(p - q)r{-rWnj + qWij>')) 



(27r)3 (27r)3g+ 



(2.27) 



In the large- limit, this factorizes to yield 

{np)r{p')) ^ 



i2TTy5\p'+p) 



1 27ci f d^r d^q 1 



iPtMl^ k J (27r)3 (27r)3 q^ 



7+ (Mp - q)r{-r)) (7'^n^(r + qWip')) 



+ 



1 2m f d'r d^q (^^(^ _ ^+ (^^(, + ^)^'^(p')) 



iPix'y^ k J (27r)3 (27r)3 q+ 



which, upon substituting (2.2), gives the gap equation: 

d^q ] 



E(p) = -2m\ j 



{2-nfp+-q 



+ ' 



i{q^)Y + ^{q)J ■ 



(2.28) 



(2.29) 



^"in integrating out the gauge field, we absorb the factor coming from the determinant of the gauge 
field kinetic operator in the normalization of the path integral. In other words, we normalize the pure 
Chern-Simons partition fmiction to 1. 
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2.1.3 Rewriting the field theory as a path integral over singlet fields 



In this subsection we will reformulate the path integral that evaluates the partition 
function of our field theory as a path integral over singlest fields. The new path 
integral is weakly coupled in the large limit (the action in terms of the new variables 
is proportional to A^). The gap equation (2.4) follows as the classical equation of motion 
of this large N action. 

While the work out of this subsection is considerably more complicated than that 
of subsection (2.1.1) it has one significant advantage; it reveals how the solution of the 
gap equation is related to the value the partition function of the theory. While the 
value of the partition function is of no physical significance for the theory on M.^, it is 
of great significance on ]R2 X (as it determines the thermal partition function of the 
theory on M^). For this reason the results of this subsection will prove very useful in 
our discussion of the finite temperature partition function in the next section. 

We now introduce some convenient shorthand notation. Let 

1 

N 



M{P,q) 



q) 



(2.30) 



M is a 2 X 2 matrix in spinor space but a singlet in colour space. While one of its 
arguments, P, is a 3 momentum, its second argument g is a 2 momentum (in integral 
on the RHS of (2.30) is over the 3 component of q). (2.24) may be rewritten as 

^^ 'J (2vr)3 

27riN'^ r d^P (Pq cPq' 1 



(]?p - 



k 



(27r)3 (27r)2 (27r)2 (g - q')- 



-Tr {M{P, q)-f+M{-P, g')7^) 



(27r)3 

TTiN^ r d^P d^q d^q' 



k 



(27r)3 (27r)2 (27r)2 (g - g')" 



-Tr [{M{P,qh+M{-P,q'h^) - {M{P,qW M{-P,q')^+)] 

(2.31) 



(the flip in sign is due to the fact that we had to take one fermionic field through three 
others). Expanding the matrix M in a complete basis of 2 x 2 matrices 

M+7+ + M_7" + M37^ 



M 



Mil 



(2.32) 



we find that (|2.31|) reduces to 
S = 



+ 



d?p — 

j^i^i-ph^Pf^i^ip) 

iriN'^ r d^P d^q d^q' 



(2.33) 



(27r)3 (27r)2 (27r)2 (g - g')" 



M_(P,g)M,(-P,g') 
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where we have used 

Tr([7-,7l7') =-4. 
Note in particular that M_(_ and M3 drop out of this expression. 



(2.34) 



We will now rewrite the interaction term (the term quadratic in M) in (2.33) in 
terms of a Lagrange multiplier field 

S = S+7+ + S/J. (2.35) 

where S will turn out to be the self energy of the fermion field. To this end we define 
the "inverse" Greens function G~^{p) by the requirement that 



;G-\p-q)- 



{2ny6\p - r) 



(27r)2 ^^(g_r)+ 
Note that is an odd function of its argument. Note also that 

^ G'-^(r-p) = (27r)V(g-p) 



(2.36) 



(27r)2 (g-r)+^ " ^'^'^^^ 

These are the only properties of G^^ that we will need in this paper; in particular we 
will never need the explicit form of the function G~^. 

Now it is obvious that 



(2.38) 



where we have chosen 



E = 2 X 



N 



d^P (Pq (Pq' 



(27r)3 (27r)2 (27r)2 
(s4P,g)-4vr^A j ^.MAP,r)^^^^G-\q-q') 

S,(-P,g')-47rzA j 



(2.39) 



(27r)2 [q' -r')+' 

Note that ii^ is a function of the two new Lagrange multiplier fields ZI_ and S/. The 
path integral in the denominator in (2.38) is simply a number of order unity and we 
will omit to write it in the equations that follow. The effective action in the numerator, 
S + E, evaluates to 



S = i 
+ 
+ 



(Pp — 

ip{-ph''p^,ip{p) 



(27r 

d^P d^q -.P „ , 

(2^(2^^^2-^)^(-^'^)^(2+^) 



(2.40) 



N f d^P d^q d'^q' 



27TiX J (27r)3 (27r)2 (27r) 



S+(P,g)G-i(g-g')S/(-P,g') 
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In the second line above we have used the fact that 

1 f dq-i 



-2 S+A-f_ + S/M/ 



-Tr^M 



N I 2.r^^^^^ = ^/^^^^ 



Using (2.34) the last line in (2.33) may be rewritten as a trace, yielding 



S = i 
+ 



d^P d^q -,P ^ , ,,P 

(2;0^(2^^(2-^)^(-^'^)^^2+^) 



(2.41) 



N r d^P d^q d^q' 
8mX J (2^(2^(2^ 



G-\q - q')Tr (7^S(P, g)7'S(-P, q')) 



The action ( 2.41^ may be rewritten as 



d^P d^n P P 

(^^( 2 - ^) {i'27r?6\P)^Yq, + S(-P, g)) ^(- + g) 



N 



^ ^^,G-^(g-g')rr(7-S(P,g)7^S(-P,g')) 



(2.42) 



TTzA J (27r)3 (27r)2 (27r) 



The dependence of (2.42) on fermionic fields is quadratic, so the later may be integrated 



out. Performing this operation yields 



Z 



(2.43) 



where 



S = -A^Trln {{27if6%P)i-f^q^ + S(-P,g)) 
N r d^P d^q d^q' 



SiriX J (27r)3 (27r)2 (27r) 



G-\q - q')Tr {^-^{P, g)7'S(-P, q')) 



(2.44) 



Notice that (2.44) is written purely in terms of singlet fields, and is multiplied by an 



overall factor of A^. (2.44) represents an exact rewriting of the partition function of 



the original theory as a partition function over the singlet fields S; this path integral 
is weakly coupled in the large N limit. 



The action, (2.44), is somewhat formal, as it is written in terms of a determinant 



over an infinite dimensional matrix. However the equivalent of (2.44) is much simpler 
for translationally invariant S configurations of the form 



S(P,g) = (27r)3<5(P)S(g). 
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The form of this action is perhaps most clearly obtained by retreating to (2.42), which 
reduces, for translationally E configurations to 



S 



(277)3"' L / ■ yr-nryr-; g^.^ j (27r)2 (27r)2 

V here is a factor of the volume of spacetime, and we have used 
in the last term of the last line. 



(2.45) 



Integrating out the fermions in (2.45) yields a very explicit special case of (2.44) 



S=-NV 



(27r) 



-Trln(^7'^g, + S(p))-^ / , J, ^J,, G-\q-q')TT (7"S(g)7^S(gO) 

(2.46) 



-KiX J (27r)2 (27r)^ 



While all terms in the action in (2.46 ) are proportional to A^, the fields in that action are 



gauge singlets. At leading order in the large expansion it follows that the free energy 



for our theory may evaluated simply by minimizing (2.46) w.r.t. S. The variational 
equation we encounter in this minimization process is 



Tr 



L(2vr)^ 



.<5E[g] 



dq^ 



1 



2tt ij^^q^ + S 
(2.47) 



In terms of the function H_ defined (2.6) 
d^q 



Tr 



L(2vr) 



d^q' 



(27r)2 8mX 



The equation (2.48) is of the form 
where 



d\ 
{2nf 



5E[g]i?(g) 





(2.48) 
(2.49) 



dq^ 



1 



2tt i'j^^q^ + S 



As 5S is an arbitrary matrix of the form (2.35) it follows that 



i.e. that 



B4q) = Bjiq)=0 



HUB) 
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(see (2.5)) Using the fact that 



i/+(iJ_(S))=4S 



and integrating both sides of (2.48) against the kernel 
property of the function G~^, it follows from H+{B)) = that 



and using the defining 



S(p) 



-2mX 



r 



7 -7 



;7^ 



i'^^'q^ + T. J {p - q) 



(2.50) 



in precise agreement with (2.4). 



The value of the Euclidean action on the saddle point, (2.46), may be rewritten as 
d^q 



S = -NV 
= -NV 
= -NV 



(27r)3 
d^q 

(2^ 
d^q 

(27r)3 



Trln {^Yq^. + ^iq)) + 



NV 



d'^q d^q' 



IQixiX J (27r)2 (27r)2 



G-^(g-g')Tr(iJ-[S(g)]S(g')) 



Tr 
Tr 



ln(^7^g^ + S(g)) + -/f4S(g))/f, 



ln[^7^g^ + S(g)]--E(g) 



i-^^'q^ + S(g) 



27/^ + S(g)y_ 



(2.51) 



where we have used the equation of motion in going from first to the second line. In 
going from the second to the third line we have used the fact that for an arbitrary 
matrix A 

Tr {H_{^)H+{A)) = -4Tr (SA) . 



2.1.4 Diagrammatic expansion for the vacuum energy 



The result (2.51), expressing the vacuum energy in terms of the exact fermion self- 



energy S(p), has a clear diagrammatic interpretation. The perturbative planar dia- 
grammatic expansion takes the form 



vacuum energy 



+ 





+ 



One may verify that, in terms of exact planar fermion self energy, the above expansion 
can be written as 
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vacuum energy 





+ 



+ 



+ 





where the shaded blob represents the exact planar fermion propagator [i^ + S) ^. In 
particular, the second term in the last line of (2.51), — |Tr[S(i^ + is precisely 



such that the correct symmetry factors are restored. 



2.2 The finite temperature theory 



In this section we study the logarithm of the path integral of our system on x 
where the circumference of the is taken to be /3. This path integral determines the 
free energy of the field theory at temperature T = 

The formulas that determine the path integral on x are straightforward gen- 
eralizations of the formulas on M^. Every equation in subsection 2.1.3 carries through 
with the replacement 



/ ^/(P3) - ^ E /( 



27rn, 



V ^V2l3 



so that 



V 



(2^ 



(27r)2 



E 



In particular the Euclidean action is given by 



ln[^7^g^ + ST(g)]--ET(g) 



where T, the temperature is and the function ^^(g) obeys the gap equation 



(2.52) 



St(p) = 
where n is an integer and 



1 



Z7^g^ + ST(g) 



27r(n+ 



1 



7^ 



1 



(2.53) 
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In order to determine the free energy at finite temperature T, we need to solve tfie 
gap equation (2.53) and plug the solution into (2.52). We take up these exercises in 
turn. 



2.2.1 The finite temperature gap equation 



As in subsection 2.1.1 it follows immediately that is a linear combination of 7''" 
and J, and that it is independent of ^3. Rotational symmetry and the constraints of 
conformality then imply 

for some as yet unknown functions f{Pps) and g{Pps)- Note that the new dimensionful 
scale /3, now allows / and g to be functions of ps, generalizing the pure numbers /o 
and Qq of the previous section. The zero temperature results of the previous subsection 
imply that 



The analogue of (2.11) is 

Psf{PsP) - Mbare = 47r^ J 



hm f{y) = /o = A 
lim g{y) = go = ->? 

J/— >oo 



(2.54) 



(27r)2 /Mn+|) 



9{PsP)p~ 



-Att- 



X 



+ gKl + ^?M) + l/M)P) 
f{qsP) 1 



^ (p - q) 







(p - qY 

(2.55) 
- (2.56) 



The summations in these equations are easily carried out using the formula 

1 



{n + \y + a^ + q^ 



\l-e 



tanh(7r|a|) 



(2.57) 



have set e to zero in every place where it will be inessential for regularization) yielding 



(this is the analogue of (2.15) in the previous section - as in the previous section we 

ielding 

(2.58) 
(2.59) 



d'^qq's 



tanh 



Pqs 



q- 



9P 



-2nX 



(Pqq 



(2vr) 



tanh 



l^qs 



'l + 9 + P 



qs{p-q)+^l + 9 + P 
f 1 



1 + 9 + pip- qy 
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where we have left imphcit the fact that the / and g are functions of Ps on the LHS of 



(2.58)(2.59), but are functions of Qs on the RHS of the same equations. 



perform the angular integrals to obtain 



In each of (2.58) and (2.59) we move to polar coordinates and use use (2.18) to 

'f3q 

'Ms 



9 



-2A 



p 

P /yl— e 



fPs = I Q 'dqstanh ( -^^1 + 9 + P . r- - 



(2.60) 



(Is 'dQs . , 

^ — tanh , 

,2 \ 2 



P. 



V^ + 9 + P 



f 



'l + 9 + P 



(2.61) 



Adding and subtracting q from the integrand of (2.60) and doing the integral on the 
trivial piece we find 



/ = A-A 



g = -2\ 



dqs 

Ps 



pI 



tanh(^Vl + (7 + /2) ^ 
V^ + 9 + P ^ 

tanh(^^/rTiT7^)/ 
V^+9 + P 



In terms of the variable 



X 



P 



'tanh (y^^l+g(yx) + f{yxy 

f[y) = A - A / dx\ ^ - 1 

'i \ a/I +5'(ya;) + /(yx)2 



'1 / tanh ( f Vl + g{yx) + f{yxY ) f{yx) 
g{y) = — 2A / xdx 



where the variable 



Equivalently 



V^ + 9{yx) + fiyxy 
y = pP- 



f{y) = A 



9{y) 



^ .ex. /tanh(fVl + ^?(^) + /(x)2) ' 
yjy "^""y ^l + g{x) + f{xf ^ 

A fy tanh(fv/l + ^7(x) + /(x)2) /(x) 



-2— / xdx 
y^ Jo 



^l + gix) + fixy 



(2.62) 



(2.63) 



(2.64) 



(2.65) 



(2.66) 
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2.2.2 The exact solution 



Quite remarkably it is possible to find tlie exact solution to (2.66). We start with (2.66) 
written in the form 



y U{y) - h) = -A / dx 



y'^giv) = -2A / xdx 
Jo 



tanh I I y/l + g{x) + /(x)- 



^l + g{x) + f{xy 
tanh (fv/l + ^?(x) + /(x)2) f{, 



(2.67) 



y/l + g{x) + f{xy 
Differentiating both equations w.r.t. y we obtain 



yf{y) + f{y) = A 



tanh(|v/l + ^7(2/) + /(y)^ 



y(?'(l/) + 2^(l/) = -2A/(i/) 



tanh(|Vl + ^7(l/) + /(y)' 



(2.68) 



Multiplying the first equation by 2/ and add it to the second equation, we cancel 
the RHS and obtain 

(2.69) 



y^{9 + f) + 2{g + n=0. 



From this we solve 



9iy) + f{yf 



y 



(2.70) 



where c is a constant. Now the first equation in (2.68) becomes simply 



yf{y) + f{y) = A- 



tanh(|. /l 



1 + 



(2.71) 



Integrating we have 



yJo \/T+lf y y 



cosh 




cosh 


2 





+ 



(2.72) 



9{y) 



f{yf 



Here c is another integration constant. 

While the solution to a differential equation depends on integration constants, the 
solution to an integral equation is unique (it does not have undetermined integration 
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constants). The appearance of c and c in our solutions above is an artifact of our 
having solved by converting the integral equation into a differential equation. The 



integral equations (2.66) are actually solved by (2.72) only for a particular choice of c 
and c. 

We first note that the function / must tend to /o at large y. This is automatic 
in all our solutions (it does not impose any constraints on c or c). However a further 
requirement is that the expansion of / about this constant value (at large y) should 



start at ^ rather than ^ (this follows immediately upon plugging (2.70) into the RHS 



of the first of (2.67); the RHS of that equation is manifestly ^ )• The requirement 
that 

f{y) = /o + 0{l/y') (2.73) 



determines 



c = 2A In I 2 cosh 



Let us now turn to the small y behaviour of / and g. At small y 



(2.74) 



f{y) 
9{y) 



+ 0{y), 



+ 0{y'). 



(2.75) 



Plugging ( |2.70[ ) into the RHS of the second of (2.67), however, we find that the RHS 
evaluates to 0{y^) at small y implying that g{y) = 0{y^) at small y. It follows that 



Plugging this relation into (2.74) yields the following equation for c 



c 

2A 



In I 2 cosh 



(2.76) 



(2.77) 



Note that c is an odd function of A. |c| is a monotonically increasing function of A, 
which diverges at |A| = 1. (2.77) has no solution for |A| > 1, indicating that the theory 



does not exist for \k\ > N. At leading order in at small A we have 



As A approaches unity we have 

~ 1 2 

c = m 



1-A) 



£ = 2Aln2 + C(A=^). 



In ( In^^ 1 + C ( Inlnln ^ 



1-A 



1-A 



(2.78) 



(2.79) 
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In order to physically interpret the divergence of c as A — )■ 1 note that the exact thermal 
propagator has a pole whenever 

In other words c has a simple physical interpretation; it is the thermal mass of the field 
i{j in units of the temperature. It follows that the fermion thermal mass diverges in 
this limit A — ?■ 1. 



Using (|2.77|) and (|2.76|) we may rewrite our solutions for / and g as 

a/c2 + ?/2 



2A 



In 2 cosh 



C2 



(2.80) 



9iy) = -,- fiyV- 



with c given by (2.77). In the large y limit 



fiy) = \\ll + ^,+0{e-y), 

9{y) = -X' + ^^ii^ + o{e-y) 
y 



(2.81) 



while at small y 



f{y) 
9{y) 



— I — — tanh 

y 2d 



+ 0{y' 



-Atanh(-)+0(|/2). 



(2.82) 



2.2.3 Free energy as a function of temperature 



As we have explained above, the path integral of our theory on the manifold x 5*^, 
with the circumference of the equal to (3, is given by e^"^^ where 

72. 



\n[tYq^. + ^T{q)]-^^T{q) 



Let us define 

= -NV2^ 



cPq 



(27r)3 

Then the partition function 



Tr 



ln[27"g^ + S(g)]--S(g) 



i-ii'q^ + S(g) 



Tre 



-m 



(2.83) 



(2.84) 
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of our system in a flat spatial box of volume V2 is given by 

In Z = So — St, 

so that the finite temperature free energy, F{T), of the theory is given by 

FiT) = 

We will now proceed to use our exact solution to the finite temperature gap equation 
to compute St — Sq. 



2.2.4 Explicit evaluation of the free energy 



In order to find an explicit expression for the free energy, we find it convenient to use 



the expressions in the second line of (2.83) and (2.84). St — 5*0 may be written as 



S^-So = -NV2j2 I ^,Tr\n{trq, + ^T{q)) + NV2 J ^Tr In (^7% + S^g)) 



W./3/^Trln.^^^'^^^ + ^-(^)) 



(2vr) 



(2vr) 



27^g^ + ET(g) 



(2vr) 



(2.85) 



The integral on the first line of (2.85) is convergent and evaluates to 



E / (0 1- (^3 + q! + c'T^) + NV,P / ^ In (^3^ + + c^T^) 



-2NVo 



-2NT^V2 

2NT^V2 
2^ 



(2vr)2 



In 1 + 



In 1 + e 



(2.86) 



/■oo 

/ dy y\xi{l + e~y) . 

J\c\ 



The second line of (2.85) has a linear divergence, which however disappears in our 



dimensional reduction scheme 



(2t) 



2|~|3 



^n^ 



(2.87) 



27r2 



dyy'-'\n(l 



2|s|3 



Qn 
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where in the last step we integrated by parts and used 



TT 



in the hmit of small e. Let us now turn to the third line of (2.85). The second term in 



the third line simply vanishes under dimensional regularization. The first term in the 
third line may be evaluated as follows 



NVo 



d\ 2qlP 



(2vr) 



ql + ql + T^c 



73. 



NV2 

4 



+ 



d\ {2qlf{Pqs) + qlg{Pqs))t^n\i , 
(2vr)2 



7.2 



d'q {2qlf\q,) + qlg{q,))t^nh 



- A2 (£2 + g2 



C2 



(2i 



(27r)2 

d^q \^ (c^ + g2) + g2 



C2 



(2vr)^ 



£2 



The first term in (2.88) is finite and evaluates to 

NV2T'^ 



Air 



151' 



A2 
6 



1 



1 



6|A| 2|A| 



(2.89) 



The second term in (2.88) is divergent. In dimensional regularization it evaluates to 

NV2T^\cf 



A2 
6 



(2.90) 



The third line of (2.85) is given by the sum of (2.89) and (2.90) and evaluates to 

NV2T^\c\^ 



67r|A| 

Putting it all together we find that the free energy is given by 



NV2T^ 
Gtt 



l~|3 
C 



1 - |A| 



6 / dy y\n{l + e-y) 



(2.91) 



(2.92) 



Although it is not manifest, the free energy can also be written in the form (1.1), and 



is an analytic function of A in the interval (— 1, l)j^ To see this more explicitly, note 

Theories with massless bosons usually do not have an analytic expansion of their free energy in 
terms of the coupling constant. Such non-analytic behavior in Chern-Simons-matter theories with 
massless bosons was observed explicitly e.g. in [t], |38| . This non-analyticity has its origin in infrared 
divergences which are absent in our theory because the only propagating degrees of freedom are the 
fermions which lack a zero-mode along the thermal circle. 
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F 




that we can write the term in square brackets in (2.92) as 

|3 



\C\ 

'\x\ 



- Icf + eJ^ dy y\n (l + e"^) +6 J^^ dy(^^ -y\n (^2cosh 



c3 



2^(3)+ A 



6 / dy y\n(2 cosh 



y 



(2.93) 



where in the last hne we used |c|^/|A| = c^/A, as follows from (2.77). Note that the 
integral in the last term is clearly an even function of |c|, and so the absolute value 
can be omitted. This shows that one may effectively rewrite (2.92) as in (1.1), which 
is manifestly analytic. Indeed its small A expansion is given by 



3C(3) 2(log2) 



An 



371 



(log 2)^ 

271 



A^ + C(A« 



(2.94) 



Note that it contains only even powers of A, consistently with parity. A plot of the free 
energy as a function of A is given in Fig. [T] We see that —F decreases monotonically 
from the free field value —F = ^^^^'^^ ({3) to zero at A = 1. 



In the limit that |A| — > 1, |c| is large and the integral in (2.92) may be approximated 



by 



POO 

/ dyy\n{l + e'^) = \c\e-\'\ + e"!'! + 0{e-^''\] 

J\c\ 



At leading nontrivial order this evaluates to 



1 - |A| 
2 
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and the free energy is given by 



F = ^ ^^-TJ - \X\) ^|g|3^3g2^^(g)] (2.95) 



where c is given by (2.79). 



2.3 Consistency of our gauge and regularization scheme 

As we have emphasized on muhiple occasions, we have obtained the beautiful result 



(2.92) by employing a rather unusual gauge (a lightcone like gauge in Euclidean space) 
and the regularization scheme of dimensional reduction, which we have assumed pre- 
serves the gauge invariance of our theory. In this section we list the independent 
evidence that the procedure employed in this section defines a sensible and Lorentz 
invariant theory. 



1. The exact fermion propagator (2.23) develops poles only when p = 0, a condition 



that is Lorentz invariant. This is a necessary condition for the Lorentz invariance 



of fermion scattering processes in our theory. The Lorentz invariance of the 
poles of the propagator is far from automatic, and is in fact violated in several 
regularization schemes. In, for instance, the regularization scheme described just 
below (2.22) we find go = but fo^O (in fact /o obeys the equation /o = ^ 



with this regulator). The poles of the fermion propagator with this regulator occur 
at J93 + p^(l + /o) and are not Lorentz invariant. 

2. In Appendix [C] we have demonstrated that the expectation value, at one loop, of 
the gauge invariant Wilson line operator is Lorentz invariant and moreover agrees 
with the result obtained in the manifestly Lorentz invariant Feynman gauge. 

3. Our results above indicate that our theory is infinitely strongly coupled at A = 1, 
and that the theory does not exist (at least as a conformally invariant theory) 
for A > 1. There is in fact a very simple interpretation of this result. It is 
well-known that the bare Chern-Simons level can acquire a finite shift at 1-loop 
order in perturbation theory. This effect is regularization dependent. It was 
demonstrated by [21] that the Chern-Simons matter theory defined using the 



^^Particle scattering is, strictly speaking, not well defined in our theory (or any conformal theory) 
due to infrared divergences. We could cure these divergences in our theory softly breaking conformal 
invariance with a fermion mass. Scattering processes in this deformed theory are presumably well 
defined. The poles of the mass deformed theory are clearly physical and must be Lorentz invariant. 
As our computation of the propagator of massless theory, presented in this paper, exhibits no IR 
divergences, it must equal the zero mass limit of the mass deformed propagator, and hence must be 
Lorentz invariant. 
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dimensional reduction scheme, employed in this paper, does not acquire a 1-loop 
shift of k. On the other hand, if the theory is regulated by the addition of a 
small Yang Mills term, the bare level, which we denote kyM, gets shifted by 
sign{kYM) N . Therefore the results of [2^ imply that the two regularizations 
yield the same physical theory provided 

\k\ = \kYM\+N. 

Let us define A = ^ as we have done in this paper, and |Ay^/| = j^^^- It follows 
that 

lA I 



1-|A| 



and 





AyM 




1 + 


AyM 



In particular as |A| — 1, we have |AyM| oo. Now a Chern-Simons theory 
regulated by the addition of a small Yang Mills term clearly exists at all values of 
I AyM I- However it follows from the discussion above that this only requires the 
dimensionally regulated theory to exist for |A| < 1. Moreover the limit |A| — )■ 1 
should be interpreted as the approach to strong coupling. Our exact result for 
the finite temperature free energy is perfectly consistent with this interpretation. 
The theory ceases to exist for |A| > 1. Moreover the limit |A| — )■ 1 displays the 
extreme thinning of degrees of freedom that is plausible in an extreme strong 
coupling limit 
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4. In Appendix F.3 we have explicitly computed the two loop anomalous dimension, 
at first sub leading order in of the scalar operator ipip. The result so obtained 
agrees perfectly with the result of the same computation performed in Feynman 
gauge. 

The points listed above make it at least plausible that our gauge choice is free of 
problems, and that our regularization scheme preserves gauge invariance and so defines 
a Lorentz invariant theory. While the direct evidence for these claims is substantial, it 
is not overwhelming. Significant additional evidence for the Lorentz invariance of our 
final theory could be obtained by demonstrating the Lorentz invariance of four Fermi 
scattering at one loop. We will not, however, attempt that consistency check in this 
paper. 



13^ 



Note that (2.95) expressed in terms of kyM takes the form 



67r V I^fmI 

This behavior in the A — > 1 hmit could be consistent with a weakly coupled boson description (see [7j). 
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3 Operator Spectrum 



In this section, we study the spectrum of gauge invariant operators in our theory. We 
focus mainly on operators whose dimension is held fixed as N is taken to infinity. 
By large N factorization, all such operators are products of "single trace" operators. 
By a single trace we mean an operator that is constructed from the contraction of a 
fundamental (derivative of ip) and an antifundamental (derivative of ip). We will argue 
below that the spectrum of single trace operators in our theory includes one scalar and 
one current of spin s for s = 1 . . . oo. The dimension of the scalar operator is given, 
at every value of A, by 2 + 0{l/N), while the dimension of the spin s current J^*-* is 
given by A = s + 1 + 0{l/N). We will also argue that the currents J'-'*-' are "almost" 
conserved; more precisely that these currents obey the anomalous conservation law 
listed schematically in (3.16). These nonlinear anomalous conservation equations carry 
a great deal of information, and may eventually prove very useful in "solving" the 
theory. 



3.1 The free limit 

In this subsection we review several well known properties of Chern-Simons theory 
coupled to fundamental fermions in the free (A 0) limit. In this limit our theory 
reduces to a theory of free fundamental fermions subject to a U{N) singlet Gauss 
law constraint; this limit has, of course, been studied in the previous literature (see 
e.g. [2^). In this section we gather those results that will help us study the theory at 
finite A. 



3.1.1 "Single trace" conformal representation content of the free theory 

In this subsection we compute the "single trace" operator content of the theory of 
free fermions. By "single traces" we mean those operators in our theory that are given 
by contracting a single (but arbitrary) fundamental field with an anti-fundamental 
field. The results of this subsection were already presented in [28j ; however we rederive 
them here for completeness. 

Let A represent the scaling dimension and J3 the z component of the angular 
momentum of any operator. We will now compute the partition function 

Trx^/i-^3 (3.1) 
^■^This subsection was worked out in collaboration with J. Bhattacharya. 
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over all single trace operators in our theory. As "single traces" are obtained by mul- 
tiplying an arbitrary fundamental field with an arbitrary anti-fundamental, the trace 
over "single traces" is simply the product of the partition function for elementary 
fundamental fields with the partition function for elementary antifundamental fields. 

The partition function over elementary fundamental fermionic fields (including ar- 
bitrary numbers of derivatives but modulo equations of motion) is simply the character 



of the (1, 2) representation of the conformal group (see Appendix D) and is given by 



(3.2) 



^1 — fix)(l — fi~^x) 

The partition function over antifundamentals is given by the same formula, and so the 
"single trace" partition function is given by Fp. 

It is not difficult to decompose this partition function into the contribution of 
primary operators and their descendents. In order to accomplish this we note that 

-\ 2 



X{fl2 + 



(1 — /ix)(l 



/i- 
1 



'X 



x\l - x)(2 + /i + /i"^) 
(1 — /ix)(l — jJi'^x) 



X 



(3.3) 



(1 — /xx)(l — /i^ix)(l — x) 

1 

(1 — /xx)(l — /i~^x)(l — x) 

00 

= X2fi{x, ^) + X] ^s+i,.(a;, ^) 

s=l 

where XA,s(2;,/i) is the character of a representation of the conformal algebra with 
dimension A and spin s and we have used the character formulae of Appendix |D} 
It follows that single trace operators are given by primaries that transform in the 
representations 

00 

(2,0) + ^(j + l,j) 
i=i 

together with their descendents. It will be important below that the only long operator 



that appears in this decomposition is (2,0) (see Appendix D); all other operators in 
the list above appear in short representations of the conformal algebra. 



3.1.2 Explicit form of the primary operators 

As we have explained above, the primary field content of the free fermion theory is 
given by the (2, 0) operator ipip plus symmetric traceless currents jjif...^^ for all s. As 
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we have seen above, the currents are primaries that head short representations of the 
conformal algebra; the shortening condition is simply the statement that the currents 
obey the conservation equation 

d'Jl:l....... = (3.4) 

Single trace currents J^^^ that have dimension s + 1, spin s and are conserved are unique 
up to a choice of scale in the free fermion theory. We have found explicit expressions 
for each of the currents J'-*-'. Following 31 we find it convenient to package these 
expressions in the form of a generating function 0{x, e) defined by 

0{x;e) = J24:L..,.^''---^'' (3-5) 
where is an arbitrary vector. As all currents are bilinear in the fermions, the gener- 
ating function is given by an expression of the form 

0(x;e)=V^F(7,a^,^,6> (3.6) 

Making a convenient choice for the overall scale of each J*-^^ we find 

F = 7.e7(9^,5,e-) (3.7) 

where 



^ exp (u ■ e — V ■ e) sinh \/2u ■ ve- e — Au ■ ev ■ e 

/(«' e) = /r,^ ^ ^ ^ (3.8) 

y lu ■ ve ■ e — Au ■ ev ■ e 

The Taylor expansion 

f = l + e(u-v) + -e^ (3u^ - lOuv + 3v^ + 2w) 
6 

2, ( 7u'^v luv"^ uw v^ vw\ 

V"6" 6~ ^ ~6~ ^ IT ~ 6" ~ ITy 

+— (10(m - vf{2w - Auv) + {Auv - 2wf + h{u - vf) + O(e^) 

(above, w = u ■ v, u = u ■ e, v = v ■ e.) implies the following explicit expressions for the 
first four currents 

+2(£9^)S^r/^3^, - 2(£9^)9^r7^3^4)V^ 

where all indices above are understood to be symmetrized. 

The existence of spin s conserved currents gives rise to a large symmetry algebra, 
called the higher spin algebra, of the free fermion theory. 
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3.1.3 The full multitrace partition function 



In this section we present a computation of the partition function of the free theory on 



S'^ X S^, i.e. the free theory on 5*^ at finite temperature, along the lines of 39 . By the 



state operator map this is the same as the partition function over all operators (single 
and multi trace) of the free theory. 

The exact partition function on 5*^ x S*^ is given by |4l|Eo] 



Z = J DUexp i-Yl 



iV2|p„|2 _ iV(-l)"+iF^(a;")(p„ + p*) 



n 



(3.9) 



where x = e ^ = e t ^ F{x) is the letter partition function of the fundamental fermions, 
U is a. unitary matrix (the holonomy of around the time circle) and 



Pn 



N 



TrU"" 



Let us first study (3.9) in the limit that x is held fixed as we take the large limit. 



In this case the integral (3.9) is very simply evaluated by completing the square in the 
Pn variables and yields 



(3.10) 



In the high temperature limit Ff{x) = 2T^ and (3.10) reduces approximately to 

InZ = 4T\(5) 



(3.11) 



(3.10 ) is simply the Bose exponentiation the single particle partition function Ff{xY . 
However we have already argued above that F{xY is precisely the single trace oper- 



ator partition function. It follows that (3.10) represents the partition function of a 



noninteracting gas of single trace operators (or particles), with single trace partition 
function Fp(x)^. 



In order to obtain the result (3.10) we shifted the variable p„ by ^''Z ^ ■ When x is 



N 



of order unity this shift is of order and is negligible. In the high temperature limit, 
however, 

Fpix) ^ 2T^ 



and the large saddle point value of the integral (3.10) is given by 

2T2 



Pn 



-1) 



n+1 



If we set 



T 



Nt 
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and hold t fixed as is sent to oo, the saddle point value of p„ is independent of 



Pn 



-1 



(3.12) 



Let p{9) denote the density of eigenvalues of the unitary matrix U. It follows that 

1 



Pn 



dep{e)e 



inO 



The saddle point eigenvalue distribution may be reconstructed from its Fourier modes 
via 



piO) 



2lT 



') 



n=l 



The eigenvalue distribution that corresponds, in particular, to (3.12) is given by 

2rrp{e) 



1 + ^ — cos ne. 



n=l 



This eigenvalue distribution takes its minimum value at 6 

V3 



t > 



TT. p(7r) is negative when 
(3.13) 



Now a negative value for the eigenvalue density function is meaningless; when the 
inequality (3.13) is obeyed what really happens is that our system undergoes a Gross- 



Witten-Wadia phase transition [29||30] to a gapped phase. 

The saddle point solution for p{6) is quantitatively complicated above the phase 
transition temperature (3.13). In broad qualitative terms, however, the eigenvalue 
distribution gets increasingly peaked as t increases. In the limit t ^ 1 the eigenvalue 
distribution tends to a 5 function. In this limit pn = ^ for all n, and the logarithm of 
the partition function is given by 

InZ = ANT^ ^ — = 3iVC(3)T2 (3.14) 

n=l " 

In this limit we recover the partition function of the free theory on (^^i^e leading 
term in ( 2.94[ )), with the volume factor V2 replaced by 47r, the volume of the unit two 
sphere. 

In summary, the partition function of our theory on an of unit radius interpolates 
between the partition function of a gas over single trace operators (3.10) and (3.11) 



(when T -C V N) to the gas of deconfined fermions (when t ^ y N) . Note that when 



hiZ is of order A^^ when T approaches yN both from below (see (3.11)) and from 
above (see (3.14)). The interpolation between these behaviours is not smooth; it goes 
through a third order Gross- Witten-Wadia transition. 
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3.2 Non-renormalization of the scaling dimension of the cur- 
rent operators 



As we have explained above, the single trace operator content, in the free limit, is given 
in terms of representations of the conformal algebra by 



oo 



(0,2) + }_^(j + l,j). 



As we have explained in Appendix D, representations of the sort (s + 1, s) are short. 
Operators that transform in these representations can develop anomalous dimensions 
only after combining with a long representation with quantum numbers (s + 2, s — 1) 



(see (|D3|)). 

Now the operators in the (2, 1) and (3, 2) conformal representations are the con- 
served currents for the U{1) fermion flavour symmetry and the fermionic stress tensor 
respectively. These operators are exactly conserved currents so cannot develop anoma- 
lous dimensions at any value of A. Let us now consider the operator that transforms in 
the (4, 3) representation of the conformal algebra in the free theory. This operator can 
develop an anomalous dimension only upon combining with an operator in the rep- 
resentation (5,2). However the only spin two single trace operator, the stress tensor, 
transforms as (3,2) for all values of A. 

We now make a key assumption that we justify in much more detail below: we 
assume that, as far as the analysis of leading large N scaling dimensions is concerned, 
we can simply ignore all mixing of single-trace and multi-trace operators. It follows 
that, at leading order in the operator in the (4,3) representation cannot develop 
an anomalous dimension at any value of A, simply because at no value of A can there 
exist a single-trace operator in the (5, 2) representation with which the (4, 3) operator 
can combine to form a long representation. 

This argument can now be repeated recursively, to demonstrate that the scaling 
dimensions of all "single trace" spin s operators are exactly protected (at leading order 
in A^) at all values of A. The only way this argument could fail is if the theory underwent 
a severe phase transition at a finite value of A, across which the spectrum of the theory 
was not continuous. 

Note that the non-renormalization theorem relies on the assumption of non mixing 
of single and multi trace operators (which we will justify in detail below) and the the 
sparsity of single trace operators in our system. As we will see explicitly below, our 
non-renormalization theorem will be violated by corrections. 

Finally note that nothing in the argument we have presented so far prevents the 
operator ipip from developing an anomalous dimension. However we will proceed to 
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argue below that this is also impossible; the scaling dimension of this operator is also 
protected ( function of A) in the interacting theory. 



3.3 Explicit form of the current operators 



In subsection 3.1.2 above we have determined the explicit form of the primary operators 
that transform in the + s) representation of the conformal algebra. In the previous 
subsection we have argued that the interacting theory has primaries with the same 
quantum numbers at all values of A. In this subsection we will determine the explicit 
form of these primary operators, in the interacting theory, in terms of the bare fermionic 
fields 

In the interacting theory at finite A let J^^^ denote the currents obtained from the 
following procedure: replace all derivatives in the expression for J^^^ in the free theory 



(see subsection 3.1.2) with covariant derivatives. Explicitly, for s = 1 ... 4 we have 



where all indices on the RHS are understood to be symmetrized. 

The currents J^^^ are not automatically tracelessness in their vector indices. The 
reason for this lack of tracelessness is that covariant derivatives, unlike their ordinary 
counterparts, do not commute. However the commutator of two covariant derivatives 
is a field strength. Now according to the classical Chern-Simons equation of motion 

(F,,)] = p^upi^'Y^j (3.15) 

where the i and j indices are colour indices. Now consider a factor of the field strength 
F inserted inside a "single trace" fermion bilinear. By the equation of motion cited 



above this insertion splits the "single trace" into a "double trace" operator divided by 
k. Further factors of F inside any of the new resultant "traces" repeats this operation. 
It follows that the spin s — 2, s — 4 etc components of the current J^^^ is given (via 



While the equation (3.15 1 was derived classically, we believe it also applies quantum mechanically, 
in an appropriate regulator scheme, as the current J^^^ is the unique dimension two, spin one operator 
in the theory. 
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the equations of motion) by 'multi trace' operators that are schematically take the 
form -^^^ where s represents any "single trace" operator so stands for an m trace 
operator. 

The fact that the operators J*-*^ are not traceless means that these currents are 
not of definite spin; they include components of spin s, spin s — 2 etc. Let us define 
the interacting currents J*^*) as the projection of J*^*-* to its spin s component, i.e. the 
projection that removes all traces from J*^*^ is, by definition, a spin s current, and 
is of power counting dimension s + 1. As we have explained above, J'-*'-' and J*^*-' difi'er 
only by "multi trace" expressions. 

Now the primary operator that transforms in the (s + 1, s) representation of the 
conformal algebra is necessarily an expression of spin s and of power counting dimension 
s + 1 (we use a renormalization scheme in which operators can mix only if they have 
equal classical dimension). The full set of such operators is easily enumerated; in 
the free theory it consists of those descendents of the primaries {j + that are of 
dimension s + 1 and spin s (clearly this requires j < s). Note that (in the free theory) 
these are all single trace operators. All spin s multitrace operators in the free theory 
have dimension greater than or equal to s + 2. 

In the interacting theory, the full set of operators of spin s and classical dimension 
s + 1, is given by replacing all derivatives by covariant derivatives in the free answer 
and then projecting onto the spin s component (i.e. removing all traces). As above, 
the projection onto spin s leaves the "single trace" part of the operator untouched, but 
adds "multi trace" operators to the mix. 

Now let us compute the divergence of the most general possible spin s, classical 
dimension s + 1 current listed above in the interacting theory. The computation of 
this divergence differs from the same computation in the free theory in three ways. 
First covariant derivatives do not commute, and that results in extra factors of the 
field strength; as we have explained above such factors modify the multi trace parts 
of the answer but leave the "single trace" part of the answer untouched. Second, as 
explained above, the current J*-*) has extra multitrace operators as compared to the 
free current. This additional complication also affects only the multi trace parts of 
the answer. Finally the fermion equation of motion could be quantum mechanically 
modified, but any such modification is necessarily in terms of 'multi trace' operators. 

In other words the single trace part of the divergence of the general spin s current 
in the interacting theory is identical to the result of the same computation in the free 
theory (after replacing derivatives with covariant derivatives). However the interacting 
divergence includes, in addition, several multi trace contributions that are absent in 
the free theory. 
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Now in subsection 3.1.2 we computed the unique operator of dimension s + 1 and 
spin s, in the free theory, that is also conserved. The interacting theory possesses no 
exactly conserved operator of this dimension. The operator J^^'^ comes closest to a 
conserved current, in that it is the unique current that obeys the schematic equation 

d- J^''^ ^yJJ+ ^JJJ (3.16) 
k k'^ 

(on the RHS of the equation above the symbol J refers either to a current or to a 
descendent of a current; the important point in this equation is that the RHS contains 
no single trace pieces). In other words J^^-* is the unique spin s and classical 
dimension s + 1 field in the interacting theory whose divergence has no single trace 
component. 

As we will see below, the operator J^^^ constructed above, may be identified with 
unique dimension s + 1 spin s primary of the interacting theory. Before proceeding we 
will first, however, give an example of how our abstract and exact construction of J^'^\ 
(as the projector onto the spin s sector of the current J*) may actually be practically 
implemented at leading order in A, for the special case of the current J^^\ For this 
purpose in the next subsection we evaluate the current J^^^ and its divergence using 
the classical (but interacting) field equations. In the subsequent subsection we return 
to the demonstration that the dimension of J^*-* is s + 1. 



7(3) <^ _ 1 



3.3.1 J^'^^ as an example 

The trace of J^^^ is given by 

Using the identities hsted in Appendix [E] 

Jj^j" = -^V^i^^'^e,..(7^7M + 7,7^)^ 

= ^(^^)(^7'^^)- 

In the last line we used the equation of motion for the field strength F (see Appendix 
|E]). Upon projecting out the trace we find 



^^The absence of higher trace operators on the RHS of this equation foUows from a consideration 
of quantum numbers. The is of spin s — 1 and classical dimension s + 2. Recall that A — s > 1 for all 
single trace operators. It follows that A — s > m for m trace operators, and so m < 3. 
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As we have explained, the currents J^^-* are not expected to be exactly divergence 
free for s > 3. On the other hand the divergence of the currents J*-^-* and J^^-* vanishes 
even in the interacting theory (this is obvious for J^^^ and is also true for the stress 
tensor T^'^ as it can be explicitly checked). In Appendix |E] we have also explicitly 
computed the divergence of J^^-* in the classical but interacting fermion theory. Our 
result is (see Appendix E.3[ ) 

IGvr 



f)l^ 7(3) 



5k 



(3.19) 



We have presented our result in terms of the scalar "current" J^^^ = ifjip and the vector 



current J, 



(1) 



The equation (3.19) has been derived classically; it could certainly receive quantum 
corrections. However all such corrections are necessarily of higher order in A; at small 
A (3.19) gives the leading order contribution to the divergence of J(3). 



3.4 Anomalous dimensions of the current operators 



We will now use the fact that the operators J^^^ obey equations of the form (3.16) to 



argue that the scaling dimension of the operator J'-*-' is s + 1 up to corrections of order 
j^. We will also explain how the knowledge of the precise form of the RHS of (3.16) 



can immediately be converted into a computation of the anomalous dimensions of J^^-* 
at leading order in and systematically in the expansion. 

Let us first review how it follows that a conserved spin s current has A = s + 1 . For 
this purpose we use the state operator map, and denote the state dual to the operator 
O by \0). The standard argument takes the schematic form 

{dJ\dJ) = {J\[K,P]\J) = {A-s- 1){J\\J) 

(see below for the argument with all indices in place). If {dJ\dJ) vanishes then it 
follows immediately that A = s + 1. In the our theory {dJ\dJ) does not quite vanish 



(see (3.16)). However as we will argue below, (3.16) determines it to be a times a 
product of known two point functions. 

In order to all see this in detail we start with some preliminaries. In radial quan- 
tization, jI^"* inserted at the origin corresponds to a state on the sphere. For 

^ ^ is) 

convenience we introduce a set of differently normalized currents, jTi , whose corre- 
sponding states obey 

{3^^\3^u^) = 5ss'5,,._ (3.20) 
where 5^^y is 1 if /i, z/ are the same set of indices up to permutation and otherwise. J'-*-' 

( \^ is) is) 

and i are related by J^j = > where the tree level result for the normalization 
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constant is computed in Appendix B.2.2. Using the inversion x'^ = x^/x^, the two 
point function in position space is related by 



dx" 



dx" 



-)V 



-2s-2~2S, 



1=1 



x^ 



(3.21) 



where 5s is a possible anomalous dimension for J^^\ In the second step we moved 
j^\x') to j^\o) at no cost, since the difference is a conformal descendant, which is 
orthogonal to |j(0)). One can analogously work out the momentum space two point 
function. 



We can also translate (3.16) into the language of states in radial quantization, 
schematically of the form 



Taking its norm (using that = K)^ we have 



(3.22) 



(3.23) 



Since is a conformal primary, it is annihilated by K^. Using the conformal 

algebra, the LHS of (3.23) is equal to 



2 



s-l 



1 = 1 



(3.24) 



= 28 ) 

In the last step we used the fact that j^i^- -^^ is traceless. The order 1/A^ contribution to 
the RHS of (3.23) comes from the disconnected four-point function, i.e. schematically 

1 . 



{jj\jj) = mm + o{ 



N' 



(3.25) 



This relates the anomalous dimension 6s at order 1/A^ to the product of two two-point 
functions of J of lower spins in the free theory. Note that this method cannot be used 
in the s = case. In Appendix E, an explicit two-loop calculation for the 0{1/N) 
anomalous dimension of the scalar operator ipilj is performed. 
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(3) 

As an explicit example, let us consider the spin 3 current Jij.Jp, which obeys the 
anomalous current conservation relation (3.19). In terms of states in radial quantiza- 
tion, we have the (tree level) relation 

(3.26)' 

where c = Taking its norm and using the conformal algebra commutation 

relations, we find, in particular. 



i^t^ul^fju) = 252|c|2 



20^3, 



(3.27) 



and so 



63, 



252 A2 
625 iV^ 



(3.28) 



to the leading nontrivial order, i.e. two-loop order. There are both planar and non- 
planar corrections at higher orders in A and in 



A relation of the form (3.23) holds to higher order in as well; we would however 
need the connected four-point functions etc., as well as potential corrections to 
the relation (3.16), in order to compute the next order contribution to 6s in 1/A^. 



3.5 Anomalous current conservation relations within correla- 
tion functions 

As we have explained above, the primary operators that transform in the (s -|- 1, s) 
representation (in the large limit) obey anomalous conservation equations of the form 
(3.16). The equation obeyed, in particular, by J^^^ (in the classical interacting theory) 



is listed in (3.19). The nonlinear equations of form (3.19) carry a lot of information 



as we will explore in this subsection. 

Let us first investigate the implication of anomalous conservation equations on three 
point functions of spin J^^^ operators at leading order in A^. We find the schematic 
equation 

d^U. J J) = I (J J J J) + ^UJJJJ) (3.29) 
k 

While the leading behaviour of the second term in ( |3.29[ ) (which comes from factorizing 
the 5 pt function into the product of a 2 and 3 pt function) is (9(1) (and so subleading in 
the large A^ limit) the leading behaviour of the first term on the RHS (this comes from 



^^We have seen in the previous section that the operators J'-'*^ are effectively conserved within two 
point functions at leading order in N 
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factoring the 4 pt function into a product of two point functions) is 0{N) and so of 
leading order. It follows that our current operators are not, in general, conserved within 
three point functions even at leading order in the 0{1/N) expansion. By equating 



scaling dimensions on the LHS and RHS of (3.29) it follows immediately that the RHS 



of (3.29 ) can be nonzero only when s > Si + S2- In other words the three point function 

{j(') jM) (3.30) 

can obey anomalous conservation equations (rather than genuine conservation equa- 
tions) only if Si, S2 and S3 violate or saturate the triangle inequality. In all the explicit 



computations we perform below, it will turn out that the RHS of (3.29) is nonzero only 
when the triangle inequality is explicitly violated (i.e. it vanishes when si = S2 + S3). 
It is possible that this is a general exact result that could perhaps be proved by an 
analysis of allowed structures for 3 point functions, but we will not attempt such an 
analysis in this paper. 



3.6 Non-renormalization of the scalar operator J^^^ 

Let us now apply the arguments of the previous subsection to the three point correlator 
{J^^\x)J^^\y)J^^\z)). According to the arguments of the previous subsection, and 



(3.19), the divergence of this three point function w.r.t the variable x is schematically 
proportional (in the large N limit) to a term proportional ((9J*^°^ J*^°'')( J'-^^ J*^^-*) plus 
another term proportional to {J^^^ J^^^) {dJ^^^ J^^^) . The weight under overall scaling of 
each of these expressions is 5 + 2Ao where Aq is the as yet unknown scaling dimension of 
the current operator J*-*^) . On the other hand the weight under scaling of the divergence 
of the three pt function is 6 + Aq. Equating these weights we find that Aq = 2. We 
conclude that, just like their spin s counterparts, the scaling dimension of the scalar 
current J^'^^ is not renormalized clS 8b function of A I I 

This non-renormalization result may be argued for more intuitively as follows: the 



LHS of (3.19) is of dimension 5. The RHS of the same expression is also of dimension 
5 if and only J^^^ has dimension 2. The slightly technical argument of the previous 
paragraph simply makes the argument precise. 
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We thank O. Aharony for discussions on this topic. 
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4 Three Point Functions 



4.1 Quantum corrected operators and their 3-point functions 



In the previous section we demonstrated that the dimension s + 1 spin s operators 
primary operators jjii...f_is given by the following construction. First replace all 



derivatives in (3.5), (3.7) and (3.8) by covariant derivatives. Then project the resultant 
currents, J*, onto the spin s subsector (i.e. we project all lower spins out of the 
currents). While this procedure is formally exact, a practical implementation of the 
projection procedure requires use of equations of motion to express F^^, in terms of 
fermion bilinears. This replacement could receive quantum corrections; moreover it 
gives rise to quite complicated expressions involving a summation over multi trace 
contributions. 

All complications having to do with the projection operation disappear when we 
focus on the particular component of J*^'^^ that is obtained by dotting all of its free 
indices with any constant null vector e. This is because all multitrace contributions to 
J^^^ come multiplied with factors of the metric with free indices and vanish when all 
these indices are contracted with e (we use (e)^ = ). Consequently, in this subsection, 
we will work exclusively with the operators obtained by contracting all free indices of 
J^^^ with a constant null vector. 

In the gauge employed in this paper there is an additional simplification if the null 
vector in question is chosen to lie in the x~ direction. In that case the expression 
for J^'^^ involves only D_ derivatives of and xp. In the lightcone gauge employed 
throughout this paper and in much of this subsection, y4_ vanishes, so that D_ = d-. 
It follows that this particular component of the spin s current is given directly by the 
expressions (3.5), (3.7) and (3.8) upon setting e = e~. Exphcitly we have 



(4.1) 



where 



-(2q+-p+) 



sin {2e \/ q'^ (p+ — 



2e-y/q+ip+-q^ 



s>l. 



(4.2) 



(see (3.8)). In this particular case, it is easy to define and compute J'^^\p\e), the exact 
quantum corrected version of the operator j^^\p\£)- We define 



J 



cP'q 



{P\e) = J j^,^'''{p-q)ji'^{q;p\eW{q) 



(4.3) 
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where J'aJiQ'iPl^) is defined to be the sum over all the graphs depicted in Fig. 2. 
Clearly, J^^p [q] p\e) obeys the Schwinger-Dyson equation 



(2^ 



r 



07 



1 



r7 



where jl^l{q]p\s) is the classical spin-s vertex. 



(4.4) 



G^,{l-q). 
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In the special case of s = 1, it is not difficult to define, and set up an integral 



equation (analogous to (4.4)) for the exact quantum corrected operator for arbitrary 



polarization vector (spin 1 is special because the expression for the current contains no 
derivatives). 



Once we have solved the integral equation (4.4), three point functions of the corre- 



sponding operators are determined simply by Wick contraction using the exact fermion 
propagator. The exact expression for two point function (O1O2) may also be obtained 
by Wick contracting the exact operator Oi with the bare operator O2 (or vica versa). 



We postpone all discussion of the solutions to the integral equations (4.4), to future 
work 




+ 




+ 



Figure 2: The quantum corrected vertex. 



4.2 Three point functions in perturbation theory 

In the rest of this section, we postpone the attempt to compute correlators exactly in 
A (at leading order in large N) but instead turn to study of three point functions of 
fermion bilinear operators/currents J'-*-* perturbatively in A. 

We will begin by writing down the three point functions of currents in the free 
fermion theory, and describe some general structures of three point functions of higher 
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spin currents as constrained by conformal symmetry. While the full tensor structure of 
the three point functions {J J J) of currents of nonzero spin is somewhat complicated, 
it is constrained by conformal symmetry to be the linear combination of finitely many 



independent structures. In fact, as shown in 36 , for conserved higher spin currents 
all parity preserving structures of {j'-'^U^'^U^'^^) can be realized by free scalars and 
by free fermions, while a parity odd structure of (J^'^i^ J*^''^) j(«3)^ exists when the three 
spins si,S2,S3 obey triangular inequality]^ The latter cannot be realized with free 
fields. One of the main results of this section is that the parity odd structure is indeed 
present in the interacting Chern-Simons-fermion theory. 

The parity symmetry of Chern-Simons-fermion theory may be restored by simulta- 
neously taking k — > —k. Therefore, in perturbation theory, the parity invariant tensor 
structures of {J J J) occur at even loop order, while the parity odd tensor structures 
occur at odd loop order. We will analyze the one-loop contribution to the parity odd 
structure of {J^^^ J^^^ J^'^^) and {J^'^\j^'^\j^'^^) , and the two-loop contribution to the par- 
ity even structure of ( J*^*^) J'^'^^-' J^*^)) for general nonzero Si, 52,^3. The results will be 
summarized in section 4.3, while the details of the computation are left to Appendix 
G. These results will be very useful in comparing with potentially holographically dual 
higher spin gauge theories in AdS4^ (as will be discussed in section 5). 



4.3 Three point functions in the free fermion theory 

As seen before, a generating function for the conserved higher spin currents in the free 
fermion theory is 

J(x, e) = V'(x)7 ■ ef{e ■ d ,e ■ d )ip{x) (4.5) 

where 

= e-ii^ (4.6) 

2Juv 



and e is a null vector. By expanding to terms, this gives the spin-s currents for 
all nonzero s. The scalar operator iljip need to be considered separately. It is useful 
to write the null polarization vectors in a bispinor form, = Aq,A^. Since and f) 
annihilate ip by the free equation of motion, we can treat (9^ as a null vector as well 
when acting on ip, and make the substitutions 

where ^a,Va are formal spinorial symbols. We can write the generating current as 

^^Recall that when the spins si, 82,83 obey the triangular inequaUty, the non-conservation of the 
currents in our theory does not show up in 3-point functions at leading order at large N. 
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The equation of motion can be implemented by r/?/' = and = 0. We may then 
make a further replacement 

Ipa Va(p, i'a (4.9) 

where is a formal object which behaves as a free massless scalar field, while keeping 
in mind that an extra minus sign comes with a fermion loop. We can then write 

(4.10) 

As a simple check, the divergence of the current is obtained from 

(dxmJi^, A) ^ [i^d.y + ivdx)'] J(x, A) = ((^ + tv)dxm " ^v)dx)J{x, A) (4.11) 

which obviously vanishes. 

Now let us compute the three point function {J J J) for free fermions. The fermion 
propagator is 



1 (^-^%/3 

47r \x — x'P 



X 



Here we are making the identification C,'g ~ irjis when acting on (0(x)0(x')), so as to 

^ — 
be consistent with d'^^ = ^'^^'^ ~ —VaVp = — d ai3- The three point functions for the 

generating current, {J{xi, Ai) J(x2, A2) ^(^3, A3)), can be then written as 



1 



En 



1 



+ (1^2) 



Ti exp 



Ti=±l 1 = 1 



-Si ■ ^ i + nV Si ■ 



(4.13) 



1 



l-^iji+l I 



+ (1^2), 



where di and act cyclically along the product. We can further rescale Ei — )■ Uei, 
expand this generating function in ti, and obtain up to normalization the three point 
functions of three currents of nonzero spins 

^ ] ^1 ^2 t3^{Jsi{Xi,ei)Js^{x2,S2)Jsa {X3,e3)). 



(4.14) 



It follows from essentially the same proof as in the appendix of 32 that an alternative 
generating function can be written in a simple closed form 



Sl,S2,S3 

2 



sinh sinh Pi sinh P2 sinh P3. 



(4.15) 



(47r)3 Ixisllxssllxail 



The latter is naturally reproduced from the tree level three point function in the bulk 



dual type B Vasiliev theory, as shown in 32 
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4.4 A few comments on the general structure of three point 
functions 



Let us begin with an example: the three point function of the stress energy tensor 
T^jy, a quantity of interest in all CFTs. There are three linearly independent tensor 
structures for {Tf^-^,^j^{xi)Tfj,2^^{x2)T^.^i^^{x3)) that are allowed by conformal symmetry 
(together with current conservation) [34], [35], [36] . In other words. 



(TTT) = cb{TTT)b + cf{TTT)f + Com{TTT)om. (4.16) 

Here {TTT) b stands for the three point function of the stress energy tensor of a free 
massless scalar, (TTT) p that of a free massless fermion, both of which are parity 
preserving; {TTT) odd is a parity odd structure that is present only in parity violating 
interacting CFTs. 

When the fermions are coupled to Chern-Simons gauge fields, both (TTT) b and 
{TTT) odd show up in the three point function of T^^. As we will see later, {TTT) odd 
occurs already at one-loop order in perturbation theory, while {TTT) b occurs at two- 
loop order (together with the free fermion structure {TTT)f). The relevant Feynman 
integrals are difficult to perform directly for generic configuration of operator insertions 
and polarizations. Instead, we make use of the conformal symmetry, and examine 
special limits in the configuration space where only a single tensor structure of interest 
survives; then all we need to do is to compute the overall coefficient of the surviving 
tensor structure, which is manageable. 

We will make use of both conformal symmetry and current conservation to con- 
strain the three point function {Jsj^{xi,ei)Js2{x2,e2)Js3{x3,e3)). At leading order in 
1/A^, this assumption is valid, up to the anomalous current conservation relation, of 
the form (9 ■ J ~ ^JJ -l- ■ ■ ■ , discussed in section 2. The relevant part of the divergence 
of the higher spin currents is expressed in terms of double trace operators, namely 
product of two currents, of lower spins. This leads to violating of current conser- 
vation in the three point function, but only is one of the three spins, {31,82,33), is 
greater than the sum of the other two, in other words {31,32, S3) violating the triangu- 
lar inequality. Also note that this divergence of the three point function is necessarily 
parity odd (due to the 1/k factor in the anomalous current conservation relation). So 
for the parity even structures, and for parity odd structures with (si, 52,53) obeying 
triangular inequality, we may proceed by assuming current conservation. It is shown 
in 36 under this assumption that {Jsi{xi,ei)Js2{x2,£2)Js3{x3,e3)) must be of the form 



13^121 ^|3;23| ^|3;3i| ^ times a polynomial of appropriate degrees in the conformally in- 
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variant objects Pi,Qi,Si. The latter are defined as 



Pi 



P9 = 



P3 = 



23 



X 



31 



X 



12 



Qi = -Ai 



a; 



12 



X 



31 



23 



+ 



+ g2 = -A2(^ + %)A2, g3 = -A:^^='^ ' 

3 ^12 

+ (A2A3)A2(^ + 1^)A3 



31 



X 



)A3, 



23 



^1 = 


i 


kl2| 




4 


2^23 2:31 


"''12"''23 


^2 = 


i 


^23! 


(A3^23i^31^3)(Ali^3lAi) 


4 


|a:3i|ki2| 


2 2 
■^23-^31 


^3 = 


i 


k3l| 


'(Ai^31^2^l)(^2A2^2) 


4 


kl2|k23| 


^31-^12 



'?^12 I ^ 



X 



12 



X 



23 



+ (AaAOAad^ + ^)A, 



' /y»2 /y>2 

•^23 -^31 



+ (AiA2)Ai 



7^31 , ^ 



X 



+ 



12 



31 



X 



)A2 



12 



(4.17) 

where again we wrote the null polarization vector in the bispinor form, in terms of 
a spinor Aj. The polynomial in P,Q,S is such that it is homogeneous of degree 2s j in 
Aj. Because all quadratic monomials in S'j's can be expressed in terms of Pi, Qi, we can 
assume the polynomial to be linear in the Si^s. 

Note that under the exchange {xi, Ai) •<->■ {x2, X2), 



Pi ^ -P2, P2 ^ -Pi, P3 ^ -P3, 

Qi -Q2, Q2 -Qi, Q3 -Qs- 



(4.18) 



In case of correlators involving currents of identical spins, the symmetry under exchange 
of those current should be imposed correspondingly. 

The L-loop contribution should transform under parity with the sign (— )^, namely 

{Jsi{—Xi, —£l)Js2i~^2, —£2)Js3i~X3, ~£^3))l-1oop = {Jsii^l, £l)Js2i^2, £2)Js3{^3, ^s)) L-loop 

(4.19) 

for si, S2, S3 > 0. The spin case is special because Jq is parity odd unlike Jg for s > 0. 
We have therefore 

{Js^{-Xi, -ei)Js^{-X2, -e2)Jo{-X3))L = i-)^^^{JsA^l^^l)-Js2{x2,£2)M^3))L, 
{Jsii-Xl,-Si)Jo{-X2)M-X3))L = i-)^{Jsiixi,ei)Jo{x2)Jo{x3))L, 
{Jo{-Xi)Jo{-X2)Jo{-X3))l = {-)^^^{Mxi)Mx2)Jo{x3))l. 

(4.20) 

For example, {JqJqJq) can only receive contribution from odd loops. Under x — > 
—x,e — > —e, Pi,Qi are invariant while Sij changes sign. So we expect the L-loop 
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contribution to the three point function of currents of nonzero spins to be of the form 

{JsAxi,ei)Js^{x2,e2)Jss{x3,£3))L = ] rr^—T, tF{si,S2,S3\P,Q), L even, 

F12|F23|F3i| 

{JsAxi,ei)J,^{x2,e2)Js3{x3,e3))L = ] n~n iy^Si,i+iGi+2{si,S2,S3\P,Q), L odd. 

|a;i2||a;23lF3i| ^ ^ 

(4.21) 

Let us now examine the structure of the three point function of the spin 2 current, 
i.e. the stress energy tensor T^^,, in some detail. F(2, 2, 2|P, Q) contain two possible 
structures which are that of the free scalar and free fermion theories. At odd loop level, 
the contribution to (TTT), if nonzero, must be proportional to 

S^{P?Ql + 5PiP^) + S2{P^Ql + hPjP!) + S3{PIQI + 5PfP|) ^^22) 

|3;i2| |2;23| I^^SlI 

We would like to see whether this structure appears at one loop (at planar level). If 
such a term is present at one-loop, it would be the entire one- loop contribution to 
(TTT) since it is the only tensor structure allowed by conformal symmetry. Therefore, 
it suffices to determine the overall coefficient by taking some limit in which the tensor 
structure simplifies but stays nonzero. 

Similarly, there is a parity violating structure for (jjT), where j is the spin 1 flavor 
current. It is proportional to 

2S,Pi + 2S2P? + S3Ql _ ^^23) 
|3;i2| |a^23| l^^ail 

Let us consider the special conflguration in which all three polarization vectors (or 
spinors) are equal, namely Aj = A (z = 1, 2, 3), so that 

= -P2 - P3, (4.24) 

and 

(4.25) 

4 |a;i2||a;23||a;3i| 

etc., where we defined 

A = A;^3i^i2^ = A;j^i2i^23^ = ^23^1^- (4-26) 



We can write (4.23) as 

I ^ PsiPi + P2 Y + 2P^P 2{Pl + P2 
a2'^23"''31 



A '^^ ' ^ ^> (4.27) 
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Now consider the limit in which xi and X2 are very close, while X3 is kept a finite 
distance away. For instance, xu = S, X31 = x, X32 = x + 6, and x is kept finite. We 
have now 



Pi 



e ■ {x + S) 



P, 



e ■ X 



x^ 



^3 = ^, A = A#A, 



(4.28) 



To simplify things further, we may choose x ■ e = 0, say by taking x = Icq along 
Euclidean time direction (eo is the unit vector in Euclidean time direction) and 5,e 



purely spatial. We can choose e to lie along light cone direction x . In this case, (4.28) 
becomes 



Pi 



6+ 



P2 = 0, A 



The parity odd tensor structure for {js(0)js{—S)Ti;{x)) then reduces to 



t6' 



where we took the 6-^0 limit on the RHS. 



z {6+Y 
4 SH^ 



(4.29) 



(4.30) 



Similarly, we can write the parity odd structure of (TTT) in this special case Aj = A 
(z = 1, 2, 3) as 

Z^ PsiPjiPl + P2? + ^PlPj) + Pl(Pf (A + P,? + ^PlPl) + P2(P|(P3 + Plf + 5P|Pf ) 



2 2 2 
•^12"''23"''31 



A 



A ^2 ^2 ^2 
^ '''12"''23'''31 



(Pi + P2 + P3)(PlP2 + P2P3 + P3P1) 



In the limit of (4.29), this becomes 



I (5^ 



4(58(^2 + 52)4 4 



(4.31) 



(4.32) 



where the RHS is the 5 — )■ limit. Note that if we had taken a generic polarization 
configuration, in the 5 — )■ limit we would expect the three point function of T to scale 
like b~^x~^ . To determine the parity odd structure in our limit, we need to go one 
order higher in 5 than the naive leading order. In the next subsection, we will compute 
the coefficient of this term. 



4.5 One-loop and two-loop results 

In this subsection, we summarize the planar one-loop and two-loop perturbative results 
for the three point functions, leaving details of the computation to Appendix G. 
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First we consider the one-loop contribution to the scalar three point function 
{J^^\xi)J^^\x2)J^^\x3)). Note that J^^^ = ijjip is parity odd, and its three point 
function can only receive contributions at odd loop order, a priori. We find that it in 
fact vanishes at one-loop, as shown in Appendix G.l. The computation is performed 
in both Feynman gauge and lightcone gauge; in Feynman gauge there are nontriv- 
ial cancellations between different diagrams, whereas in lightcone gauge each diagram 
vanishes individually. 

In fact, we suspect that {J^'^\xi)J^'^\x2)J^'^\x3)) vanishes identically at planar 
level. This is suggested by the holographic dual, since J*^°^ is dual to a bulk scalar 
of mass —2 in AdS units, and such a scalar should have vanishing cubic coupling, 
otherwise if one imposes the alternative boundary condition in which the dual scalar 
operator has dimension 1, the tree level three point function would diverge. This is 
in particular the case for the general parity violating version of Vasiliev's theory. In 
field theoretic terms, one may be able to see this by considering the UV fixed point 
obtained from a double trace deformation by J^^\J^^\ and relate {J^^^ J^^\j^'^^) to the 
corresponding scalar three point function in the UV fixed point theory. 

Next, we analyze the one-loop parity odd structure in two examples, (jjT) and 
(TTT), where j = J^^^ is the U{1) flavor current. Since there is only a single parity 
odd tensor structure allowed by conformal symmetry and current conservation, we 
only need to extract the overall coefficient. The computation is mostly easily done in 
temporal gauge, and by taking the special limit described in the previous subsection. 
In both cases, we find the answer to be nonzero. The results, expressed in terms of the 



conformally invariant objects Pi,Qi,Si defined in [36] and in the previous subsection, 
are 



UiXi, €i)j{x2, e2)T{x3, es))i^ioop 



2iX 2SiP^ + 2S2P? + S-sQ 



3 



|a;i2||a;23||a;3i| 



{T{xi,ei)T{x2, e2)T{x3, e3))i-«oop (4.33) 



mX S,{PlQ\ + 5P|P|) + ^2(^1^2 + 5P|^f ) + S^^PlQl + 5Pf P 



|a;i2||a;23||a;3i| 

Finally, we compute the two-loop parity even structure for ( J^'^) J^^^) j(s3)^ ^■■^^■^ 
general spins. As remarked earlier, this three point function is conserved to leading 
order in if the three spins obey triangular inequality, and when the triangular 
inequality is not obeyed, the non-conserving structure is necessarily parity odd and 
does not show up at two-loop. So at two-loop order the result is a linear combination 
of the free scalar and free fermion answer. The computation in Appendix G.3 focuses 
on extracting the coefficients of the free scalar structure {JJJ)b, by working in the 
temporal gauge and taking a special limit in which the free fermion structure {JJJ)p 
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while the free scalar structure survives. Our results may be summarized in the following 
form: 



(4.34) 

where the currents are normalized by their two point functions, with 

c. ^ (4.35) 



A nontrivial feature of this result is that cb is independent of the spins Si, 52,53. As 
will be discussed in section 5, this property agrees with what one expects from the tree 
level three point function in the parity violating version of Vasiliev's higher spin gauge 
theory in AdS4^. The agreement also suggests that Cp^^^^ should be independent of the 
spins as well, though we have not computed this explicitly at two-loop order. 



5 The Formal Solution at Large 
5.1 Exact effective action for i/j and S 

In the large N limit, the exact quantum effective action for the S field, obtained after 



integrating out the fermions, is given by (2.44). The field S referred to above is an 



auxiliary field introduced as a computational device. The effective action that we 
are really interested in, in our problem, is the effective action for ip. This quantum 
effective action contains all information needed to compute scattering of ip quanta, and 
also allows the exact computation of correlators of some operators (see the previous 
subsection). In this subsection we outline how the exact quantum effective action for 
the fields ip cind E may be computed in the large N limit. 



All interactions in the action (2.41) are of the form ipTjip; consequently interactions 
always involve a E field. However the propagator of the E field is (9(-^). It follows that 
the only interactions that survive at leading order in large are those that compensate 
for every E propagator with a fermion loop. With this in mind it is easy to convince 
oneself that the exact quantum effective action of our theory has the following structure 
at leading order in large A^. 

• 1. The part of the effective action that involves only ip and ip (i.e. is independent 
of E) is quadratic in fermion fields. The quadratic form is the inverse of the 
fermion propagator. 
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2. The bare ipHip interaction in (2.41 ) is not renormalized. There are no additional 



interactions between S and the fermion fields. 



• 3. All S self interactions are obtained from the determinant in (2.44) (see the 
previous subsection). The expansion of this determinant in powers of S yields S 
self interactions of the form for all k > 2. The coefficient of the S'^ term is 
diagrammatically described by a one loop graph with k sigma lines sticking out 
of it. 

The net upshot is that, in order to compute the full quantum effective jointly for the 
ip and S fields, we need only compute the exact fermion propagator (we have already 



done this) together with the effective action (2.44) o the S fields. 

While the quantum effective action of our system is evaluated rather easily, it is not 
straightforward to determine the S propagator by inverting its quadratic action. This 
process is complicated by the fact that S is a bilocal field, and relative momentum 
(the difference between the two momenta dual to the two insertion points of S) is not 
conserved along S propagators. The quadratic term in the action for S should be 
thought of as a "matrix" in relative momentum space, and the inversion of this matrix 
requires some work. We turn to this problem in the next subsection. 



5.2 Integral equation for the four Fermi vertex 



20 



Rather than evaluate the E propagator through its effective action it is actually 
more straightforward to set up the Schwinger Dyson equations that govern the off shell 
four fermion scattering graph. This scattering process proceeds via the exchange of a 
S propagator, and so is precisely the S propagator with two fermions sewn onto each 
of the legs of this propagator. The four fermion scattering process is, in fact a more 
useful object to compute than the S propagator itself. The Schwinger Dyson equation 
that determines this object is easily derived; we find 

d'pd'qd'k^^^^^^^^ g; kW^i-k - p)ij/{p)^^\k - q)^,\q). (5.1) 



(2vr; 



9 



The axial gauge version of the integral equation described in this section was worked out in 
collaboration with J. Bhattacharya and S. Bhattacharyya. 
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J J 



1 1 



+ 



J J 





The four-fermion vertex. 



Va/3'ys{p,q', k) obeys the Schwinger-Dyson equation 

dH 



/d 



1 



z(^-f) + S(g 



t5 



(5.2) 

Substituting in the hghtcone gauge gluon propagator, and after a Fierz rearrangement, 
we obtain the integral equation 



k^ + ie 



X 



1 



i{^ + lt-^) + j:{p + k-£) 
1 



1 



7- 



7-T 



i{^ + ^-^) + E{p + k-e) 



^(^-/) + E(g-£) 
1 



aS 



aS 



(-l + ie 



(5.3) 



5.3 Effective action for fields alone 

In this brief subsection we describe an exact way of computing the scattering of an 
arbitrary number of quanta (or correlators built out of fields) working only with 
■0 fields alone. The building blocks of our construction are the exact ^ propagator and 
the exact four Fermi interaction vertex. 
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The basic rule is simple. One first draws all the Feynman diagrams that contribute 
to the process in question using the exact quantum effective action involving both E 
field and ■?/', described above. In computing a scattering process, for instance, we must 
restrict ourselves to tree graphs with the exact effective action. Graphs for correlators 
have some loops, involving the fermions that explicitly appear in the operators under 
study. 

The diagrams described in the previous paragraph will in general have some E 
internal lines, as well as some E 3 or higher point vertices. We then proceed to blow 
up every E interaction vertex into a fermion loop. Once we have done this, the only 
appearance of E in graphs is via propagators. Each propagator is sewn, at each end, 
to two fermion lines. As a consequence each E internal line carries contributes a factor 
of the fermion four point function, described in the previous subsection. The exact 
answer is obtained by enumerating all diagrams that occur in this process, writing 
down the expressions for these diagrams, and then doing the loop momenta (these 
loops arise from the blow up of E interaction vertices, as well as possible loops involving 
fields inserted in external operators). At every stage, the fermion propagator that 
participates in this process is the exact propagator. 



5.4 Exact evaluation of correlators 

The effective action of the previous section allows us, quite easily, to compute the 
correlators singlet operators, each of which is bilinear in the fermion field (see the 
previous section for an extensive discussion on operators whose correlators we evaluate). 
In the previous subsection we have already outlined a method to evaluate 3 point 
functions of the appropriate operators. The methods of this subsection generalize 
to the evaluation of arbitrary n point functions. We will make no use of the exact 
operators defined in the previous section, but use the effective action outlined in the 
previous subsection. 

In the rest of this subsection we describe in some detail how we can evaluate the 
three point functions of fermion bilinear operators using the effective action method 
(our description generalizes relatively straightforwardly to four and higher point cor- 
relators). 

In order to perform this computation we only need the fermion effective action 
expanded to sextic order in the fermions. The exact fermion propagator is determined 



by (2.2) and (2.3). The 4 Fermi interaction is determined by (5.2) or (5.3). And the six 
fermion vertex is given by a single fermion loop diagram with three attached fermion 
bilinears. The fermion propagators to be used in this loop are the exact propagators 



determined by (2.2)-(2.3), while the 4 Fermi vertex factors associated with the attached 
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fermion bilinears is the exact Fermi vertex (5.2). 



Given these ingredients, the exact three point function of the three singlet fermionic 
bihnears is given by the sum of four types of diagrams. The first is a one loop triangle 
diagram, with the exact fermion propagators, where each vertex in the triangle is one 
of the external operators. 




The exact planar three-point function. 



The second and third diagrams involve one or two four-fermion vertices. The last 
diagram involves a single sextic fermion vertex; equivalently, it consists of three four- 
fermion vertices contracted along a fermion loop using the exact propagator. It follows 



that, if we could solve both (2.3) and (5.2) we would practically be able to compute 
all "single trace" three point functions of operators that are quadratic in the fermion 
fields. 



6 Light-Cone Hamiltonian and Woo Algebra 



One of the most important features of the Euclidean field theory analysis in section 



2.1.3 of this paper is the exact representation of the path integral entirely in terms of 
the bilocal field S. The use of bilocal fields in field theory and also in the formulation 



of the large limit has been discussed by various authors, see e.g. 44 , 45 , 46 , 47 



Also, it naturally appears in discussions of non-commutative field theory 48 , 49 , 50 



See also 51 , 52 , 53 for applications to holography and higher spin theories. These 



variables, besides highlighting the inherent non-locality of the theory, are natural to 
discuss the large N limit. The large saddle point equation is the Schwinger-Dyson 



equation (2.3) for the propagator. One can then perform a systematic 1/A^ expansion 



around the saddle point solution. 
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In the following, we will briefly explore the Hamiltonian formulation of the theory 
in light-cone time and in = gauge. Such a treatment of the path-integral allows 
one to express the theory in terms of 'equal time' bilocal variables M{x, y) satisfying 
a Woo algebra. These operators are also constrained by a quadratic relation = 
M, that has appeared before in discussions of non-commutative fleld theory. The 
Hamiltonian formulation throws a different light on the structure of this theory, and 



the A dependence in the exact fermion propagator (2.23) arises in a simple way when 
evaluated at equal light-cone time. 

We will flrst construct a light-cone Hamiltonian in terms a bilocal fleld M{x, y) 
which is deflned on a flxed x"*" slice, and then use this Hamiltonian to derive an action 
for M(x, y) using the method of coadjoint orbits. Here will serve as a time coordinate 
and the Hamiltonian that we will construct is conjugate to this variable. Our analysis 



closely follows 47 , applied to the large N limit of QCD2. 

We pause to remark that light-cone quantization of any theory has a number of 
important issues, and (a priori) is not necessarily equivalent to usual time-like quanti- 
zation. See, e.g., [54]. While we leave careful consideration of these issues for future 
investigation, we will observe that, applied with appropriate qualiflcation, the light- 
cone quantization appears to give some simplifying insights into the dynamics of the 
theory. 

We will be working entirely in Minkowski space-time for this section. We use the 
following explicit representation for gamma-matrices: 

The fermion is a two-component Dirac spinor which we write as 

These conventions are similar to the Euclidean conventions with a few minor differ- 
ences. The lightcone coordinates are a sum of x° and x^ and the non-light cone direction 
is (instead of ). In Minkowski space, it is natural to deflne = {±x^ + x^)/\/2, 
which happens to equal x^^^^ , as used elsewhere in the paper. 

The action in written in Minkowski space is 
S = I ii( A ^dA + AAAAj-fd^x [i^j'-j^'D^^ + im^^) (6.3) 
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We obtain the Hamiltonian for the theory in the gauge — 0. As ahcady men- 
tioned is the 'time' coordinate, and and arc the space coordinates. The Gauss 
constraint arising from the equation of motion of is: 

- ^ - d2A'L - r'"'A''_Al) + V2iljlrtlj_ = (6.4) 

The momenta conjugate to -0+ and vanish - this reduction in the number of 
degrees of freedom is a consequence of the choice of hght-cone direction x'^ as the time 
variable - the corresponding constraints arising from their equations of motion are: 

iV2d^-^+ + id2ip- - imilj_ + v^A^T''^+ + A^T^'iP^ = (6.5) 
-i\/2a_^i - id2iljl + imipl + i)\^/2A''_T'' + V^lyl^T" = (6.6) 

The momentum conjugate to A" is 

= = -4^-^5 (6.7) 

The momentum conjugate to -0- is 



7T^i_ = ^ = iV2iljl (6.8) 



The Hamiltonian then takes the following form 

n = i7pld2ip+ + imTplijj^ + 'iplA'^T''ij;+. 

Using the residual gauge symmetry to set A_ = 0, we can explicitly solve the 
constraints to express A2, ip+ and ip^. in terms of ip- and its conjugate momentum ■j/^l. 

The resulting (spatially non-local) Hamiltonian, expressed in momentum space is: 



k J An"^ Att"^ Att'^ \ P-Q- 
271 f (fp (Pq (Pp' f ip2 — iq2 — m 



k J Att^ Att"^ Att^ \ g_ (p_ — q. 
-Att^ f <fp <fq <fp' <Pq' <fq" 



^L(-p) (^1(? -p')V'-(p - ct)) Mp') 



r f 

^ ~W J An^ 47r2 An^ An^ An^ 

(6.9) 
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Note that, because we are using light-cone time, and we must choose some prescrip- 
tion to define the inverse of the hght-cone momenta when one solves the constraints 
(6.5) and (6.6) in momentum space. Choosing the principal value prescription of reg- 
ulating the infra-red divergence effectively eliminates the zero mode of The only 
dynamical degree of freedom in the theory is then As pointed out in 54 in the 
limit of large A^, using the principal value prescription in the definition of infra-red 
divergent quantities like in the related theory of two-dimensional QCD2, one 

can obtain the correct 't Hooft equation that determines the meson spectrum. 



6.1 Bilocal variables satisfying a Woo algebra 

Let us define the following bilocal operator 

M{x,y) = ^^]_{^)^iiy) (6.10) 

These operators, which would be the gauge invariant Wilson lines along the x~ direction 
in the A_ gauge if the points x and y were not separated in the 2-direction, constitute 
the quantum phase space of the light-cone quantized formulation of the theory. 

As we will be working in momentum space, it will be more convenient to use 

M(g,g') = J^^U-^)^~(^') (6-11) 

It is easy to calculate the commutation relations using the basic fermion anti- 
commutation relations: 

{^l,(fc),^L(fc)} = -^SpT^rS'ik + k'). (6.12) 

These imply that the M's satisfy the 'quantum' VFoo algebra 

[Mip,p'),Miq,q')] = (M(p, g')5'(p' - q) - M{q,p')5\q' - p)) (6.13) 

Notationally, we may often represent M as an operator acting on a single particle 
Hilbert space, M, with 

{x\M\y) = M{x,y) (6.14) 

though M(x, y) is itself an operator acting on the Hilbert space of the quantum field the- 
ory, i.e., multi-particle states. Such constructions are familiar from non-commutative 
field theory. 
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In terms of variables, the Hamiltonian, using the principal value prescription, 
can be easily expressed as a 'cubic polynomial' in M: 

-NX^ I dbd\d'q''^^M{p,q')M{q' -q,p-q) 
Att J p-q- 

-NX— f d''pd\d\' '^^ ~ ~ "^ Mip, q')M{q' -q,p-q) 
47r J q-{p- - q-) 

- j d^pd^d^d^'dY^Y^Mip' + q",p)M{q + q",p')M{p -q',q- q') 

(6.15) 

The inverse of p_ in the above manipulations is assumed to be defined using the 
principal value prescription. 

6.1.1 Constraint on the bilocal fields 

We now restrict the space of allowed states to be singlets under global gauge transfor- 
mations i.e. the generators of global SU {N) gauge transformations constrain the wave 
functions as follows: 

J dx-dx^ (^^IV- - 1^) = (6.16) 

We also note that the baryon number operator 

^/2 



B - 



^ j d'p^i\-p)r-{p) (6.17) 



N{2i 

commutes with the Hamiltonian and hence it is conserved and takes a definite value on 
physical states. This fact leads to a very important operator constraint in the phase 
space. Prom the definition of M, it follows that 

M^{p,p') = j dqM{p,q)M{q,p') = {1 + B)M{p,p') (6.18) 
B ^ [ d^pMip,p) (6.19) 



This completes the light-cone Hamiltonian formulation. The (light-cone quantized) 
phase space operators arc the Af 's. Their commutation relations satisfy the Woo Lie 
algebra. Since the generators of this Lie algebra are represented in terms of the un- 
derlying fermions, they are not all independent and satisfy a constraint which is a 
consequence of gauge invariance and a conserved baryon number. It is clear that we 
are dealing with a W^o Hamiltonian system, whose representation is fixed by the con- 
straint — {1 + B)M. We will restrict ourselves to the B — sector of this theory. 
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6.2 The path integral and the large- A/^ classical solution 



From the light-cone Hamiltonian one can construct an action (and the path integral) 
for the theory using standard techniques of geometric quantization viz. the method of 
coadjoint orbits of the (Woo) Lie group. Equivalently the path integral can be derived 
using a basis of coherent states of the M operators between adjacent time steps. The 



action one obtains in terms of the 'c-number' bilocal fields is (see 47,49 and references 
therein), 

^classical = / dsdxHiM[dsM, d+M] - j dx+ H (6.20) 

The first term of the action is the Kirilov-Berry phase. It is written using an 
auxiliary parameter s besides the time t. If t is periodic then s parameterizes a unit disc. 
We need to impose the boundary condition M(s = 0,t) = and M(s = l,t) = M{t). 
The classical equations of motion are obtained by varying the action consistent with the 
Woo group law and using the above mentioned boundary conditions. These equations 
are by construction the Hamilton equations of motion: —id^M = [H, M]. 

The path integral measure over M is constrained by the covariant condition = 
M, which ultimately arises from evaluating the operator constraint in the space of 
coherent states (see |49]). In particular, to obtain the solution (M) that corresponds 
to the ground state expectation value of M, we solve the time-independent equation 
[H, M] = 0. 

From the commutation relations of the Woo variables, one can see that any trans- 
lationally invariant (M) 

{M{p,q)) = 6'{p-p')f{p) (6.21) 

is a solution provided the constraint (M) = (M) is satisfied. This implies that f^p^ = 
f{p), which being independent of A, is most naturally solved by f{k) = 6{— ^^^"^ ) ~= 
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It is important to note that this solution is the exact solution in the large- limit. 
It allows one to read off the equal-time correlation function for tp-^s, valid to all orders 
in A, 

(^;_(g')^Ug)) =ivi=(2vr)2(M(-g',g)) (6.22) 



which implies: 



j-y.^^-y.^, V2 
1 



mq)^UQ')) = N—e{q^){2'nf5\q + q') (6.23) 

This is, in fact, one of the components of the equal- "time" fermionic propagator. Using 
the constraints, we can use this solution to obtain other components of the exact 



'^We thank J. David for discussion on this point. 
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propagator. The fermionic constraints in this large- A?" solution read: 

i/j+iq) = {iq2 - (m + \h{-q))) (6-24) 

where h{-q) = ^ J <fpj^M{q-p). 

Let us discuss the regularization of (m + /*(— g)). First, let us note that if q_A{q) = 
J{q), then A{q) = + h{q2)5{q-))J{q) is a solution for any /i(g2), because 
is effectively equal to zero, (assuming J{q) 7^ at g_ = 0). This is ultimately an 
infrared ambiguity in the inverse of 9_. 

Second, note that we will also regulate the vacuum solution M{q—p) = lima.2^0 (^{q— 
pyk2x^ so that / £pM{p,p) = B ^ 0. 

Taking both the infrared ambiguity and the regularization of M into account, we 
have 

h(q-) dp2dp_ (^^ + h{p2)5{p_)^ e'^^--^->'e{q^ 

= 7^ / dp2dp_h{p2)5{p_)e'^'i^-P^>^e{q_-pJ) 
27r J 

= / dp2h{p2)e-'^''"e{q_) 

Note that although h{p2) is not specified, it must be dimensionless. Let us take h{p2) — 
hoe-^^^^''^' cos{p2X^). Then 

I^{q_) = ^e^^^'^y dp2e-^P'''"^" cos{p2X^)e-'P-''''e{q_) 

j rfpse-^^^^')' cos2(p2^^)^(?-) 



27r 

- {l/x^ + iq2)e{q_) 

where, in the last line, we took the limit x'^ — > 0, and adjusted the numerical constant 
/lo- 
in the limit — > we get a divergent term independent of momentum, which we 
cancel with the mass counterterm, and a finite term remains, iq2] hence (m-|-/*(— gr)) = 
iq2. Regardless of the choice of /i, by dimensional analysis the finite part of the integral 
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lim^2_i.o e^'i'^^^ J dp2h{p2)e~''^^^^ would have to be proportional to q2, so our answer would 
be the unchanged up to a redefinition of A by a numerical factor, assuming h{p2) is 
real and reasonably well-behaved. 



Using (6.24) and (6.25), the other components of the exact propagator are: 

-^2(1- A) 



{^P+iq)^|JUq')) 
i^PM^iiq')) 



2g_ 

-g2(i + A) 
g2^(l-A^) 

2V2q^ 



Ne{q^){27ry6\q + q') 
N9iq^)i2nfS'iq + q') 
Ne{q.){27,f5\q + q') 



(6.26) 
(6.27) 
(6.28) 



Let us compare with the equal time propagator derived earlier in section [2| obtained 

dkA 



from (2.23) 



lim 

x+-s>0+ / 271 



hm I ^ ' ^fj^^^^^ + ''^^''\ \2nmk + k') 

r+^o+ 271 k^ + ie ! y J K ) 



(6.29) 



Comparison with ((6.23)-(6.28)), immediately implies that fo+go = 0, and /o ~ A. 
We see that, because of the A independence of the vacuum, all A dependence arises 



from the fermionic constraints (6.5), (6.6). 



Note that in this discussion we did not use dimensional regularization which was 
used in the Euclidean calculation and hence this agreement is significant. 



6.3 Large N perturbation theory and a formal large solution 

Now that we have the classical solution at large N, we can construct a Woo covariant 
perturbation theory by parameterizing 



(6.30) 



in terms of the fiuctuations W. 



The W^s satisfy an algebra inherited from the commutation relations of the M's, 
and break up into subalgebras when expanded around Mdass. exactly as in f47|, where 
a similar treatment was carried out for QCD2 that gave rise to 't Hooft's well known 
equation for the spectrum of mesons. 
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The fluctuation equation obtained from substituting (6.30) into —id^M = [H,M] 



is a second order integral equation for W. It takes the following schematic form: 
d+W{p,p';x+) =- ( pIIHI. _ dial.) W{p,p';x+) 



P- 



+ X J d^qi...W) + X^ J d\d\' (...W) 



(6.31) 



where W{p,p'; x~^) = W{p,p')9{p^) 

Any single-trace operator (Jq, Ji, J2 etc.) can be expressed in terms of M's using 



solution to the constraints, (6.4), (6.24) and (6.25), and Hamilton's equation for d+ip 



hence, any two-point function of J's can be expressed in terms of two-point functions of 
M's. Because {M{p,p + k)) = for nonzero k, the Green's function for the fluctuation 



equation (6.30) is needed to calculate the leading nonzero contribution to a generic 



two-point function of M's in the 1/N expansion. Thus, solving equation (6.31) for 



PFoo-covariant fluctuations around the classical solution would, in principle, enable one 
to calculate exact two-point functions of all single trace operators in the large limit 
(subject to, of course, subtleties associated with the light-cone Hamiltonian formalism.) 



7 Comments on the Holographic Dual 

7.1 Vector models are dual to higher spin gauge theories 

The non-renormalization of single trace currents in U{N) Chern-Simons-fermion theory 
at leading order in 1/A^ indicates that its holographic dual has a free spectrum of purely 
higher spin gauge fields, of spin s = 0, 1,2,3, ■ ■ ■ , and their bound states (which are 
dual to multi-trace operators). The divergences of the single trace currents are mixed 
with double and possibly triple trace operators, and the mixing is suppressed by . 
By comparing two and three point functions of the currents, we see that N~2 should 
be identified with the bulk coupling constant. 

Classically, the bulk nonlinear equations of motion must therefore preserve the 
higher spin gauge symmetries as well, for otherwise they would introduce extra lon- 
gitudinal degrees of freedom for the higher spin fields even in the zero coupling limit. 
This means that the dual bulk theory is a pure higher spin gauge theory. On the other 
hand, the higher spin symmetry in the Chern-Simons-fermion theory is broken by 1/A^ 
corrections, and we expect the HS gauge symmetry to be broken in the bulk by bound- 



ary conditions through loop effects and the HS gauge fields become massive through 



22 



This has been shown to occur, for instance, in Vasihev's A-type theory in AdS4 with the A = 2 



boundary condition, dual to the critical 0{N) model 33 . 
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the mixing of its longitudinal mode with bound states. 

Our argument of non-renormalization of singlet bilinear currents extends to general 
vector models coupled to Chern-Simons gauge fields, such as the various supersym- 

metric extensions of the Chcrn-Simons-fcrmion theory discussed so far. We therefore 
anticipate all such theories to be dual to some higher spin gauge theory in the bulk. 
The bulk theory is generally not parity invariant. 

7.2 Parity violating Vasiliev theory 

More precisely, we conjecture that the three-dimensional critical U (N) Chern-Simons 
theory at level k coupled to a single fundamental fermion in the 't Hooft limit ^ 1, 
A = N/k finite, is dual to a parity- violating higher spin gauge theory in AdS4 that 
contains one gauge field for each non-negative integer spin s = 0,l,2,---. A class of 
such theories has indeed been constructed by Vasiliev, generalizing the parity preserving 
A-type and B-type theories (which have been conjectured to be dual to free/critical 
U{N) scalar and fermion vector models in the singlet sector). We will refer to them as 
parity- violating Vasiliev theories in AdS^. The bulk interactions of the parity violating 
Vasiliev theory are governed by a single function 

/(X) = l + Xexp(i^(X)) (7.1) 

where 0{X) is a real, even function in X, but otherwise unconstrained (in the classical 
theory). Roughly, when expanded in powers of X, the order X'^ term in f{X) governs 
the {n + 2)-th order coupling of the bulk theory. 

Vasiliev's theory has the feature that the bulk scalar field has mass — — 2/£^^5, 
and the dual scalar operator (with one choice of boundary condition) has classical 
dimension 2. This does not follow from the linearized higher spin gauge symmetry, 
and is rcqTiircd in order to construct higher spin gauge invariant interactions. This 
feature is in precise agreement with the non-renormalization of the scalar operator ipijj 
at the planar level in the Chern-Simons-fermion vector model, which we have argued 
using the anomalous current conservation relation and verified explicitly at two-loop. 

Now let us discuss the parity transformation in Vasiliev theory. The function f{X) 
transforms to its complex conjugate f{X) under parity. One may either assign even or 
odd parity to the bulk scalar field, which corresponds to assigning even or odd parity to 
the variable X. In the former parity invariant theory is given by f{X) = f{X), 

or f{X) — 1 + X; this is the A-type theory. In the latter case, where X is parity odd, 
a parity invariant theory is given by f{X) — 1 + iX; this is the B-type theory. The 
B-type theory has been conjectured to the dual to the free limit of the U (N) Chern- 
Simons-fermion theory, where the Chern-Simons gauge fields decouple and their only 
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effect is to restrict the operator spectrum to gauge singlets; parity is restored in this 
free hmit. 

It is thus natural to suspect that the Chern-Simons-fermion theory at large but 
finite 't Hooft coupling A = N/k, is dual to the parity violating Vasiliev theory governed 
by a function f{X) = 1 + X exp(i6'(A, X)), where 

^(A,X) = ^ + /3(A,X), (7.2) 

is real, even in X, and depends on A in some way. For instance, one should have 
/3(0, X) = 0, so that in the A = case the conjecture reduces to the duality between B 
type theory and free fermions. While the CS-fermion theory violates parity symmetry, 
parity is restored if we flip the Chern-Simons level k (i.e. A — >■ —A) together with the 
standard parity action on the fields, under which the scalar operator -xfjilj is parity odd. 
So if the Chern-Simons-fermion theory is dual to the parity violating Vasiliev theory, 
we must also have /3(— A,X) = — /3(A,X). 

While there is a priori an infinite family of parity-violating Vasiliev theories in 
AdSii governed by the function /(X), we may speculate that there is a particular 
choice of /(X) determined by the 't Hooft coupling A in the above form, that describes 
the bulk theory dual to critical U{N) CS-fermion theory at large N. It is then an 
interesting question to determine ^(A,X) as a function of A. In principle, this can 
be done perturbatively by comparing the n-point functions of the HS currents in CS- 
fermion theory and the bulk Vasiliev theory. We may expand the function /3(A,X) 
as YlT=o f^kW^'^'' ■ Then may in principle be determined by comparing with the 
{2k + 3)-point function of the higher spin currents in the CFT at leading order in 
1/N. In particular, the leading phase (3o{X) can be determined by comparing with the 
three-point functions of currents. 

The simplest check of such a conjecture would be the comparison of bulk tree level 
three point functions in parity-violating Vasiliev theory, and the planar three point 
function {J J J) in the CS-fermion theory. The former computation, performed along 
the lines of [32], will be presented elsewhere. The result obtained using the gauge 
function method of ^2] ('W = gauge") indeed agrees with the parity-even part of 
{J J J) in CS-fermion theory up to two-loop: in particular, the property that Cb in 



(4.34) is independent of the three spins is a prediction from the bulk theory! However, 



we have failed to produce the parity-odd structure of {J J J) from the bulk Vasiliev 



theory. This could be either due to the singular nature of the 1^ = gauge in 32 and 



that the parity-odd contribution is simply missed by the regularization of 32 , or it 
could be that the correct dual bulk theory to the CS-fermion CFT is not described by 
Vasiliev's system, but rather some other higher spin gauge theory with the same field 
content in AdS^. We shall leave the resolution of this puzzle to future work. 
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Let us comment on higher spin symmetry breaking as seen from the correlation 
functions. While at tree level in the bulk, we expect the higher spin symmetry to 
be preserved, there can be subtle boundary contributions in computing correlators. 
Examples of similar nature in the bulk dual of critical 0{N) model have been discussed 
in [33]. Recall the anomalous current conservation relation in the boundary theory, of 
the form 

d^jjflr^-^ E i/"^] ® [/"^] + ^ E [/"^] ® [/"^] ® [/"^] (7.3) 

On the right hand side, we have listed the conformal families that can appear in the 
divergence of the current j The dependence on 't Hooft coupling A is omitted in this 
expression. As discussed in section [3} though these multi-trace operators appear to be 
suppressed by powers of they can contribute to the corresponding divergence of 
correlation functions at the leading nontrivial order in 

For the three-point functions (j^^'^^j^^^^j^^^^), the naive current conservation equation 
at leading order in namely 0{N~^) when the two point functions are normalized 
to 0{N^), 

d^Hjl7^.j'-''h^"'^) = (7.4) 

is necessarily satisfied if gi < ^2 + -SsF^ 

On the other hand, if si > S2 + S3, then there is the possibility of a nontrivial 
0{N~2) contribution, 

^'^^(jt^.b-^^^^J^^^^) ~ -^{J^''^J^''^){J^''^J^''^^). (7.5) 
V i\ 

So even at leading order in the naive current conservation relation is violated by 
the three point function. Naively, the boundary-to-bulk propagator of the higher spin 
fields in AdS^ is divergence free with respect to the position and polarization of the 
source current, and it may seem that such a non-conserved three point function cannot 
result from a tree level bulk computation. As pointed out in [33], however, this is not 
necessarily the case, due to subtle boundary contributions to the bulk integral. 

^■^In the explicit computations we have performed for the classical divergence of currents for spins up 
to 5 we in fact find that this the current is effectively conserved in correlators provided si < S2 + S3 + 2. 
It is possible that this relation is more generally true and can be proved from an analysis of structures 
allowed by conformal invariance in 3 point functions. If this is the case the three point function will 
be divergence free with respect to all three currents if the three spins si, 82,33 obey this less strict 
triangular inequality. 
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7.3 A higher spin gauge theory as a string theory 



A long standing question is whether the pure higher spin gauge theory in AdS can be 
embedded in string theory. We have seen that quite generally vector models coupled 
to Chern-Simons fields are dual to higher spin gauge theories. On the other hand, a 
large class of supersymmetric quiver Chern-Simons-matter theories are dual to string 
theories in AdS^. A particularly interesting class of examples is that of the ABJ 



theory 12 : the M = Q U{N)k x U{M)_k Chern-Simons theory coupled to a pair of 
bifundamental hypermultiplets. It has been conjectured to be dual to type IIA string 
theory on AdS 4, x CP'^, with a flat 5-fleld obeying 

1 f N-M 1 ^ ^ 

M\ < k. The need for the shift by | in the equation above was explained 



when 


N 


in 55 
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While the dual string theory is weakly coupled in the 't Hooft limit, the duality is 
believed to hold exactly for all A^, M and k. In the special case M = 1 (or generally 
small M), ABJ theory reduces to a supersymmetric Chern-Simons vector model, in the 
following sense. The most straightforward A/" = 3 supersymmetric generalization of the 
vector model we considered so far is the A/" = 3 f/ (A^) or SU {N) Chern-Simons theory 
coupled to a fundamental hypermultiplet. Most features of large A^ vector models 
discussed in this paper still hold when a number of fundamental flavors are introduced 
(the number of flavors does not scale with A^ in the large A^ limit here). Now consider 
Af = 3 U{N)k Chern-Simons theory with 2M hypermultiplets. By essentially the same 
argument as in section 2, at large A^ this vector model is dual to a higher spin gauge 
theory in AdS4^ with U (2M) Chan-Paton factors, To obtain the Af = 6 theory, one 
needs to further gauge the U (M) flavor symmetry with A/" = 3 Chern-Simons coupling 
at level —k (as opposed to level inflnity). This however corresponds to simply modifying 
the AdS boundary condition on the bulk U{M) spin-1 gauge flelds, from the purely 
magnetic boundary condition to an electric-magnetic mixed boundary condition [62^. 

Let us comment briefly on the distinction between U{N) and SU{N) gauge groups 
here. Take the ABJ theory with M = 1 as an example. The gauge group is U{N)k x 
U{l)-k. Denote by a and a the diagonal U{1) of the U{N)k and the f/(l)_fe gauge 
flelds, respectively. Note that a has Chern-Simons level Nk as a U{1) gauge fleld. The 



^''We thank O. Aharony for bringing this to our attention. 

^^Note that Vasihev's system can be generahzed to a theory of nonabehan higher spin fields by 
introducing U{Nf) Chan-Paton factors, i.e. replacing the * algebra by its tensor product with the 
algebra of Nf x Nf matrices. 
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matter fields couple to 6 = a — a only. The Chern-Simons action for a, a takes the form 

k f k f 

{Na Ada- a Ada) = — / {Na A da - {a - b) A d{a - h)). (7.7) 



47r y 47r „ 

Since only h couple to the matter fields, we can integrate out a exactly, which sets 
a = —jfzib. The resulting Chern-Simons action is that of the U{1) gauge field b at 
level —j^k. So the same theory can be expressed as SU {N)^ x U {l)-k' Chern-Simons 
coupled to 2 bifundamental hypermultiplets, where k' = j^k. 

We are now naturally led to the conjecture that the U{N) x U (M) AB J theory in 
the vector model limit, i.e. large N and fixed finite M, is dual to an A/" = 6 version of 
the parity violating Vasiliev theory with U (M) Chan-Paton factors. Furthermore, we 
conjecture that the type IIA string field theory on AdS^ x CP^, with a fiat -B-field as 



in (7.6), in limit 3> M, is a higher spin gauge theory. Note that in the vector model 
limit A^ 3> M, the 't Hooft coupling N/k is restricted to be less than 1, and the bulk 
IIA string theory is always in the highly stringy regime. 

In the ABJ theory, the bulk single string states are dual to single trace operators, 
of the form 

Ti{ABAB ■■■) (7.8) 

where A and B denote fields in the (□, □) and (□, □) representations respectively. In 
the M = 1 case, the single trace trace breaks into the product of a number of bilinears 
which are singlets under the SU (N) gauge group. So we learn that the general single 
string state should correspond a multi-particle state of higher spin fields. It suggests 
that the higher spin gauge fields are in fact more fundamental than strings in this case, 
as the latter are composite objects made out of the former! Note that the coupling 
constant g of the higher spin gauge theory is mapped to ~ 1/ \/N of the vector model. 
Now Xbuik = g'^M ~ M/N is the bulk 't Hooft coupling. We see that the vector model 
limit A^ ^ M is precisely where the nonabelian higher spin gauge theory is weakly 
coupled. As M/N becomes of order 1, the higher spin fields interact strongly, and we 
expect their bound states to form the string excitations. If the dual higher spin gauge 
theory is indeed described by Vasihev's system, it would be very interesting to find 
the explicit map from Vasiliev's master fields to the (closed) string fields, which could 
allow for a precise formulation of the closed string field theory in this limit. 



8 Discussion 

This paper has been devoted to the study of U{N) Chern-Simons theory coupled to 
a single massless fundamental fermion. We set out to to solve this theory in the 't 
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Hooft large N limit and have been partly successful in this task. In particular, we were 
able to compute the partition function on of the theory at finite temperature at 
all values of the 't Hooft coupling. In order to perform this computation we employed 
a lightcone gauge, and used a dimensional reduction regularization scheme. It would 
be useful to perform further checks of the consistency of our gauge and regularization 
scheme. 

In the course of our analysis of this theory we have encountered the fact that the 
higher spin currents obey anomalous conservation equations. In the classical theory 
we have computed the explicit form of these conservation equation at low values of 



the spin (see e.g. (3.19)). These nonlinear anomalous conservation equations contain 
a large amount of information; for example they encode the anomalous dimensions of 
the spin s currents in a ^ expansion (see subsection 3.4). It would be interesting to 



determine the explicit form of these anomalous conservation equations, if possible, as 
a function of A. 

We have also argued that the spectrum of local gauge invariant operators in our 
theory, at dimensions of order unity, is not renormalized as a function of the 't Hooft 
coupling at leading order in large A^. A class of such operators (those that may be 
written entirely in terms of x~ derivatives) are dressed due to interactions only by 
rainbow graphs in the large N limit. The "quantum dressed" operator obtained in 
this way is given by the solution to an integral equation both at zero and at finite 
temperature. Exact expressions for all two and and three point functions are obtained 
by sewing these quantum dressed operators together using the exact fermion propagator 
computed in this paper. Such a solution would allow us to find exact expressions, for 
example, of fermion number current and stress tensor two point functions at finite 
temperature as a function of k and u. As these quantities are of obvious physical 
interest, the corresponding integral equations deserve serious attention. 

In this paper we have analyzed in detail the spectrum of local operators. But clearly 
the theory also admits interesting non-local gauge invariant observables, namely the 
Wilson loop operators specified by a choice of a closed contour in space-time and a 
representation of the gauge group. In pure Chern-Simons theory, these are topological 
observables which compute certain knot invariants [3]. When Chern-Simons theory 
is coupled to matter, the Wilson loops are no longer topological. In a theory with 
fundamental matter, however, at leading order at large A^ their expectation values are 
the same as in pure Chern-Simons theory, since diagrams with fundamental matter 
loops are suppressed. It may be interesting to explicitly compute the leading 1/A^ 
corrections for some Wilson loop in our theory. In standard examples of AdS/CFT, 
Wilson loop operators are dual to macroscopic open strings ending on the closed contour 
at the boundary. In this paper we have proposed that our theory should be dual to some 
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higher spin gauge theory in AdS4. It would be interesting to understand if Vasihev's 
theory, or a deformation thereof, may possibly contain extended objects which could 
be dual to the Wilson loops of the gauge theory. 

It has recently been realized that the partition functions of supersymmetric Chern- 
Simons theories on are often exactly computable via the techniques of supersym- 
metric localization 56 59 . This quantity appears to play the role of a c function under 
the renormalization group flow and so is of clear physical interest. It seems possible 
that the techniques employed in this paper will allow the exact computation of the 
partition function of our non-supersymmetric theory in the large limit, as a function 
of A. 



As we have explained in subsection 3.1.3 above, in the limit A — 0, the partition 
function of our theory on S"^ x undergoes an interesting phase transition at a temper- 



ature of order 0{\/N) 39 . The low temperature phase may be thought of as a gas of 
single traces. At temperatures that are held fixed when A^ is taken to infinity, the single 
trace anomalous dimension non renormalization theorems of this paper demonstrate 
that the free energy of this phase is not renormalized as a function of A. At tempera- 
tures at of order \/ N the free energy of of the low temperature phase presumably does 
get renormalized as a function of A. And the free energy of the high temperature phase 
of the partition function (which be thought of as a gas of the fundamental fermions 
rather than their mesonic bound states) certainly does get renormalized a function of 
A. In this paper in particular we have computed the coefficient of the growth of the 
free energy in the very high temperature limit, and have shown that this coefficient, 
which is nonzero at generic A, vanishes in the limit A — > 1. It would be very interesting 
to understand how the competition between the two phases changes as a function of A, 
especially in the limit A — )■ 1. Does the vanishing of the free energy at extremely high 
temperatures, for instance, suggest that the phase transition temperature between the 
low and high temperature phases diverges as A — )• 1? 

It would also be interesting to investigate transport properties in the high tempera- 
ture phase. The following argument suggests that the mean free path of any particular 
fermionic quantum, is of order A^. Consider, for instance, a fermion of colour i scat- 
tering off an anti fermion of colour i to yield a fermion - anti-fermion pair of colour 
j. This scattering amplitude is of order as a consequence of which the scattering 
probability is of order The number of potential scattering partners is of order 
unity in the thermal bath, while there are A^ possibilities for the colour j of the final 
fermion pair. It follows that the probability that any given fermionic quantum un- 
dergoes a scattering process is of order giving rise to a mean free path of order 
A^. This result should apply at every value of A. All this very similar to the known 
behavior of transport properties in the critical 0(N) model (see e.g. 22 ) and very 
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different from theories with matrix degrees of freedom in the t' Hooft hmit where the 
mean free time between colhsions of a single parton is independent of f60] . The fact 
that the fermion gas effectively behaves as a free gas as far as transport properties are 
concerned, suggests that fluctuations about the high temperature phase do not decay 
over a time of order unity. This discussion suggests that large vector models in the 
't Hooft limit thermalize much less efficiently than their matrix model counterparts. 
For this reason the dual description of such a high temperature phase has properties 
that are very different from a conventional classical black hole. |^ 

It should also be possible to generalize our discussion of the finite temperature 
behaviour of our theory to include a chemical potential for fermion number. The finite 
A behaviour of such a system describes an interacting Fermi sea in three dimensions, 
and so may be of interest for various condensed matter problems. 

It would be interesting to generalize our computation of the finite temperature free 



energy to Chern-Simons theory interacting with fundamental bosons 37 , and also to 
the supersymmmetric theory with different numbers of supersymmetries. 

An outstanding question about the theory studied in this paper as well as the gen- 
eralizations described above is "what is its bulk dual description?" We have argued 
above that this dual description is a higher spin theory. However we do not yet have a 
precise conjecture for the nature of this dual; this is clearly an issue of immediate inter- 
est, especially since the maximally supersymmetric version of our theory is embedded 



into string theory as a limit of the ABJ construction (see subsection 7.3). 

In this paper we have studied only the 't Hooft large A^ limit of our theory. It is 
however clear from the example of ABJM theory, that large A^ Chern Simons theories 
admit a wider class of large A^ limits (e.g. A^ — )■ oo, A — )■ oo, y fixed in the ABJM 
context). It would be interesting to investigate whether the theory studied in this paper 
admits an analogous double scaling limit, in which A is scaled to unity in a manner 
coordinated with the large A^ limitP^ The dimensions of the higher spin currents could 
be renormalized as a function of coupling in such a limit. The fact that we have been 
able to obtain exact results as a function of A perhaps makes such an investigation 
feasible. 

To end this discussion we recall we could very easily argue that the spectrum of 



^^When the ranks, iVi and N2 of the ABJ theory are both comparable, the theory undergoes a 
deconfinement transition, dual to a Hawking Page transition, at a temperature of order unity. The 
high temperature phase in this transition is a black hole in Einstein gravity. It follows from the 
discussion above, on the other hand, that when A^i is lowered to a number of order unity with N2 
still large then the deconfinement phase transition temperature is much larger (it scales like \/N2). 
Moreover the dual black hole appears to change its nature qualitatively; it appears to stop absorbing. 

^^We thank O. Aharony for discussions on this point. 
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'single trace' operators in our theory is not renormalized as a function of A, in the 't 
Hooft large limit, by combining conformal representation theory with the sparseness 
of the single trace spectrum in vector models. The same non-renormalization argument 
does not apply in theories with adjoint or bifundamental matter fields as the single trace 
spectrum of these theories is not sparse]^ On the other hand, as we have explained 
in the introduction, non supersymmetric effective large Chern-Simons fixed lines 
are easily constructed with matter fields in adjoint or bifundamental representations. 
The scaling dimension of single trace operators in such theories is protected neither 
by conformal nor supersymmetric representation theory. In these theories all single 
trace operators baring the stress tensor and currents for global symmetries can, and 
presumably do get renormalized. It is at least conceivable that some theory with adjoint 
or bifundamental fermions admits a strong coupling limit in which all but a finite 
number of single trace operators are infinitely renormalized, and the dual description 
of the theory is Einstein gravity in AdS4, coupled to a minimal number of additional 
fields. It would clearly be very interesting to identify any theory with this property. 

As is clear from the discussion above, the computations presented in this paper 
have merely scratched the surface of a large and potentially very interesting area of 
investigation. We hope to report on some of the topics discussed above in the future. 



Acknowledgements 

We are grateful to Ofer Aharony and Rajesh Gopakumar for crucial discussions that 
led to the formulation of some of the questions addressed in this paper. We would like 
to acknowledge Sayantani Bhattacharyya, Jyotirmoy Bhattacharya and V. Umesh for 
collaboration at various stages of this project. We thank O. Aharony and J. Bhat- 
tacharya for comments on the manuscript. We would also like to acknowledge useful 
discussions with C.-M. Chang, K. Damle, J. David, A. Dhar, S. Dutta, D. Hofman, 
N. lizuka, D. Jafferis, S. Kim, S. S. Lee, G. Mandal, K. Papadodimas, M. Rangamani, 
S. Sachdev, A. Sen, C. Sonnenshein and E. Witten. S.G. is supported by Perime- 
ter Institute for Theoretical Physics. Research at Perimeter Institute is supported by 
the Government of Canada through Industry Canada and by the Province of Ontario 
through the Ministry of Research & Innovation. The work of S.M. was supported by a 
Swarnajayanti fellowship. The work of S.M., S.P., S.T. and S.W. was supported by the 
DAE, government of India. S.G. would like to thank Princeton U., the Simons Center 
for Geometry and Physics, Mc Gill U., the Weizmann Institute, the organizers of the 

^^In supersymmetric theories with such matter content, supersymmetric indices and representation 
theory give rise to a new source of non-renormahzation theorems for single trace operators even in 



theories with bifundamental and adjoint matter, see e.g. the paper 61 . 



77 



PhyMSI conference in Cargese and of the Simons Summer Workshop in Mathematics 
and Physics 2011 for hospitality during completion of this work. S.M. would like to 
thank the Weizmann Institute, Cambridge University, KIAS, ICTP, the organizers of 
GR19, PASCOS 2010, and the Indian Strings Meeting (ISM) 2011 for for their hos- 
pitality while this work was in progress. S.P. would like to thank the organizers of 
ISM 2011, and the Centre for High Energy Physics, Indian Institute of Science, Ban- 
galore for their hospitality. X.Y. would hke to thank the organizers of GR19 and of 
ISM 2011, Tata Institute of Fundamental Research, Berkeley Center for Theoretical 
Physics, Simons Center for Geometry and Physics, Fields Institute, and Aspen Center 
for Physics for their hospitality. The work of X.Y. was supported by the Fundamental 
Laws Initiative Fund at Harvard University and in part by NSF Award PHY-0847457. 
S.M., S.P., S.T. and S.W. would also like to acknowledge our debt to the people of 
India for their generous and steady support to research in the basic sciences. 



Appendices 

A Conventions for Propagators and Gauge Condi- 
tions 

The Euclidean action for our theory is 

5 = ^ y Tr (^AdA + ^A'^ + J ij^D.^J, (A.l) 



where 
We have 



D^iP^d^iP-iA'^T'^i;. 



ei23 — e — i, 



Y = a\ (z = 1...3) (A.2) 

where cr* are the ordinary Pauli matrices. Note then that all 7'* are Hermitian. This 
implies that 'ipip and 'ip^'^il) are real, while J dx'ip^^di^ip is imaginary. The gauge field 
will be taken to be A"'Ta where Ta is the fundamental generator normalized so that 
TrT^ = |. Note that 

= b'j;- (a.3) 
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A.l Lightcone gauge 



Let us define 



± ix^ 



V2 ' 

y/2 ' (A.4) 



pI=pI + pI = ^p-^p-. 

We define the lightcone gauge in Euchdean signature by the condition A_ = 0. This 
can be obtained from Wick rotation of the standard Ughtcone gauge in Lorentzian 
signature. However, in this paper we often think of x^jx"^ as purely spatial coordi- 
nates. In particular, in the finite temperature calculation, the thermal time direction 
is orthogonal to the complex lightcone direction. 

Note that under a rotation in the 12 plane A_ — )■ e*"A_. Consequently rotations 
in the 12 plane commute with the condition A^ — 0. 

Defining the momentum space fields 

r{p) = {M-p)r, 

the momentum space action in lightcone gauge is 

/ -^'^M-P)P-A4p) + J -^^M-P) irP, + M,are)m 

r d^v r (fa - ^^'^^ 



where A = A'^Ta. 
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It follows from this action thalP^ 



k 



(A.7) 



(HpU^i-q)) = X i2nY5ip - q), 

Atti 1 

{Al{p)A\{-q)) = -{Alip)A^,i-q)) = - — - x (27r)='5(p - g)5'^^ 
Adopting the notation 

{A;{p)Al{-q)) = {2nf6ip-q)G,,{p)6'^\ 

we have 

G+3(P) = -Ga+lp) = (A.8) 
In coordinate space, the gauge propagator reads 

k J [zirj-^Pg+te k J (zvrj^pi + ^e 



(A.9) 



A. 2 Temporal gauge 



The temporal gauge is defined by the condition ^3 = (Wick rotating the Lorentzian 
temporal gauge Aq = 0). In this gauge, the gauge field propagator is written in position 
space as 



{A{x)A,m 



2m 



eijsign(x )6 (f). 



and in momentum space. 



{A,ip)A,i-q)) 



271 

4:71 



1 1 

+ 



P'^ + ie p'^ — ie 
,3 



{27Ty6\p-q) 



p- 



k (p3)2 + e2 



(27r)V(p-g). 



^^The propagator for a theory whose EucUdean action is 

1 f (fp 



S ■ 



M-p)Q^\p)Mp) 



is given by 



{<^a{p)U-q)) = {2nfS{p-q)Q-,\p). 



(A.IO) 



(A.ll) 



This rule is correct for both bosons as well as fermions. In the case of the gauge field we have 
— — n+3 — -tkp- ^j^g ^^gg ^£ fermionic field we have Q*^^ = ipfj,j^ + M^are and 



47T 

■are • 
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A. 3 Feynman gauge 



We may add to the action a covariant gauge fixing term of tlie form 

SF = ^j d^x^Tiid^A'^f. (A.12) 

Tlie Feynman gauge is obtained in tlie limit ^ — )• oo, in whicli case tlie propagator for 
An becomes simply 



{AUp)Ki-q)) = -y'^'^Vp^ i2nmp - q). (A.13) 



B Perturbative Analysis of Fermion Self Energy 

In this Appendix we give some details on the perturbative computation of the fermion 
self-energy in the three different gauge choices defined in|Xj the Feynman gauge, tempo- 
ral gauge, and lightcone gauge. The Feynman gauge has the advantage of being Lorentz 
covariant, but the cubic Chern-Simons interaction makes computations beyond one- 
loop complicated. The cubic coupling of the gauge field disappears in temporal and 
light cone gauges. The temporal gauge has the advantage of having a very simple form 
of the gauge field propagator, which makes explicit perturbative contributions easy. On 
the other hand, we will see below that it introduces unphysical logarithmic divergences 
in the fermion self energy at two loops (which should disappear in gauge invariant cor- 
relators). The light cone gauge does not suffer from this log divergence problem. While 
explicit perturbative calculations of correlation functions of gauge invariant operators 
are sometimes more conveniently done in temporal gauge, the lightcone gauge allowed 
us to partially solve the planar limit of the Chern-Simons-fermion theory (to all loop 
order) . 



B.l Feynman gauge 

In Feynman gauge, the 1-loop free energy is 

V/ A o \ f -z(|^ + ^) ^ 

f d^q qP{p + qY ^ 
J {2TxYq^[p + qf 

.47rA I g • + g) ^ OttA ^ {p + qf + q^-p^ 



{2'Kfq^{p + qy J (27r)3 q^{p + qy 
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The 1-loop corrected fermion propagator {ipiaip)'4^Bip)) is then proportional to 60,^3/ 




B.2 Temporal gauge 

In the temporal gauge A3 = 0, the one-loop fermion self energy is 



E(p) = 2n\ 



d^'q ■ — q^ 



3 3 

p — q 



(27r)^ 



■00 



dQQ I dq 



p — q 



(p3 _ ^3)2 + ^2 (^3)2 + Q2 



(B.2) 



(p3)2 + Q2 

(a-||p'I 



After subtracting off the linear divergence by tuning the bare mass of the fermion, we 
end up with the renormalized 1-loop self energy 



Si(p) = -|a|/|. 



(B.3) 



We can carry on to the two-loop self energy in temporal gauge, which only receives the 
contribution from the rainbow diagram. It is given by 



p3 _|_ g,3 



4 

= iX^^' 



r 



»^ + q^y + e 
p3 _^ g3 



Si(p)7^ 



dQQ 



{271 f (p3 + g3)2 + g2 q2 

„3 In Ml! 
P ni 02 



(B.4) 



12 



in A 



If 



We believe that the log^ divergence is an artifact of the temporal gauge and should 
drop out in gauge invariant observables. A systematic treatment of regularization in the 
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temporal gauge appears to be rather complicated. In the one and two loop corrections 
to three point functions of gauge invariant operators that we will perform later, such 
divergences are avoided by special choices of position and polarization configurations 
of the operator insertions. 




B.3 Light cone gauge 

The 1-loop fermion self energy in the light cone gauge is (note the Gamma/Pauli matrix 
identities = 1 + 7^, = 1 — 7^) 

''''' J (27r)3 y + ie^ ^ + ie^ ) {p - q)1 + ie 



= SttA 



(27r)3 52 + ie(p- q)l + ie 



We will now use dimensional reduction to regularize the integral, which will turn out 
to be convenient at two loops and higher. The integral can be done by separating out 
two lightcone directions q'^ , q~, and the remaining 1 — e dimension. After performing 
the 1 — e dimensional integral, wc obtain (constant factors of the form 1 + 0{e) are 
ignored, as there is no logarithmic divergence here) 

47rA ' "^'^^ ^^"^^ 



(27r)2|g,|i+^(p-5)2 + ie 
-27rAfl + 6) C dxil-xY'^+i I ^'g^ (g + xp)+(g-(l-x)p)- 

-7rA(l + e) / dx(l-x) 2 + 2 / 7-^ — ^ ^ ' , , ^ ^ 

A/" , \— 1 111/ \- 
dx [1 — X) 2 — 3|PsF2 (1 — a;) - 



2 

= Mps\- 

So with dimensional reduction, the bare fermion mass is zero, and we find the renor- 
malized 1-loop self energy 

Ei(p) = A|p,|. (B.7) 
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It is useful to write the 1-loop correction to the fermion propagator in position space, 

_A/^Ji!lLe'.-.-^^ (B.8) 



(27r)^ + ie An \x\ 

The only potentially nontrivial planar contribution to the 2-loop fermion self energy 
comes from the rainbow diagram. It is given by 



2 + / (Pqs 1 ijp-q) _ .,2 - 



where the integral in the last step is performed as in section 2.1.1. From this we find 
the two loop contribution to the fermion self energy, 

S2(p) = -tX'p+l'-. (B.IO) 

As shown in section 2, the planar fermion self energy is in fact two-loop exact in the 
dimensional reduction scheme. 



C Lorentz invariance of the Wilson line at one loop 

30 

The exact fermion propagator computed earlier in this paper in lightcone gauge 
is not Lorentz invariant. As the fermion propagator by is not gauge invariant, it is 
not physically observable. In this section we will compute a gauge invariant physical 
observable that is closely related to the fermion propagator, and demonstrate that this 
observable is Lorentz invariant (at one loop), and moreover takes the same value in 
lightcone gauge as in Feynman gauge. Our results may be taken as evidence that the 
lightcone gauge employed in this paper defines a Lorentz invariant theory, atleast at 
first order in A. 

Consider the Wilson line 

W = l^{0\[Pe'^^y'^-'f^^{x,)n^P{xfnO) (C.l) 

where m,n are colour indices, we do not indicate the spinor indices of '4){xi),il){x f) 
explicitly (the Wilson line is a matrix in spinor space, in the same way that the prop- 
agator of the previous section is). The expectation value of the Wilson line operator 
described above is Lorentz invariant if and only if it takes the form 

{xf - Xi)^'^^A{\xf - Xi\) + / B{\xf - Xi\). 
SQrpj^jg appendix was worked out in collaboration with Sayantani Bhattacharyya. 
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In this subsection we will demonstrate, at the one loop level, that the expectation value 
of the Wilson loop in lightcone gauge really does take this form (at one loop we find 
that A vanishes while B is nonzero). 

At leading order, the Wilson line defined above is simply given by 

d^P _^i{xi-xf).p_ 1 



and is manifestly rotationally invariant. 

At next to leading order (or leading non-trivial order in A) the Wilson loop receives 
contributions from two graphs. The first arises from the leading self energy correction 
to the ipijj propagator. In the lightcone gauge it is given by 

As we have seen above this diagram is easily evaluated, however we will find it useful 
to leave the integral unevaluated at this stage. 

The second graph has one insertion of the "Wilson line factor" A ■ dl along the 
y direction (direction from Xi to Xf) ; this graph evaluates to 

W, = — f -f^^ie^P-i^^-^f) - e^p'-i^'-^f)) P"7p7V^7. ^ , _ . . 
2 J (27r)3(27r)3^^ ^ - p') ■ ^ ^ ^ ^ 

(this result is correct both in lightcone gauge and in the Lorentz invariant gauge we will 
employ below). In lightcone gauge, it is possible to show that this diagram evaluates 
to 

d'^q d^p e'^^'^^^~^f^ q.y 



^""^'(271)3(271)3 p2^2 ^(p-q)+ {p-q)-y^ 



In order to obtain (C.4) we have assumed that our regulator respect the 3 to —3 flip 
symmetry as well as rotational invariance in the 12 plane. These properties are both 
obviously true of the dimensional regulator employed in this paper. 

It is obvious that neither Wi nor W2 is rotationally invariant by itself. However 

W = W,^ = 4.A / ^J ^'-"-:'\ (C.5) 

J (27r)3 (27r)3 p^q^ [p - q) . y 

Note that a rotation of x/ — Xj can be undone by a rotation of the dummy variables 
p and q establishing that the Wilson line depends only on |x/ — establishing rota- 
tional invariance in any regularization scheme - like the dimensional reduction scheme 
employed in this paper - that respects rotational invariance. 
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It is interesting that the calculation of the Wilson loop proceeds along very similar 
lines and yields exactly the same answer in Feynman gauge. In this gauge, once again, 
we receive contributions from two graphs. The self energy graph gives 

(Pq (Pp e*P'(^'~^^) Q ■ (p — q) 



W[ = (47rA) 



(27r)3 (27r) 



p2g2 



ijp - if 



(C.6) 



The second graph is given by (C.3) in which we must use the gauge propagator appro- 
priate to the Feynman gauge, and evaluates to 

d^q (Pp e'^-i^^-^f) q ■ {p — q) 



-4nX 



q-y 



(27r)'^ (27r)3 p^g^ (^p ^ qy (jj^q^ y 
Clearly W[ + W!^ = Wi + W2 = W where W is listed in dOSj ). 



) 



(C.7) 



D Unitary Representations of the d = ?> Conformal 
Group 

Unitary representations of the conformal group are labeled by the spin s and a scaling 
dimension A of their primary states. When s > 1 these labels are subject to the 
inequality A > s + In this case representations that saturate the inequality are short; 
the null states in this representation fall into a (long) representation with A = s + 2 
and spin =s — 1. 

In the special case s = | it turns out that A > 1. The representation with A = 1 is 
short and its null states fall into a (long) representation with A = 2 and spin =^. The 



later is the representation of a free fermionic field, and the character (3.1) -Fp(x, /i) of 
this field is given by 

(1 — /iX)(l — /i ^x) 

Finally, when s = we have A = or A > ^. The representation with A = has 
a single state. The representation with A = ^ has null states in a representation with 
A = I and s = 0. The states of a free scalar field fall into this representation; the 
character Fs{x,fi) of this representation is given by 

We now present character formulae for all unitary representations of the 3d confor- 
mal algebra. Let us define 

s 

j=-s 
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XsifJ') is, of course, the SU{2) character at spin s (here s is a positive integer or half 
integer). Let us also define 

Ga,s is the partition function over states that are obtained by acting on an SU{2) 
primary of spin s with an arbitrary number of derivatives, and so yields the character of 
any long representation of the conformal algebra (i.e. any representation with A > s+ 1 
for s>lorA>lfors = ^orA>| for s = 0). 

Characters xi^^fJ') of short representations of the conformal algebra are obtained 
by evaluating the character of a hypothetical long representation of that algebra and 
then subtracting out the character of its null states. It follows that, for s > 1 



For s = ^ 



while for s = 



Xs+i,s(a;, fi) = Gs+i,s{x, /i) - Gs+2,s~i{x, fx) (D.5) 

{x, f^) = G^^i{x, fi) - G^^iix, fi) = Ff{x, fi) (D.6) 
Xioi^, /^) = ^10(0;, fi) - Gs o(x, fi) = Fs{x, /i) (D.7) 



E Primary Operators in Free and Interacting Fermion 
Theories 

In this appendix we present the details of some slightly tedious computations involving 
free and interacting fermions. 

E.l The generating function of conserved currents for free 
fermions 

As we have explained in section [3] above, the single-trace primary operators of the free 
fermion theory satisfy the condition that their scaling dimension, e and spin s satisfy 
the relation e = s + 1, with the exception of ■?/'■?/', which has scaling dimension 2 and 
spin 0. Here, we will determine explicit expressions for the corresponding primary 
operators. 

How can we produce an operator of spin s and dimension s + 1 in a theory of free 
fermions in c? = 3? We are interested in operators built out of fermion bilinears (with 
colour indices contracted). As ip and ip each have unit scaling dimension, the operator 
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of interest must contain exactly s — 1 derivatives. As s — 1 derivatives can give rise to at 
most s — 1 free traceless vector indices the remaining index (to make our operator spin 
s) must come from a 7 matrix. |^ Consequently if we define the generating function 
F such that 

0{x; e) = ^F(7, 5^, K = 4:L..,s^'' 
then F must take the form 

F = ^.ef{t,K^ (E.2) 



Following 31 we denote the arguments of / as vectors u and v, so that u = d and 
V = . The fermion equation of motion gives 

M.'y = {7.7 = = f ^ = (E.3) 

It is convenient to change variables to y = u — v and z = u + v. Then we have: 

z-j = y- j = y- z = 0, —if = = 2u ■ v (E.4) 

Terms in / will be of the form: 

f = A + Be-y + Ce- z + Du-ve-e + ... (E.5) 

where each coefficient is a number. 

If we define w = {u ■ v)e ■ e, z = e ■ z, and y = e- y, then / can be thought of as a 
function of three variables f{z,y,w). 

The condition that each current is conserved can be expressed as: 

iK + dl)^F = (E.6) 

which translates into the following condition on /: 

zd^f + dj = (E.7) 

The condition that each current is traceless can be expressed as: 

^ ^ F = Q (E.8) 



del" de^ 



As 7(;i7iy) — rj^iu we cannot have more than one of the current indices come from a 7 matrix. Also, 
the equations of motion teh us that a 7 matrix contracted with a derivative vanishes in its action 
on a fermion. The contraction of a derivative with a 7 matrix, sandwiched by other 7 matrices, can 
also be reduced to a form with a single 7 matrix and derivatives using the equations of motion, (e.g. 
lip.lp^^'^lv) = ~7(/it^''7i/)7p + l{nd'^'nv)p where the first term on the RHS vanishes by the equations of 
motion and the last term has a single 7 matrix). 
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which translates into: 

(59^ + 2wdl - d^y + dl + 2zd,d^ + 2ydyd^) / = (E.9) 



X 1 2 2 

The solution to (E.7) is f{w,y,z) = g{y,w — y). If we define t = w — y, the 



equation for g{y,t) is 

{4d, + 2td^ - dl + lydyd^) g{y, t) = (E.IO) 
The general solution satisfying 5f = latt = y = 0is: 



2k^/2t + y^ 

where k is any constant, which we take to be 1/2. 
The final form for / is thus 

^ . exp (u ■ e — V ■ e) sinh \/2u ■ ve- e — Au ■ ev ■ e 

fi.u, v,e) = ^ ^ ^ ^^-^2) 

yj 2u ■ ve ■ e — Au ■ ev ■ e 

Expanding the above expression in a power series around e, we obtain: 

f = l + e(u-v) + -e^ (3^2 - lOuv + 3v^ + 2w) 
6 

+^3 ^ ^ ] 



6 6 6 3 6 3 ; 

+ Y^e^ (10(m - vf{2w - Auv) + {Auv - 2wf + 5{u - v)^) 

(above, w = u ■ v, u = u ■ e, v = v ■ e.) which yields the currents reported in (|3.9[). 



E.2 Two-point functions of primary operators in the free the- 
ory 

In this subsection we explicitly compute the two-point function of conserved currents 

{O{x-e,)O{0;e2)) (E.13) 

determined in the previous subsection, and demonstrate that the two point functions 
of currents of different spin are orthonormal. We take the two-point function of the 
basic fermionic fields to be given by 
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E.2.1 Simple examples 

To set up notation and get intuition we first work out some simple examples. 
As a first example, consider the two-point function of two scalar currents: 

Using Wick's theorem we rewrite this as: 

Ti{-i){0)ip{x)){i){x)ip{0)) 
where the trace is over gamma matrix indices. We then compute it to be: 



(E.15) 



(E.16) 



a 



2 ^y^t^ n au 

V"/- 3,6 



- n 



Tr7^7^ 



To evaluate more complicated two-point functions, we make use of the identity 



(E.17) 



Note that (E.17) is symmetric under interchange of p and a. 



We next consider the two-point function of two spin-1 currents: 

{ij{x)-f''ij{x)ij{y)-f^'ij{y)) 

We have: 

{^.(rc)y'V'(rc)'i/:'(»)7'''/'(»)> =Tr{-^.(9)V^(a:))7''{Vj(rc)'i/:'(!/))7'' 



(E.18) 



2x'^x'' 



x" 



^p,u 



X'' 



As above, we will often set ?/ = in the last line. 



E.2.2 Results for all spins 

Let e be a null polarization vector. The two point function of the generating operator 

be 



0{x\ s) = ^ jjii--p.s^'^^ ■ ■ ■ is evaluated in the theory of free complex fermions to 



{O{x;e)O{0;e)) 



N 



327r2x2 



Ae ■ X 



x^ 



(E.19) 
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Expanding this in £, we have 



(E.20) 



where J''^\x; e) is the spin-s part of 0{x; e). The spin case is special, where we have 

(jW(a;)j(°)(0)) = ^|a;|-l 

We have defined the set of currents j^^^ with a different normahzation convention, 
namely normalizing the norm of the corresponding state in radial quantization. The 
relative normalization between J^^^ and j^^^ can be determined as follows. If we define 
e) = j\^l..^X^)e^^ then 



(/-)(a;;e)/-)(0;e))=2- jyj;^ , (E.21) 



for s > 0. In the spin case, we have {j^^\x)j^'^\0)) = \x\-^. Prom this we deduce 



2' 



,327ri s\ 



2 



(E.22) 



/2N 
4-) 

theory, receives quantum corrections. 



for s > 0. In the spin case, Qq = In the interacting Chern-Simons-fermion 



E.3 Explicit computation of the divergence of J^^^ 

In carrying out all our manipulations below, we use the fermion equation of motion 

D^j^iP = D^iPj^" = 0. (E.23) 

Some useful identities are: 

YD^YDui' = 0, 



(E.24) 
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Note our convention is such that [^^, '^u]'^ — —iFn^i/j. We now use the equation of 
motion for F^,, = F^i.T"', namely 

eM-PF,, = -J^ or f;, = -e^^.J'^", (E.25) 
where = J^T"-, = ■ipji^T'^ip. It is also useful to have 

^-\hpi^J^'^+ripi^4'^), (E.26) 

and 

^J;T« = -- (j(o)V^7p + Ji'^H^') ■ (E.27) 
To derive these relations, we used (r«)}(r")J„ = and the 3d Fierz identity 

We can now proceed to explicitly compute d^j'^Xi^i- First, consider the current J^^^ 
which is not traceless, 

Using the identities above we find that (before subtracting the trace) the divergence is 
given by: 

+ iPY' (6(fe'fe.2 + ^V2^') - 10(fe'5,, + fe.,^')) ^ (E.28) 
We now substitute for D"^ and further simphfy: 

e^'^^S., = -#7^ (32^,,F,,^ + 32F^,,:^,, + 2r],,,, (pxF^^ + F^^t^)) 

- i^Y' {l^i^'^F,,^ + F^.,:^")) t/j-i^ (2^u,F,, + 2F,,:^,, - 6(F',,F,,)) V, 

(E.29) 

where F^ = e^^pF'^P. 
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Now substituting in F and using Fierz identity, we have 

327r 



k 



^F + F 



2tt 

T 
Stt r 



- ill) (2D 



Gtt 

T 



(^7^2^)^1.1 (^^) + 5r^i(^7i^2^)^^ - e;^2AM(V''-C''^i7''^/')(^7V) 

(E.30) 



We now use the identity 



•(-4 



•(-4 



e.2AMV'(- .,7'^ + D '^7,JV'(V'7» = -77.i.29^(V'V')(V'7aV') + 9,,(V'V)V'7.2V' (E.31) 
to obtain the following total for the above sum: 

^d^j'iXu, = y [-977.,.29m(V5V')(V57''V') + 30a.,(V;V')V57.2V' - 189.,(V;7.2V')V5V'] ■ 

(E.32) 



Subtracting the trace we obtain: 



d'J^X., = d^J^X., - ^ [v.,.,{i'l'^)d^{m + 2a.,(V^7.2V^)V^V^ + 2(V^7.2V')5.i(V'V^)] • 

(E.33) 

The indices on the RHS are understood to be symmetrized. 



E.4 Higher spin examples 

While we do not have a general explicit formula expressing the divergence of the current 
in terms of double and triple trace operators, a straightforward but tedious computation 
was carried out using Mathematica to determine the double trace terms. We find that 
d'^jjfili-Hs-i involves the product of two currents of spins Si and S2, with si + S2 = 
s (mod 2) and S1 + S2 < s — 2. It is convenient to write the results in spinorial notation. 
For instance, the anomalous current conservation relation for J^^^ takes the form (up 
to an overall normalization factor) 



(E.34) 
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For spins s = 4 and 5, we find: 



+59, 



(E.35) 



J(0)j(2) 

a3a4Q;5Q;6 



^9(aia2'^a3a4a5a6) 



and 



7(2) 



— ifin,^ f) 7^^) 7(2) _ AOr) 7(i)/3 7(2) 

_Lfi4 7(i)/3^« r? 7^^^ _iSfi7W/3-9 ^ 7(2) +o^fi, F) f) 7(0)7(1) 

^ "^"'(ai ^/3a2 ^0:30:4 -^as-as) (ai ^a2a3^a!4a5'^a6a7a8)^ ^ ^^^(aio:2^a30!4^a5a6'^ "^aras) 

(E.36) 

In writing the above expressions, we have omitted an overall normalization constant 
as well as possible triple trace terms. 



F Perturbative Analysis of Anomalous Dimensions 

F.l Vanishing anomalous dimension of the scalar primary at 
planar level 

In this subsection, wc check explicitly the vanishing of the anomalous dimension of 
the scalar operator ipip at two-loop planar level, as previously argued based on the 
operator identities relating divergences of currents to double trace operators. The 
simplest way to do this computation is in the lightcone gauge, where there is no cubic 
gauge interaction vertex, and the planar two-loop fermion self-energy is finite. There 
are two diagrams that contribute to the logarithmic divergence. 
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The first one (which comes with an additional factor of 2) is computed as 
_4;rA^ f \Ps\{7^H^ - l^^l^) P~ ^ ^^^2 f d^P \Ps 



(27r)3 (p2)2 p2 J (27r)3(p2)2 



= 47rA 



, i?(i,-|) lnA 
5(2, -i) 27r2 

= -lnA. 
2 

The second diagram (the rainbow correction to the vertex) is computed as 

d^pd^q 2p~g~(7+|^7+7'^|^7^ — 7^|^7"'"7^^7"'") 



- (47rA)2 I 
= (47rAf| 
= (47rA)2 dx J 



= -SttA 



(27r)6 (p + g)2(p2)2p^^2 

d^pd^q 4:p^q~ 

(27r)6 (p + g)2(p2)2^2 

d^pd^q 4p'^{q~—xp~) 

^-^y (27r)6 (p2)2[^2+^2^(l_^)+^2^]2 

2 / C^^P bsl 



/ 



(27r)3 (p2)2 

It indeed cancels twice of the previous diagram. 



(F.2) 



F.2 Anomalous dimension of currents at order 1/N 

It has been argued that the higher spin currents Js{x) as well as the single trace scalar 
operator 'tpip receives no anomalous dimension in the infinite N limit. They do acquire 
anomalous dimension at subleading orders in For the currents of nonzero spins, 

this can be understood in terms of the mixing of the divergence of the current with 
double trace operators. In this subsection, we perform a direct two-loop computation 
of the anomalous dimension of these currents at order X'^/N. There arc three diagrams 
that contribute, all of which involve a matter-loop-corrected gauge field propagator 
(which is suppressed by 1/iV because the matter is in the fundamental representation). 
These diagrams are hsted below: 
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(a) ib) (c) 

In the case of the scalar operator ipilj, only the first two diagrams (a), (6) contribute. 

To begin, let us compute the matter loop correction to the gauge field self energy. 
This is given by 



d^q Tr(7 



+ 



2' ' 
32g 



(27r)3 + 



(F.3) 



The factor | is due to our normalization convention for the gauge group generators. The 
contribution to the fermion self energy by the matter-loop-corrected gauge propagator, 
involved in diagram (a), is given by 



(2^ {p + qf 



A 



Z}^,(g)Sg^(g)D.,(g) = z-ai|^ In A + finite, (F.4) 



where D^y{q) is the classical gauge propagator. Define the matter loop contribution to 
the gauge propagator as 

G'^M=D,M)^'^il)D.u{q). (F.5) 
Let us consider for now the anomalous dimension of ipil). The diagram (6) is computed 



as 



[ d\ D,,iq)%-iq)D^^iq) _ N [ , . 

J (27r)3 g2 2 J (27r)3g2^'^ 



as In A (F.6) 



The anomalous dimension is then given by 



Ki^i^) = + 02). 



(F.7) 



The details of the computation of ai, a2 is given in F.3 (using Feynman gauge, temporal 
gauge, and lightcone gauge respectively). The result is ai + 02 = |. 
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F.3 Two-loop non-planar anomalous dimension of the scalar 
primary 

F.3.1 Feynman gauge 

In Feynman gauge, we have 



47r qP 

and so 



Its contribution to the fermion self energy is 

Nti' f d'q q^l^{i> + i)l,-i{i> + i)i 



)3 q'(p + qr 

1 -2g • (p - 

iVvr^ /■ ci^g — q^ 



^ Ntt' f d'q -2q-{p + q)i 

Ak^ J (27r)3 q^{p + qY ^ ' ' 



4A;2 J (27r)3g3(p + g)2> 
The logarithmic divergent part is 
Nti^ f (fq i 



— i- 



4A;2 J {2nYq{p + qy 



- -i^^— C dx a - x)-^ I ^-^^ (F 11) 



-|^ln A. 



The diagram {b) is computed as 

Ntt^ [ d^q 2 N 



4A;2 J (27r)3g3 4A;2 

Therefore we find 



InA. (F.12) 



12 4 3 

F.3. 2 Temporal gauge 
In temporal gauge, we have 



ai + a2^^ + - = \. (F.13) 
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and 

Its contribution to the fermion self energy is 



^^^^ Uk^ J (27r)3 ' (p + g)2g[(g3)2+'e2] 



4A;2 y (27r)3 + g)2^[(g3)2 + ^2] 

. Ntt^ r d^q {{p + qf - (p3 + ga)^ - p^)^^ - (g^ - + (p + g)V) - 2(g3)^(p3 + q^)-f^ 

* 4A;2 y (27r)3 (p + g)2^[(^3)2 + ^21 

Nn' f d'q r I ^^_(p + g)0^3 

4A;2 7 (27r)3 |_g((g3)2 + ^2) + (p + g)2g 

i ips + qsY+p^s _ \q\ ^s + ip + (i)^r 



(F.16) 

The first term gives zero. The integrals of the remaining three terms give logarithmic 
divergences 

Adding these together, we find 

A^ln A 

E'(p)=i^^j^ + finite. (F.18) 

The diagram (6) is computed as 

Nn^ f dMq^r + q^ 1 ^ ^_ ^^^^^^ 



4A;2 7 (27r)3 q^ (g3)2 + ^2 
The result agrees with the Feynman gauge computation. 



F.3.3 Lightcone gauge 
In lightcone gauge, we have 



D^W, = -DM = f-f^^ (F.20) 



and 

'GUq) G',4q)\_2n'' {q'f f-iq^f q'q^ 



GVo(?) G'+4^)J q{ql + iey \ q'q+ q, 



(F.21) 
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Its contribution to the fermion self energy is 



E'(p) = -i 



X 



A;2 J (27r)3' 



_.7V7r2 r d^q _ q+]^ + ^ _ 2(p + + 2g^[(p + g)+7=^ + (p + g)°7^] + 4g~(p + g)+7"^ 

(F.22) 



With a httle manipulation, this integral can be written as 

cP'q 



(27r)3 



{i-xy 



(1— x)"^ (1 — x)^ — 2p^7^ 



—I 



+ 



(6^J*-;^fV)lnA 

(1 — x)~i 4g~7+ 



(F.23) 



In 



A 



4A;2 V67r2 



1 



1 2 

^P^7^ ^P"7^ ) In A 



While this expression is not Lorentz invariant, it is not a gauge invariant quantity. 
Inserting it into a two-point function computation gives 



1 1 2 
6 ^ 3 ^ 3 



5 

24' 



(F.24) 



The diagram [h) is computed as 

ATtt^ r d^q {q+q-y 



F J {2T:fq^{ql + ief 



N 



In A. 



(F.25) 



We have 



5 1 1 

"^ + "^=24 + 8 = 3- 



(F.26) 



Once again, this agrees with the result we found in Feynman and temporal gauge. 
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G Perturbative Analysis of Three Point Functions 



G.l Vanishing one- loop contribution to the three point func- 
tions of scalar operators 

In this subsection, we consider the three point function of the scalar operator Jo{x) = 
■^■ilj{x)ip{x). The normahzation is such that its two-point function does not scale with 
N. By conformal symmetry, 

(MxMxMx,)) = -^ ^^^yf^J^^ (G.l) 



Parity requires C(A) to be an odd function of A. So if the three point function of Jq is 
nonzero, it can only receive contribution from odd loop order. A priori, for instance, 
there could be a one-loop contribution of order 0{\) to {Jo Jo Jo)- We will perform this 
computation explicitly and see that it is in fact zero. 

It suffices to take the limit Xi2 while keeping xis finite. In this limit, we expect 

{Jo{x,)Jo{x2)Jo{xs)) ^ -^3^- (G.2) 



There are four distinct diagrams to consider (as shown in Fig. 3): 

(Al) 1-loop self-energy correction to the fermion propagator {ilj{xi)il!{x2)). 

(A2) 1-loop self-energy correction to the fermion propagator {iJj{xi)i!j{x^,)) (or {ilj{x2)'ip{xs))). 
Such diagrams do not contribute to terms like x]^2Xi3, and will henceforth be ignored. 

(Bl) a gauge propagator correcting the vertex of Jq^x^). 

(B2) a gauge propagator correcting the vertex of Jo{xi) (or Jo{xi)). 

3 3 3 3 







(G.3) 



12 12 12 

Al A2 Bl B2 

Figure 3 

We will work in position space. Let yi, 1/2 be the two ends of the gauge propagator. 
In the limit 2:12 — > 0, the relevant contribution comes from yi, ^2 ~ 0{xi2)- From now 
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we will set Xi = and X2 = x, while taking 0:3 to be large. The contribution from 
diagrams (Al), (Bl) and (B2) (note that (B2) contribution should be multiplied by a 
factor of 2 since there are two diagrams giving identical contributions) are proportional 
to (the overall X3 dependent factor |a;3|~^ is factored out) 



/(Al) = J d%d%{A^{y^)A,{y2))T:T [(^(0)^(?/i))7''(^(?/i)^(?/2))7'^ (^(1/2)^(0;))] 
I(B1) = [ rfVA2(^M(?/i)^^(?/2))TV[7'^(V'(z/i)VS(0))(V'(0)VS(^))(V'(^)VS(z/2))7l , 



I{B2) ^ J £y,d'y2{A,{y^)A,{y2))Tr [r{^{y,)^{0)){^{0)^{y2))Ymy2mx))] • 

(G.4) 

In the following two subsubsections, we demonstrate the computation in Feynman 
gauge and in light cone gauge, respectively. 

G.1.1 Feynman gauge 

The propagators are 

{A,{x)A40)) = -^W^, (^W^(O)) = (G.5) 
The results for the loop integrals are 

Indeed they cancel, as claimed. 
G.1.2 Light cone gauge 

We have seen that the 1-loop correction to fermion propagator in light cone gauge is 
proportional to |a;°|/|a;p. Without loss of generality, we can take the Euclidean time 
direction to be orthogonal to the plane spanned by Xi,X2,X3, and so diagrams (Al), 
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(A2) simply vanish in the hght cone gauge. (Bl) and (B2) are then computed as 



yh \yi\^y2 - x\^\x\^ 



(47r)3 y yf^ \yi\^\y2 - x\^\x\^ 

A 2 [ d^z (fz dt ^ ^^'^ ' - x) + zi ■ Z12X ■ {z2 -x) - Z12-X [zi ■ {z2 -x)+ t^] 



(47r)3 J Zl^ (^2+i2)|((^2_^)2 + ^2)i|^|3 

rdafd^Z2dt^'''^^'''^'-^'^-'^+'''' 



(47r)3 Jo J {azi + t^)l{{z2-xy + t'')2\x\^ 

= 0; 

^(52) = /■^3^3^(?/?2)l/r2l>{(7=^^lM+-7+^lM')(^2-;^)} 



(47r)3 J yf2 \yi\^\y2\^\y2 - 

^ ,2ld'y,d'y2^^y'^^'^'-y'^y''-^y''''^ 



(47r)3 J yf^ \yi\^\y2\^\y2 - 

^ ;2fd-zAdt-^ (.,-., + t^).i,-(..-a:) 



(47r)3 J zj2(zl + t^)Hzl + t^)H{z2-xy + t^)l 

= 0. 

(G.7) 

We see that the contribution from each diagram vanishes individually in the light cone 
gauge. 



G.2 One-loop parity violating structure of {J J J) 

We will work in temporal gauge, which is particularly convenient for this computation 
since the position space gauge field propagator is very simple. With our special choice 
of null polarization vectors £1 = £2 = ^3 = £, all diagrams involving a one-loop fermion 
self-energy vanish (Figure 4(a)). Similarly, the diagrams with a gauge propagator 
coming out of a vertex and ends on an adjacent fermion line also vanish (Figure 4(b)). 
We further set X3 = x = tio, ^ = 7+, Xi = 0, X2 = —S purely spatial, as explained 
earlier. There are three diagrams that a priori contribute to {J J J) in this special case 
(J is assumed to be a current of nonzero spin): 
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Figure 4: diagrams that do not contribute to {J J J) in the special hmit. 




(G.8) 



Figure 5: type of diagrams that contribute to {J J J) in the special limit. 

(I) A gauge propagator coming out of the vertex J{x2) and ends on the fermion line 
between Xi and x^,. This is shown in figure 5(a). Two analogous diagrams, in which 
the gauge field propagator comes out of J{xi) or J{x3), vanish identically because 
one encounters the structure ^^gj^^ = ^tj^^ = in the loop integral, for the special 
configuration of positions and polarizations we are considering here. 

(II) A gauge propagator correcting the vertex 7(0:3). This is shown in figure 5(b), 
where the vertex on the left is J{xs). 

(III) a gauge propagator correcting the vertex J{xi). The diagram is the same as 
in figure 5(b) except that the vertex on the left is J{xi) in this case. Note that the 
analogous diagram with the gauge propagator correcting J{x2) vanishes identically, 
because once again it involves the structure ^^^^^ — ^t^^^ — 0. 

The diagrams I,II,III are exhibited more explicitly in the limit \5\ <S |t| in figure 6. 



3 3 3 




(G.9) 



12 12 12 

I II III 



Figure 6 

G.2.1 {jjT) 

To begin, let us consider the one-loop contribution to {j£{xi)js{x2)T^{x3)). For this 
correlator, diagram I does not contribute as does not appear explicitly in j{x2)- 
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Diagram II is given by 



/ 



A2(^i(yi)^,(y2))IV [7^(V^(t/2)V'(x3))^(£ • V,3)(V^(x3)V'(yi))y(V'(yi)V'(^l)) 

x^{'il^{xi)i){x2))i{'il^{x2)i){y2))\ 

(G.IO) 

At our special configuration xi = 0,^2 = — a^s = x = ico {S purely spatial), ^ = 7"*", 
the loop integral is proportional to 

c?Vc?'2/26(|/?2)5'(fi2)e^,Tr [^^{j, - - ^,h%^+h^{-^, 

1 / 1 ^ 1 



\yi\^\S\^\6 + y2\^ \\x-y2\'^ \x - yi\^ ^ 
= -16i^ y (fzdtidt2€{ti2) [{z'^)%{t - t2) - z+{z + Syti{t - ti)] 

X ^ / ^ V ^ 

(^2 + t2) i + ^)2 + t2] i I [(^ -t2)^ + Z^] I ' [(t - + ^2] i 

In the small 5 limit, the leading contribution to the integral comes from two domains: 

D2:{h-t)M-0{5). 
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The contribution from region Di is easily seen from (G.ll) to vanish at order 0{t ^). 
The contribution from D2 is 



- m 



(fzdt[dt2- 



Stt2 ^ 



1 ^ 

-T O , 



1 



+ [{z+)%t-z+{z + 5yt\t\ 



(t2 + 



d 



(t2 + 
1 



{t'^ + z^y^ 



{t'^ + z^)i 
.5^ 



l + 0{ 



t2 



d?zdt\dt2 



x{[-(.+)^t^-2z+(z + 5)+tf]+3(z+)2t^} 



l + 0{ 



32i 



\6\H^ I 5(1,2) 



[{6 + zy + tp\z\ 



l + 0{ 



6' 



daa^{l — a) / d^zdt2 



[4ti(-a)] 



[22 + a(l_a)52 + a^2]| 



da a 2 / cPzdt2 



-2a2(l - a) 



[z"^ + a(l - a)(52 + atl 



l + 0{ 



5\ 

t2 



-1287r? 



t2 



Diagram III is given by 



rf3yid3y2(A.(yi)Ai(y2))Tr[y(V'(2/i)^(a;i))^(^(a;i)V5(y2))7''(^(2/2)^(x2)) 



X^(V'(X2)V^(X3))^(£ ■ 9 ^3)(V^(x3)^(?/i)) 

The loop integral is evaluated as 

d'y,d^y2e{yl2)5\m2)e^,1i[l%l^{-i2)l\ 



1 ^ 1 

•9 



16z5+ y d^ zdt^dt2e{tx2)\{z^ft\ + + - ti)] 



(G.13) 



(G.14) 



(z2 + t\) I {z^ + tl) ^{z + 6y + tl] i 1 [(t - + ^2] i 



(t2 + 52)i 

(G.15) 
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In the small 5 limit, the integral is again dominated by the contribution from the 
regions Di and D2. The contribution from Di is 



(5+ 

— l&i—^r I (f zdtidt2e{ti2)tiz~^ {z + S)' 



-S{z + 5) 



+ 



= 96i 

= 96i 



{z^ + q)Hz^ + ti)mz + sy + tiY^ 



(fzdU 



f B{2 

in' 



f daa^(l — a) [ d^y 



^5- = 96i _ 



{-a){l-ay 



[y2 + (52a(l-a)]i 



— 2567ri 



5H^ ■ 

The contribution from D2 is 



(G.16) 



-16i 



= 64i 
= 64i 
= Siri 



|t|3 



d^zdt[dt2 



[z+)Hi + z+{z + s)n[it-t[: 



-3z+ 



{Z^ + {t- t;)2) 5 (^2 + ^2) i [(^ + 5)2 + ^2] i (^/2 + ^2)1 



d^zdt\dt2 



d^zdt2{z^f 



[z+)Hl - z+{z + + 3z+(;2 + -3z-^ 

2^1 



zni + z\z + 5y 

{z^ + tl)l[{z + 5Y + tl]l\z\ 



(tf + z 



6^ ■ 



1 + 



^2 



i6 

5+e{t) 1 



3—3- / (iaa2(l — a)2 f d'^zdt2{z~^y 

2' 2) •''0 

-3- / (iaa5(l — a)5 
> 2) -^0 



z+tl + z^{z + 5)+ 



[{l-a)z^ + a{z + Sy + tl]^\z\'^ 



l + 0{ 



I 



X / d:'z 



+ 



3{z+y{z + 5y 



|z|4(z2 + 2a5-z + a52)f ^2(^2 + 2a5-z + a^^)! 



^2 



X 



B(i,2) ' ^ 



62 



Bill) 



[3{zmz+sr] 



= Svri 



[6(^2 + 2a5-z + a6^) + (1 - 6)^2] 5 



52 ■ 



1 1 1 /-i f 65 

^ / daa2(l -a)2 / rf6 <^ — 

2' 2) jQ y 



1)2 



bHl-b). ^,3 



[3(a6)'(l - a6)] 



X / d'z- 



[z^ + ab{l - ab)5^]-2 



52 ■ 



= 1287ri 



l + 0{ 



t2 



(G.17) 
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Putting back in normalization factors associated with the propagators, we obtain 



{JJT) 



l—loop 



27r3 54^7 



(G.18) 



in our hmit. 



G.2.2 (TTT) 

Next, we consider the one-loop contribution to {T^{xi)T^{x2)T^{x3)) . Diagram I is 
given by 



X Tr 



f{i^{y)i^{x,Me ■ 'S',,){i^{xi)i^{x2)W{i^{x2mxsMe ■ ,,){i^{xs)^{y)) 



(G.19) 

At our special configuration xi = 0,X2 = —S,X3 = x = ti-Q {S purely spatial), ^ = 7"*", 
the loop integral is proportional to 

J d'ye{y')5\y- X2)e,,eiT. - ^,)^^^,,^^^,,^^{^, - ^)] 



X 



\y-xi\^ \xu\^ J \\x23\^ ""'Ixs-yl^ 



1 ^ 1 



- I £ye{y')6\y + 6)e,jS^Tr[Y^^ 
9 



X 



y+ 6+ 
+ 



\y\^\S\^\x — yl'^lx + \y'^ S"^ J \{x — y^ {x + S) 



y+ , 5+ 



-9 • m 



dy'- 



1 



1 



1 f 



X 



- y0)2 + 52 (a;0)2 + ^2 



= 9 • 16i 



dy' 



e(j/°)(|/°)3(2x°-i/°) 

+ (1/0)2)1(52 + (3.0 _ ^0)2)1 



(G.20) 

In the small S limit, the integral is dominated by the contribution from (x°— y°) ~ 0{6), 



9 • m 



dy' 



e(xO) 



WW J (52+ (xO -?/0)2)5 

and from y^ ~ 0{5), 



52 ■ 



-192i 



l + 0{ 



52 



9a6^(^/.,o^(^^ = -192.r)^ 



10 



((52 + (yO)2)i 
107 



'\5\H^' 



t2 

(G.21) 



(G.22) 



The former is an artifact of temporal gauge, and will cancel against diagrams II, III. 
The latter is subleading and does not contribute to the overall coefficient of the parity 
odd (TTT) structure. 

Diagram II is given by 

d?y,d?y2{A,{y,)A,{y^))i:i [^^ {ij{y2)i,{x,me ■ ,,){i^{x,),l^{y,)W {i;{y,)i,{x,)) 

^ ^ _ ^ y _ 



(G.23) 



The loop integral is computed as 



d%,Sy^e{yl^)5\y,^)e,,Ti [^^{i, - - i {-j ^ 



\5' + y2\ 



■16i(5+ j (fzdtidt2e{ti2) [{z+)H2{t - 12) - z^{z + (5)+ti(t - ti)] 



[{t-t2f + Z^Y 



[(t-tl)2 + ^2]3 



<5'=<5 



— (1^ ) ^ 

{tl + z^)V- ■ - ''\tl + {z + 5'f]i 

(G.24) 

In the sm all 5 limit, the contribution is again dominated by regions Di and D2 as 
in (G.12). The contribution from D2 is of the form ^'^"^g^g'-*^ 1 + , and should 

cancel in the end just as in the {jjT) case. The contribution from Di is easily seen to 
vanish at order ^p^, and does not contribute to the one- loop coefficient of the parity 
odd (TTT) structure. 

Diagram III is given by 



d'y,d'y2{MyMAy2))T^^ ^ ■ d .,)(^(xi)^(i/2))y (^(l/2)^(x2)) 



^ y _ ^ y _ 



(G.25) 
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The loop integral is computed by 



1 



X 



5 X- 



5 <5- 



■16i5+ y zdtidt2t{ti2)[{z^ftl + + - ti)] 



1 



X 



[(t-tl)2 + z2]i 



(t2 + 



1 



[ii + (^ + 5')^]^J,,_, 

(G.26) 

In the small 5 limit, the contribution from region D2 should cancel against diagrams 
I, II. We only need to compute the contribution from region Di, given by 



-16i(5+ j d^zdtidt2t{ti2)z'^{z + 5)+tit 



1 ^ 

o , 



X 



(t? + ^2)f - - (t2 + ;,2)f 

-3(2 + 5)- 



(t2 + z2)f 



.5^ 



-162— / d^zdtidt2t{ti2)z^{z + 5)+ti 



5+ /" 

27-162— / d^zdtidt2t{ti2) — 
t J {tl 



{t\ + z^)'i{tl + z^)'i 

zm{z + 5y)H,{ti-ti) 



9{{z + 5)+f 
[tl + {z + 5)2] 



Using 



dtit{t 



^ + z^)2{tl + z'')2[tl + {z + 5Y]^ 
hitj-tl) 4 1 



12) 



We can write the RHS of (G.27) in a Lorentz invariant form 

.6^ 



27-16.- /A^^^^^^^^^^ 



-9 - 2567rz 



(G.27) 
(G.28) 

(G.29) 



This is the only nontrivial one-loop contribution to the parity odd (TTT) structure. 
Putting back in the normalization factors, we obtain 



{TTT)i_ 



loop 



(G.30) 



109 



G.3 Two-loop parity preserving structure of the three point 
functions 

The goal of this section is to compute the (parity-preserving) two-loop contribution to 
three point functions {J J J). In particular, we will demonstrate that the "free scalar 
tensor structure" of {J J J) receives a nontrivial contribution at two-loop. 



G.3.1 (TTT) at two-loop 



Let us begin with three point function of the stress-energy tensor (TTT). T^^yi^x) is 
given by 



T 



(G.32) 



(G.31) 

or when contracted with a null polarization vector e, 

The three point function {Ts-^{xi)Ts2{x2)Ti;g{x3)) takes the form 

(Te^(xi)Te2(x2)T£3(x3)) = CBF^22i{^i^^i}) + CF^2^2({3^i' ^i}) + parity violating term. 

(G.33) 

where F^2{{^iy^i}) and F^2{{^ij^i}) are proportional to the three point function of 
stress-energy tensors for a free scalar and a free fermion, (TTT) b and (TTT) p. Our 
normalization convention is 



gi^^(- Qi+Q2+Q3 ^ sinh Pi sinh Pg sinh P3 



Fi2|F23F3i| 

cosh( 2i±^±^) cosh Pi cosh P3 cosh P3 



F12I |a;23F3i| 



(G.34) 



A:Ja|A| 



The parity- violating term in (G.33) has been described in previous subsections; it 
only receives contribution at odd loop order, and its one-loop contribution has been 
evaluated explicitly. 



The constant coefficients cb and cp in (G.33) have the following perturbative ex- 
pansions in A in the planar limit: 



5b = 9oN{cf\^ + ■■■), CF = g^N{l + cfX' + ■■■) 



(2)x2 



(G.35) 



where is an overall constant having to do with the normalization convention on the 

f2) 

stress We will now compute the two loop contribution . Since we are only interested 
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in extracting the coefficient of (TTT) b tensor structure, it suffices to consider a special 
configuration of at which {TTT)f vanishes and {TTT)b simphfies. 

The stress-energy tensors of a free complex scalar (p and a free complex fermion 

are 

Tf = 2\e-d<P\'- ^-{s ■ df\<P\\ Tf = ^p^e ■ V^^. (G.36) 



With this normalization, the three point functions are related to ^-^d (G.34) 
by 

432 ^'^■^^'> 



Now, restrict to the special case ei = €2 = = e. Further, at the special configu- 
ration Xi3 ■ e = (the labeling of the three operators is an arbitrary choice), it can be 
seen that {T^TfT^) vanishes, while (TfTfTf) does not vanish. 

To extract cb, it suffices to take the limit X12 = 5 — ?■ 0, while keeping X13 finite and 
subject to the condition xis -6 = 0. In this limit, we have 

31 ^ . ;r^6 

(Tf (xi)Tf (x2)Tf (X3)) = -^s^^,^ {Tf{x,)Tf{x2)Tf{x,)) = 0. (G.38) 

We find it the simplest to work in temporal gauge Aq = 0, when the CS cubic 
coupling disappears and the gauge field propagator is 

9'TTi 

{A{x)MO)) = —e,jsign{x')5\x). (G.39) 

Without loss of generality, we will choose the temporal direction to be orthogonal to e; 
in other words, e is assumed to be a purely spatial complex null vector. Note however 
we can not assume all of Xi^2,3 to be purely spatial. We can proceed by assuming that 
Xi,X2 are purely spatial, while X3 is a general vector. As we have seen in the two-loop 
fermion self-energy, there can be unphysical divergences in the temporal gauge, which 
are expected to drop out in gauge invariant correlators. 

There are potentially many two-loop diagrams that could contribute, but most of 
them do not contribute to terms of order 

(G.40) 



in the limit 



Ei = e, = e ■ cq = e ■ X13 = 0, 5 = X12 — > 0, X13 finite. (G.41) 
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For instance, any diagram with a bare propagator {iIj{xi)i(j{xs)) vanishes, because this 
propagator is inserted in a trace Tr y{ip(xi)ip{x3))^ • • ■ ] oc x^ -^Tr ■ ■ ■ ] = 0. Shghtly 
less obvious is that diagrams involving bare propagators {ip{x2)'4'{xi)) or {ip{xs)ilj{x2)) , 



thought not identically zero, do not contribute to (G.40) in the limit (G.41) either. 
The former comes with at most two derives acting on e ■ S/\S\^ from the structure 
Tr ^^{ip{x2)ip{xi))^ ■ ■ which at most goes like \S\~'^ in the small S limit. The lat- 
ter involves at least one derivative e ■ dx^ acting on the factor ~ l/|a;i3|^ from the 
denominators in the two fermion propagators connected to X3; in this case there are 
no additional factors that scales like X13 in the numerator, and so the result cannot 
contribute in the large X13 limit at order (9(|a;i3|~^) as in (G.40). Further, if one of the 



fermion propagators between two T vertices, is corrected by a one-loop 



self-energy diagram, such a diagram vanishes in the configuration (G.41) because one 
encounters Tr(^I^---) in the trace over Gamma matrices. Therefore, the only dia- 
grams we need to consider are of the form shown in figure 7. In fact, there is only one 



such diagram that contributes in the limit (G.41), namely the one with T[x^) being 
the vertex on the left, in figure 7. 




Figure 7 



This diagram is computed as 



2 

i=i 



X Tr 



{ip{zi)i){x^))i^£ ■ d x-A'^{x^)il){z2))'^k{ij{z2)i){x2))i^e ■ d X2{i^{x2)i){y2))lie 



ki 



(G.42) 

In the limit ( [a4l| , ( |a42| ) is given by the following loop integral multiplied by a 



normalization factor — ^(27riA)^(47r'' ^ 
3 



Fl3 



d^yi(fy2d^zi(fz2 ^r—j^ 



i=l 



1 ^ 
7T^£ ■ O 



Xi Zi 



Fi - yi\\ 



(G.43) 
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Using the assumption that e is purely spatial, we have 7j^7je*'' = 0. The trace of 
product of Gamma matrices in the integrand can be simplified as as 

= 4Tr - ^,h%f, - - mh - k)] 

where denotes the spatial part Xi ■ ^,yi ■ 7. Define i/i = (— yi), i.e. the vector 

with the time component flipped. Taking into account the restriction = yf due to 
the delta functions in the temporal gauge propagators, we can write the above trace as 



4Tr 



32z^z^e ■ y2ie ■ {xi - yi)e ■ {x2 - ^2). 



(G.45) 



The loop integral (|G.43|) can be then written as 
27-32 



X 



Ylsign{y^ ~ z^)6\yi - z. 



z^ 



(G.46) 



I /} 1 , fY» .151 . />'■ . I 5 

After performing the integral over z^, Z2, which turn out to give very simple results, 
(G.46) reduces to 

3-29 



d-'yid'y2 



\yi2\ 



3-29 r(f) 
|xi3|«r(3)2r(|) Jo 



X / d^yid% 



da 



l-a 



3 

2 1,2 



\yi - x\\^\y2 - X2\^ 



dha%\\- a-hy- 



\e ■ {yi2 + 6)e ■ yie ■ y2X 



[ayl + hyl + (1 - a - h){yi2 + 5^] ' 



(G.47) 



3-29 r( 



15^ 
2 ' 



\xv,\'m'm Jo 



da 



l-a 



dha%\l - a - 6): 



X / d%d% 



{u/a)\u/hf{l - u/a - u/h)\e ■ 6)' 
[ayf + hyl + (1 - a - h)yl2 + ^ 
where in the last line we defined 

ah(l — a — b) 
a(l — a) + 6(1 — h) — ah 
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u 



(G.4J 



We can further shift y2, and calculate the integral (G.47) as 



|xi3|6r(3)2r(|)yo Jo ^ ' 

3 ,3 yu/af{u/hf{l - M/a - u/hf{e ■ 5f 



[ayf + %i + (1 - a - &)z/?2 + 



15 
2 



3-29(e-5)6 r(f) 87r3 /-i /•^-'^ ^ a^^^^l - a - 6)^ (u/a)2(M/6)2(l - ^/a - 



-1 yl-a ..2^21 

, , da db 

ix.slW r(3)2r(|) 1287 70 Jo 



9 

ti2 



a6{l— a— fe) 



3-29(£-(5)6 ri 



15' 
2 ' 



in" 



da 



l-a 



db- 



3 7 3 

a^b^ 



\xu\'\S\^ r(3)2r(|) 1287 y„ J„ [a{l-a)+b{l-b)-abf 
Putting back in the normalization factors, we found 



{Te{xi)Ts{x2)Te{x^))2-loap 



32 



21% |Xi3|6|5|9' 
(Tf (xi)Tf (x2)Tf (0:3) )xi3-e=0,xi2=5^0 



^ irrBf^ \rrB(^ \rrB 



This determines the two-loop coefficient in (|G.35l) 



(G.49) 



(G.50) 



.(2) 



TT 

32' 



(G.51) 



G.3.2 General spins 

Now let us consider the two-loop contribution to the three point function of higher spin 
currents, of the form (J^^ J^^ Jsg). Once again we will take all e^'s equal to £, restrict to 
the case e ■ X13 = 0, and take the limit 0:12 = 5 — )■ 0. As before, the free fermion {J J J) 
tensor structure vanishes in this case, and we can extract the coefficient of the free 



scalar {J J J) tensor structure by taking the limit (G.41). More precisely, we expect 



[Js^{xi]e)Js.,{x2]e)Js^{x^]e)) Cf_^ 



''a^s |^^g|2s3+2|^|2si+2s2 + l 



(G.52) 



in the limit of (G.41), and shall now determine the coefficients C^sj^a' ^ generating 
function for the higher spin currents in CS-fermion theory (with a null polarization 
vector e) is 

J(e) = V^^/(e-^,£-:^)^ (G.53) 
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where 



f{u,v) 



(G.54) 



2^/uV 

Define fn{u,v) to be the part of f{u,v) with total degree n in {u,v). 

As in the discussion for the correlator of stress-energy tensor, the only two-loop 



diagram that contributes to in the limit (G.41) will be the one with two gluon 



lines connecting fermion propagators adjacent to Js^ and Jsj. We can then effectively 
replace by (9^ in the current insertions. The loop integral associated with the 
diagram in figure 7 is 



i=l 



X Tr 



x-l, 



i- 



■{■ip{x2)ij{y2))lee''\^{y2)ij{yi)H{-ip{yi)ij{xi))^fs,-i{e ■ d^^,e- d ^^){ilj{xi)-^{zi))-fje'^ 

(G.55) 

Using the following relation involving the generating function f{u,v), 



\x — z^^ 



\x - z + e ^ ^ (2n + l)\ \ J \ J \x - y 

n=0 



y-£\ 



E 

n=0 



(2n + 1)! 



(2n + 1)! 



{e-{x-z)r {e-{x-y)y 



X — y — e\^\x — z + e\ 



1 + 



X - Z + e|3+2n \x - y - £:|3+2n 

Ae ■ {x — y)e ■ {x — z) 



\x — y — eY'\x — z + e\ 



(G.56) 



we have 
1 



1 



X — z\'^ \x — y\'^ \s — \ J \x — z\'^\x — y\-^ 



and in particular 



e ■ {x — z) e ■ {x ~ y) 
\x - 



s-l 



\x — z\ 



(G.57) 



1 



|X3 - Z2[^ 



\X3 - Zi\^ 



^ 2^3 



-1 



2 ^ |X3|»|^1|,|^2| 



I Zj 

2^,-1 



2 ,(^-(^.-i/.)r^-M(2/°)'-k°)T^"' 

Si - I. 



Inf. . 2s j -|- 1 L-vj . 2sjH-l 

1^2 1^2 i/2 I 



{yi = Zi). 



(G.58) 



115 



The integral (G.55) then reduces to 



^ ' |xi3P^3+2 y ..... . I^^^l 

i=l 



(G.59) 

where 

^s^r-'(~}X (G.60) 



Integrating out zf in (G.59) now gives 



r(Si + |)r(s2+|)|Xi3P^3+2_/ |?/l2nXi-l/i|2^i+2|3,2_^2|2s,+2 

Jo Jo 



r(si + |)r(s2 + ^)|xi3|2^3+2 r(|)siS2 

X / d-'yid-'y2 



,3„. ^3„, (e-(i/i2+5)r^(e-yir^(e-i/2r^ 



K + W2 + (1 - a - + 



X / d^yid^y2 



S1+S2+S3 



[u/ay^{ulhy^{l - u/a - u/bY'je ■ 6) 
[ayl + byi + il-a-b)y!, + 6'ur^'''-''^ 

(G.61) 
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Putting back in the normalization factors, we obtain 



C. 



B 

1S2S3 



X 



N(2m\f 
2(47r)7 



\Si+S2+S3 + l 



?>2'kUs,Us,K, r(gi + S2 + I) 

Tisi + \)Tis2 + \ 



r(|)siS2 



da 



\-a 



NX' 
• 1 

X 



{u/ay^{u/by^{l - u/a - u/bY'^ 
[ayj + byl + il-a-b)y!, + ur^''- 
,^+,^+,3+1 327rW.,Ar,,A/;3 r(gi + g2 + i) 



da 







l-a 



dba'^b'^{l- a-bY 



i{u/aY'{u/by^{l -u/a-u/bY^ 



ab{l—a—b) 



NX' 



\Sl+S2 



327i^Af,,Afs,Afs3 r(si + S2 + 



r(si + i)r(s2 + |) r(|)siS2 



X 



da 



l-a 



rf6 



-(abYHl-u/a-u/bY' 
1 — a — 

1 327rW,X,2A43 r(si + S2 + |) 



2«i+«2+s3r(^5g + 



r(si + i)r(s2 + ^) r(|)siS2S3 
Lr(si + S2 + 1) 



Sl!S2!S3! 



.2 ■ 2' 

(G.62) 

Our normalization for the higher spin currents is such that the two point functions of 
currents for a single free complex scalar and a free complex fermion are given by 



-{Jf{x,e)jf{0,e)) = (Jf (a:,.)jf(0,e)) = ^^^^^ 



1 2^'-^n-"2T{s + l){e-x] 



2s 



\X 



4s+2 



(G.63) 



With this convention, the three point function coefficient of the currents for a free 
complex scalar in the limit of (G.52) is IsT] 



f^free scalar 

S\S2S3 



2-1+.2+.3 r(si + S2 + |)r(s3 + \] 



471^ 



[siy.{s2y.{ss} 



(G.64) 



We see that C*^S2S3 is C^^^^^,^'^"'^"'^ times an Sj-independent constant. This constant has 
been computed in the previous subsection in the spin-2 case. 
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